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PREACH 


In this book I have tried to give a fairly complete 
account of the practical methods of analysing any 
given curve or set of observations into its harmonic 
constituents, if it be a periodic function, and of detect- 
ing any periodic component in it if it is not a periodic 
function. 

The book will, I hope, be readily intelligible, and 
readable, to anyone familiar with the elements of 
trigonometry and the calculus: for this purpose, an 
introductory chapter has been added dealing primarily 
with imaginary quantities, but including also, for 
convenience, a few theorems on the integrals and sums 
of trigonometric expressions that are required later in 
the book. 

I have ventured to use the term “ artificial function ” 
to signify a function which is defined by one mathe- 
matical expression over one part of its range and by 
another expression over another part of its range. 
This may displease some modern writers on pure 
mathematics, who will contend that there is no 
necessity for such a term. But such functions are of 
such wide occurrence in nearly all branches of mathe- 
matical physics that it seems desirable to introduce 
such a term to cover them; and the difference 
between such functions and analytic ones is very 
obvious to all elementary students and _ practical 
workers. Fourier’s theorem, in fact, may be regarded 
as the method par excellence of representing an 
artificial function in the form (or semblance) of an 
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analytic function. And it is partly because of this, 
as well as its property of analysing functions into their 
harmonic constituents, that it is of such great value 
in mathematics. 

The view taken of Fourier’s theorem in this work 
is rather different from that sometimes adopted. 
There are many known ways of expanding an arbitrary 
function in a definite series of analytic functions which 
expansion is valid for all points between two definite 
limits, and Fourier’s theorem may be looked upon as 
merely one of such cases. In the other cases, how- 
ever, we neither know nor care what the series repre- 
sents outside these limits; while the Fourier series, 
outside these limits, represents continual repetitions of 
its behaviour between the limits. We have therefore 
taken Fourier’s theorem as applying only to periodic 
functions, either artificial or analytic. However any 
function behaves over any finite range, we can, by 
repeating this behaviour an indefinite number of 
times to both minus infinity and plus infinity, obtain 
an artificial periodic function from it. 

One advantage in looking upon Fourier’s theorem 
in this way is that the Fourier’s series (apart from 
the constant term) of any artificial periodic function 
can be directly written down from an examination 
of the discontinuities of the function and its succes- 
sive differential coefficients. Moreover, the Fourier’s 
series representing an analytic periodic function can 
be deduced in the same manner as the limit of the 
case when the discontinuities are infinitesimal ; 
although this is rarely the easiest way of obtaining 
the expansion. For instance, the harmonic analysis 
of the function 2cos?¢# can be deduced from the 
analysis of the artificial periodic function of period 2z, 
and represented from o to 2m by 2 cos? gt, by putting 
q= 1. This artificial function possesses, both for 
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itself and for all its differential coefficients, discon- 
tinuities when t= 0; and from the formule on 
p- 59 the Fourier constants are readily found to be 
given by 


pena ids? 
4g — 1 
in2 
and oa EE 
nm — 4g 


Putting g = I gives zero for all values of » save when 
n = 2 for which the expressions become indeterminate. 
If in this case we put g = 1 + 6 where 6 is very small, 
we easily see that a4, = 1 and 4, = o, and thus that 
the Fourier’s series representing 2 cos? ¢ is cos 2¢ plus 
a constant. 

Chapter III deals with the harmonic analysis of 
various simple artificial functions, for the most part 
built up of straight lines, parabolas or sine curves. 
This chapter will afford the student excellent practice 
in the evaluation of trigonometric integrals. I have 
been surprised at the very small proportion of students 
who can evaluate such integrals correctly after having 
been through an ordinary course of the Integral 
Calculus. This chapter accordingly will form a useful 
item in the student’s education. In fact, the whole 
book has been written with the intention of making 
it constitute a valuable chapter in the education of 
any student of applied mathematics, as well as of 
making it valuable to the actual worker engaged in 
any kind of harmonic analysis. 

Chapter V deals with the harmonic analysis of any 
periodic function defined by the co-ordinates of a 
number of points and will be, for many readers, the 
chapter of most practical importance in the book. 
I have tried to include therein all methods that I 
think of practical value, while giving full attention to 
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the most satisfactory ones. Attention may be drawn 
to what is, I think, a new method of harmonic analysis, 
which consists in joining the given points by arcs of 
parabolas of the second or third degree and then 
obtaining the exact harmonic analysis of the artificial 
function defined by this series of arcs from the dis- 
continuities. 

Chapter VI, on theoretical considerations, I have 
added for the sake of completeness, and also because 
I find that most students of physics and engineering 
do show a considerable interest in the proofs of the 
various mathematical methods and operations taught 
them. Very few, in the course of their education, 
simply want to be given mathematical results that they 
can use and betray no interest in the establishment of 
the results. 

‘In Chapter VII I deal with Fourier’s integral 
theorem, the chief application of which is in the 
science of Optics; but since it is the limiting case of 
Fourier’s theorem when the period is made infinite, 
it could not logically be excluded, and it is, more- 
over, wanted for the next chapter. 

Chapter VIII deals with the detection and deter- 
mination of any periodic components in any non- 
periodic function of which observations have been 
made over some finite range of time. The method 
adopted has been to apply the theorem of Chapter VII 
to the artificial function defined to be equal to the 
observations over the region for which these extend, 
and to be zero for all values outside this region. 

The detection of periodic components in pheno- 
mena not themselves strictly periodic is required 
in analysing many meteorological and astronomical 
phenomena. 

Finally, a very brief historical survey of the subject 
has been added in Chapter IX. 


PREFACE 1x 
I need hardly say that I shall be grateful to any 


reader who will be kind enough to point out any 
mistakes or misprints. I shall also be grateful for 
any suggestions, especially from experienced practical 
workers, for incorporation in a further edition of the 
book should one be called for. 
ASE, 
The University, 


Manchester. 
May, 1925. 
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A PRACTICAL TREATISE ON 
FOURIER’S THEOREM AND 
HARMONIC ANALYSLS 


INTRODUCTION 
ON IMAGINARY QUANTITIES 


Brrore proceeding with the proper subject-matter 
of this book, it is desirable, we think, to give a very 
brief account of the so-called, and very badly called, 
“imaginary quantities,” as we shall frequently have 
to make use of them. It is almost impossible, in 
mathematical work, to deal with sines and cosines 
without imaginary quantities inevitably presenting 
themselves ; so it is very important at the outset to 
get clear ideas about them and about their relations 
to the trigonometrical functions. 

Let O and A represent any two points in a plane and 
eevetepresent the: vector’ OA, 1.¢..V represents a 
step of magnitude OA in the direction OA and in the 
sense from O to A. This being so, it is evident that 
2V will represent a step of magnitude equal to twice 
the distance from O to A in the direction from O to 
A; and that in general if a is a positive quantity aV 
will represent a step of magnitude a times OA in the 
same direction and sense. Further, —V will clearly, 
with the ordinary meaning of the minus sign, represent 
something which will cancel or neutralise V, and 
therefore represents a step of magnitude equal to OA 

B 


2 A PRACTICAL TREATISE ON 


in the sense from A to O; and similarly — aV repre- 
sents a step of a times OA in the direction and sense 
of AO. Hence the effect of the factor (— I) is to 
reverse the direction of a vector without altering its 
magnitude. But reversing the direction of a vector 
is equivalent to rotating it through two right angles 
or x radians. Similarly, (— 1) times — V represents 
— V rotated through z radians or V rotated through 
2x radians. Proceeding in this manner, we see that 
if m is a positive integer (— 1)” V represents V turned 
through an angle of ma radians. Hence we may say 
briefly that the factor (— 1)” represents a rotation of 
nm radians. 

Now, exactly as in elementary algebra when dealing 
with indices, we define them first for positive integral 
exponents and then find we must interpret them for 
fractional or negative values because our mathematical 
operations give rise to them, so we must do the same 
here. If(— 1)”+™” is always to represent a rotation 
of (m+n) a radians whenever m+ m is an integer 
(— 1)” must represent a rotation through m2 radians, 
whether 7 is positive or negative, integral or fractional. 

Hence we take (— 1)” as representing a rotation of 
nz, to be universally true. Writing mz = @ or n= 6/z, 
we see that the factor (— 1)®%* applied to any vector 
rotates it through an angle @. In particular, the factor 


(— 1)? or ,/—1, which is generally denoted by 3, 
. MLA . . 
represents a rotation of 5? 2: Totation through a 


right angle. If we denote the vector V turned 


through an angle @ by Ve, we have the result that, in 
general, 


Ve= (— 1% Vi 
and in particular, 


Via aV. 
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We will now obtain another expression for Vo in 
terms of V. Let OA in Fig. 1 denote the vector V 
and OB denote Ve. Draw BC perpendicular to OA, 
then OC denotes the vector V cos 6 or (cos 0).V. 
Similarly, CB is a vector whose magnitude is equal to 
that of (sin @).V, but since the vector CB has been 
rotated through a right angle away from the direction 
of V it is represented as a vector by i sin @.V. But the 
vector OB is equal to the vector OC plus the vector 


4B 


CB, which is merely asserting that the step OB is equal 
to the step OC plus the step CB. Hence we have 
Vere cos 6 ..V — 4 sin 0. V 
or Ve = (cos 8+ 27 sin 6).V. 
But we had previously found that 
Ve= (— 1)9* V, 


hence we see that the two expressions (— 1)9/" and 
cos 9 + 7 sin @ are equal and both represent rotation 
through an angle @. If is any integer and we replace 
6 by n0, we have 
no 
(— 1)* = cos nO + 7 sin 16, 


no 
but (— 1) = {(—1)**}" = {cos 6 + 7 sin 6}. 
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We thus arrive at the result which is known in 
trigonometry as De Moivre’s Theorem, viz. that 


(cos 6 + 2 sin 6)” = cos nO + 1 sin n@. 


De Moivre’s Theorem is generally regarded as belong- 
ing to the advanced part of trigonometry, yet, as we 
have seen, it is nothing more than the expression of 
the fact that a rotation of 7@ is equal to m rotations of 6. 

We can, by De Moivre’s theorem, easily obtain the 
power series for the sine and cosine. If we write 
0 = x/n we get 


eos x Sas n 
cos ¥ + isin x (cos * + isin *) 
n nN 
Now let 7 tend to infinity, while is finite; so that 


cos —-+ 72 sin — represents a rotation through an 
nN n 


. . ° Oe « ° 
infinitesimal angle ~, and thus reduces to 1 + 7~ in 
n n 


the limit. We then have 
cos.* + 7 sin x ={(- + ix)’ ey a) 
n 


when » is infinite. Expanding this by the Binomial 
‘Theorem, we obtain 


cos 4 +-.2sin & = 1 + ax (x — )ita2 + 
n 


[2 
(=D 
3 |4 


when 7 is infinite. 

If we remember that 22 = — 1; 73 = — i, etc., and 
that, when this is done, the terms devoid of the factor i 
represent displacements in the line of V; and that 
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those which contain the factor i represent displace- 
ments at right angles to V, we see that 


cos ¥ = I -(: = y+ — G-Z)G-3 xt 
nN |2 nN nN nN 14 


and 


sine = x—(1—1)(1-2 aio pays 
n nt} \3 


when 7 is infinite. 

Now both these series, when u is very large, are 
such that for any given value of x a finite number of 
terms will give their sums to any desired degree of 
accuracy. Hence we may safely put n= 0 or 1/n=o0 
in these series. For in any term of finite order there 


are only a finite number of the factors (: ~~ ‘) wen 
n 


(: — *), etc., and these give exactly unity when 1/n 
n 


is exactly zero. Hence we have the important results 
that 


x? gt x 
cos 4 = 1 + a+. (2) 
2 4 [6 
and 
3 5 7 
se ies a ee eet tee pene 018) 


B 
It is shown in books on Algebra that the limit, when 
n is infinite, of the expression 


( :. z 
nN 


is ¢2, where ¢ = 2:71821 . . . and is the base of the 
Naperian logarithms, hence we see by (1) that 
COSMAS att ee Cl ee) G- (4) 


This is consistent with the previous value, viz. 
(— 1)*/", found for the left-hand side since it only 


6 A PRACTICAL TREATISE ON 
implies that (— 1) = e™=cos a +7 sin 2. If we 
change @ into — @ in the equation cos @ + 4.6in 0 =e 
we get COS)O =e? SIO een ma 


adding and subtracting these equations we have 


ei8 e—i8 
COs i) == = a 
D 
and 
’ i — p—i0 
sin 0 = 
21 


These results are known as the exponential values 
of the sine and cosine. Their meaning can easily be 


D 


Fic. 2. 


grasped from the diagram, Fig. 2. The factor e°, 
since it is the same as cos @ + 12 sin 0, represents a rota- 
tion through an angle 8. Hence in Fig. 2 OA, OB and 
OC represent respectively V, e° V and 2a" V seClear a 
the sum of the latter two vectors is a vector equal 
to 2 cos @ times V. Similarly OD represents the 
vector — e=* VY and the resultant of OB andsODa- 
clearly a vector of magnitude equal to 2siné@ V, but at 
right angles to it, and hence represented by 2 7 sin 6 V; 
this is equivalent to the equation e? — ¢— = 27 sin @. 


Any expression of the form a + bis called a Complex 
Quantity ; @ is said to be its “real” part and 1, or 
briefly }, is said to be its “imaginary” part. Much 
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better and more appropriate names would be “ axial ” 
part and “non-axial” part. From what has been 
said above it is clear that (a + bz) V represents a vector 
with a component aV along V and a component )V 
at right angles to V.__It thus represents a vector equal 
in magnitude to ./a* + 4? V, and making an angle 
tan—1b/a with V. But this vector can be represented 
by e°,/a? + 5? V, where 6 = tan-! b/a. 
Hence we see that 
a+ bi = re®; where r = /a?+ } 
and 6 = tan~!b/a; a result that can be obtained 
immediately by writing a + bi in the form 
el a bie | 
BEA Ramin prank 
The quantity 7 is called the ‘‘ modulus” of the 
complex quantity a+ bi, and the quantity @ its 
“ argument.” ‘The modulus of any quantity Q is often 
denoted by | Q|; thus if Q is a real negative quantity 
| Q| denotes its numerical magnitude taken positively. 


SomMeE TRIGONOMETRICAL THEOREMS 


It is convenient to place here some results which 
we shall require later so as to prevent the argument 
then from being interrupted by lengthy digressions. 

The formule for the resolution of the product of 
two sine or cosine terms into the sum of two terms are 
so frequently required that we will first state them 
Inere for reference. ‘Ihey are 

cos A cos B = $ cos (A + B) + $cos (A — B), 
sin A sin B = $cos(A — B) — $cos (A + B) 
and sin AcosB = $sin (A + B) + #sin (A — B). 


If n is any integer we have the well-known results that 


an an 
\ Costas == O1..... (A) and { sin #xvdX = Ol. «2 (A) 
° 


fe) 
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. es : sin nx 
For the indefinite integrals are respectively ht 


and — ©°S”* and both these expressions take the same 
n 


values when x = 22 as when x = o. From these results 
we can easily deduce that, if m is also an integer, 


27 
\ cos mx cos nxdx = 0 if m + n and 
@) 


—nifm=n . . . ~. (B) 
27 
\ sin mx sin nxdx =o if m -— nm and 
. = if mimi is x 1 eee 
27 
and \ sin ma cos nxdx = 0 in-allvcases 3) 3.) 
Oo 


To prove the first, we write cos mx cos nx = 
% cos (m+ nx +4 cos (m— n)x; and, since 
m + n and m — mare both integers, the result follows 
at once from (A) unless m =n. In this case, the 
last term is $ and the integral of this from o to 
2a is. ‘The result (C) follows in an exactly similar 
manner. For (D) we write sin mx cosnx = Fsin(m 
+ n)x + 4 sin (m — n)x, from which it is seen by 
(A) that this is zero in all cases, even when m= 1n; 
for in this case the second term vanishes. 


Some summation results, of which the integral 
results just given are the limiting case, will now be 
given. If p is any integer we have 


e082 caste eee ae etente eres ald C= 
iP p 
or as it is professionally written, 
P 
Sas aay (E) 
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The g under the ¥ indicates which letter the 
summation is with respect to, and the 1 below and 
the p above indicate the limits of g. This result is 


obvious geometrically, for ~ is the external angle of a 


regular polygon of p sides ; so that the sum in question 
represents the sum of the projections of the sides, each 
of unit length, of a closed regular polygon on its base, 
which we have taken as the zeroth or pth side. 
Similarly, 
p 
yin noe Hi et ari) 
gente? 
for this is the projection of the same sides on a line 
at right angles to the base. 
Further, if 7 is any integer, 
5 age, Sins dir, hea (G) 
g—t Gi Pp 
unless 2 is a multiple of », when the sum is p. For if 
m is prime to #, the first sum still represents the sum 
of the projections of the sides of the same polygon on 
its base line, but taken in the order mth, 2th, 3mth, 
etc., and in this sum each side is included once. If 
nm and p have a common factor 7, it can easily be seen 
that this sum will mot include all the sides of the 
polygon, but that every rth side will be taken 7 times. 
But the p/r sides that are taken would, by themselves, 
form a closed regular polygon of p/r sides, so that the 
sum of their projections is still zero. The only 
exceptional case is when 1 is a multiple of 9, then 
each term in the summation represents the projec- 
tion of the base and therefore the sum is p in the 
first case and zero in the second. 
We can now prove that if m and m are integers not 
greater than p/2 both 
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p . 
3S cos a oe (H) 
q=1 P iL 
bP. 20mg . 20ng 
d > sin sn 0 
i SER Y, p 2 


if m + n, while both sums are equal to p/2 if m = 1, 
save in the exceptional case when m= = p/2, 
when the sum in (H) is p and the sum in (I) is zero. 
These results can readily be obtained by replacing 
the products by the sum or difference of two cosines 


of the form cos ain we The restriction on m 


and # is necessary to prevent m+ » becoming equal 
to p or a multiple of it when m + n. 

Similarly, we can prove that with the same limita- 
tions on m and n 
2am 200n 


sin tao. . 


Lastly for this Introduction the definite integral 


in all cases. 


oe) 

sin bx 5 ; ; ; 
pene dx is very important in the theory of Fourier’s 

(e) 


series. It may be evaluated as follows : 


The indefinite integral \e*dx = oe of course, 
c 
whether ¢ be real or complex. Let us write c= — a 
-+ 1b, where a and b are positive, then we have, 
oe) Sayasep 
\e — (cos bx + i sin bx)dx = Ee Al zi gu say 
5 —a+ib —a+tib 
re) 


a+ 1b a+ 1b 


CERCA cone 
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Equating real and imaginary parts we have, 


\e-# COs ba dx ye 

fe) a* =e 52” 

and Woe sin Oviax = Meas! 
fe) a = b 


In the first result, the integrand is a “‘ function ” of 
b and the given result holds for all values of b. Let us 
suppose, then, that we write down this result for 
an infinite number of values of b between o and B 
each differing by an equal infinitesimal amount db 
from the preceding value and add the whole set of 
equations after multiplying by db. ‘The expression 

B 


under the integral sign then becomes \¢~ 4 cos bx db, 
te) 
4 — ak of 
which = =" sls 


, and the right-hand side becomes 


B adh 

{ae Sy ge Which = tan—" B Ja. 
fe) 
Hence we have 


(ea be ae ep a 


o x 


If now we make a = o, the integral is still conver- 
gent and the result becomes 


[o-e) 
sin Bx Tt 
\ doi — tan P00. 


oO 


Example.—By integrating both sides of this result 
with respect to B from o to 24 obtain the result 


Paes b 
sin” 0X IU 
12) a 2 


CHAPTER I 
PERIODIC FUNCTIONS IN GENERAL 


Classes of Functions. Mathematical functions 
may, for the purposes of the practical mathematician, 
be divided into three classes, viz. (a) What are known 
as “ Analytic Functions ” which we may perhaps call 
“Natural Functions * such as, 92%) sin) w.e- (0ce. = 
etc., and combinations of them. ‘These functions exist, 
not only for all real values of x, but also for all complex 
values of x of the form x + 1y.* 

(b) What we may call ‘“‘Artificial Functions.” These 
are functions which are arbitrarily defined as being 
equal to one analytic function in one part of their 
range and equal to another one in some other part. 
For instance, the function defined by 


f(x) =o when «x << — 1 


f(*) =bwhn—-1<xe<4+1 
and fiw) = owhenx >+1 


is an artificial function which is hereby defined for all 
real values of x.f Artificial functions must only be 
taken as existing in the regions for which they have 
been defined ; and thus must never be taken as existing 
for complex values of « unless they have been specific- 


ally defined for such complex values. (c) The third 


* Some complicated analytic functions only exist over definite 
regions of the wy plane, but this does not affect the argument. 

{ Strictly speaking, we have not defined the function for the values 
*=-+1. We could have got over this by saying f(x) = 5 when 
—I<*=<+ 1, where < means is “less than or equal to.” For 
practical purposes, however, this is an unnecessary bit of hair splitting. 

2 
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class are known as “ Empirical Functions ”; these are 
functions which we arrive at as the result of some 
experiments or observations where we obtain the 
result in the form of a graph or a set of points. In 
either case, the ideal function which we had in mind, 
and which it was the object of the experiment to deter- 
mine, is only known to a limited accuracy owing to 
necessary experimental errors. We are supposed not 
to know the “natural” or “ artificial ” function which 
represents this ideal function, as otherwise we should, 


Fic. 3. 


of course, use this representation, instead of the experi- 
mentally obtained values, in any mathematical opera- 
tions we had to perform on it. 

It is almost entirely to artificial and empirical func- 
tions to which we desire to apply Fourier’s Theorem 
and the method of Harmonic Analysis. 

Discontinuities of Artificial Functions. The 
most important distinction between an artificial func- 
tion and an analytic one, is the possession by the former 
of what are called “‘discontinuities,”’ the nature of which 
we must now briefly study. The heavy broken line in 
Fig. 3 represents an artificial function which possesses 
discontinuities at the points where « = a, b, ¢ or d. 
At x = a there is an abrupt change in the magnitude 
of the function represented by A, A, ; at x = 3, there 
is an abrupt change in slope but no change in magni- 
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tude, while at C there is an abrupt change in the rate 
of change of slope, or curvature, but no change either 
in magnitude or slope. ‘Thus at any point of discon- 
tinuity of a function there is an abrupt change, either 
in the magnitude of the function, or of some of its 
differential coefficients. In general, if there is an 
abrupt change in the magnitude of any of these quan- 
tities, there is also an abrupt change in the magnitude 
of all the higher differential coefficients at the same 
points unless these are zero. It is accidental, so to 
speak, if a function changes in magnitude, but not in 
sloperasiat rate lias. 

The discontinuities are measured by the magnitude 
of the abrupt changes at them. ‘The ordinate of the 
point A, is denoted by /(¢ — 0) and that of A, by 
f(a + 0) and the magnitude of the change, f(a + 0) — 
f(a — 0), we shall denote by Iu so that Ia stands for 
the magnitude of the abrupt change in f(«) at x = a. 
Similarly the slope of the arc A,B at B is denoted by 
f'(6 — 0) and the slope of BC at B is denoted by 
f(6+ 0), thus f(b + 0) — f(b — 0) is the magnitude 
of the discontinuity of f(«) at « = b, and we shall 
denote it by I. Similarly, the magnitude of the 
discontinuity of f’(#) at x = c will be denoted by 
I”, and so on. 

Periodic Functions. In the great majority of 
practical cases the independent variable in these 
functions is either the time or an angular magnitude. 
It will thus be more convenient to denote it by ¢ than 
by «. A function f(#) is said to be a periodic function 
when a constant 7 exists, such that f(Z + 2) = f(2) 
for all values of t. Changing ¢ into J + t we have 
K2T + 4) = f(T + 4) = f(a); and in general f(nT + 2) 
= f(t) where m is an integer. ‘The least value of T 
for which the equation {f(T + 4) = f(é) is always 
satisfied is called the period of the function. Any 
multiple of the period is said to be a period of the 
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function. ‘The behaviour of a periodic function over 
the range from t= — © to t= + is thus a 
continual repetition of its behaviour over the range 
from o to f. 

It is generally convenient, in discussing periodic 
functions, to change the scale of the ¢ axis so that the 
period is 27 as this will be found to make a great simpli- 
fication in all the expressions. Actually, of course, in 
a diagram, we should mark the period divided into 
360° and not into 2m” radians. Sometimes it is 
advisable to take an aliquot part of 27 as the period. 
For instance, the length of the radius vector from the 
centre of an ellipse has an angular period of x, and it is 
better to keep it such than to change the variable. 
But we shall suppose in future that the period of a 
periodic function is always 2% unless otherwise expressed. 

The student must grasp clearly the difference 
between a natural (or analytic) periodic function and 
an artificial one. If a periodic function of period 
2m is represented by one and the same analytic expres- 
sion over a whole range of 27, the function will never- 
theless be an artificial periodic function, unless the 
same analytic expression also represents it for all 
values (to + 0) outside this range. For instance the 
periodic function, of period 22, represented from 
t= oto t = 2a by e” is an artificial periodic function 
with discontinuities at t= 0, + 2%, + 4m, etc.: so 
too is a function of the same period and represented 
from — 2 to + x by sin pt unless p 1s an integer ; for 
it is only in this case that the periodic function we 
have defined is also represented outside the range from 
—z to +2 by sinp#. As another illustration, the 
function of period 2x defined by 

I 


OO aera 


is an analytic periodic function ; and so, it will easily 
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be seen, is the function defined by any rational 
expression which only involves sines or cosines of 
integral multiples of ¢. 

Four Special Classes of Periodic Functions. 
Besides the general periodic functions, which have no 
symmetry, there are four important special classes 


Fic. 4. 


characterised by different degrees of symmetry which 
enables them to be easily distinguished from one 
another by the eye. We have: 

(a) Class I called ‘‘ sine-harmonic functions.” <A 
periodic function belongs to this class when it 1s 


sere -7 0 vis 27 3c0r 


Fic. 5. 


possible to choose the origin so that f(t) = — f( —?2) for 
all values of t. ‘The function is then an odd function. 

Writing a + ¢ for ¢ in the equation f(4) = — f( —d,, 
we have fa + #) = — f(—a — 2d, = — fa — 0, 
since 2 is a period, and hence the function is also an 
odd function when the origin is shifted half a period. 
An example of this class of periodic function is illus- 
trated in Fig. 4. 
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(2) Class II, called “ cosine-harmonic functions.” 
A periodic function belongs to this class when it is 
possible to choose the origin, so that f(t) = + f(— 2) 


for all values of t. The function is then an even func- 


Eire. 16: 


tion. Writing a + ¢ for ¢ in this equation, we have 
fe +4) =f(—2—- ),=—/@ — 2), since 27 is a period ; 
and hence the function is also an even function when 
the origin is shifted half a period. An example is 
illustrated in Fig. 5. 


0 Tv i) 7 
(2) 2) 
Bean / + 
0 a 20 (e) we 
(e) @) 
Fic. 7 


(c) Class III, called “ odd-harmonic functions.” A 
periodic function belongs to this class when Kt + 7) 
— — f(?) for all values of t, when z is half the period. 
This relation is independent of the origin. An 
example is illustrated in Fig. 6. 

These three classes of functions can be obtained by 

c 
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repeating any given half wave in different manners, 
see Fig. 7, where (a) is supposed to represent a half 
period wave. 

If the second half wave be obtained by rotating the 
first half through 180° about the point 2, we obtain a 
Class I function as in Fig. 7 (2); while if it is repeated 
by reversing it fore and aft as in Fig. 7 (¢) we have a 
Class II function. If the repetition is obtained by 
rotating the half wave about the axis of #, and then 
displacing it half a period along this axis, we have a 
Class III function, as in Fig. 7 (d). A mere repetition 
of the half wave, without any turning about at all, 


BiGaeos 


would, of course, result in the period being no longer 
an, bUt sm. 

(d) Class IV, called‘ bi-symmetric functions.” These 
are the periodic functions with the highest possible 
degree of symmetry. Each quarter period is geo-. 
metrically similar to every other quarter period. An 
example is illustrated in Fig. 8. It will be seen that 
the functions of this class belong to all the three 
previous classes. With the origin as shown the func- 
tion is an odd function, if the origin were shifted to 


7 . . 
. the function would clearly be an even function; 


moreover, the function is also obviously a Class III 
function. ‘The student should prove to himself that 
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any function that belongs to any two of the first three 
classes necessarily belongs to the other one also, and 
is a Class IV function. 

It is important to be able to recognise to which of 
these four classes any periodic curve that occurs in any 
problem belongs. Let us consider, for instance, the 
open circuit voltage wave form of an alternator. 
Since in practice a “north” field pole is made 
geometrically similar to the adjacent “ south” one, 
and since both are symmetrical about their centre 
lines, it is clear that the voltage wave form must be a 
Class IV function. This is not true, however, when 
the machine is on load, for then, as is well known, 
the field due to the armature current weakens that 
due to the field poles on the sides at which the armature 
is approaching the poles and strengthens it on the 
leaving sides. Nevertheless, it is clear that after half 
a period, the induced voltage will still be numerically 
equal, but of opposite sign, and thus that the voltage 
wave form is a Class III function. 


Example 1.—Prove that any periodic function can be represented 
as the sum of an odd and an even periodic function. 


We have f(Z) acid heel ier t) Tele. au (Se ?) 


a identically. 


The first term merely changes sign on changing ¢ into — ¢ and the 
second one is unaltered thereby. Hence the first term represents a 
Class I function and the second a Class II one. 


Example 2.—Prove that any periodic function can be represented 
as the sum of a Class III function and another of half the period. 


fli) + flt+ 7) 
2 


She: 
eae fo le a nea 
The first term merely changes sign on changing ¢ into ¢-+ ™, since 

2m is the period. The second term is unaltered by this change and 
hence ~ must be a period of this term. 


Example 3.—Prove that any Class III function can be represented 
as the sum of two Class IV functions. 
For if the function be represented as the sum of an odd and an 
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even function by Ex. 1 it will be seen that both these functions are 
still Class III functions and must hence be Class IV functions. 


Example 4.—(a) Prove that ay + 4, cos t+ ay Cos 2¢-> dg Cos 3t 
+... represents a Class II function of period 27. 
(b) Prove that b,sin ¢+ b,sin 2t+ bg sin 3t-+ . . . Tepresents 


a Class I function of period 2n. 


(c) Prove that a, cos t-++ ag cos 3t-+ 45 cos St . . « b, sin ¢ 
+ b, sin 3t-+b, sin st+ . . . represents a Class III function of 
period 2r. 


Each of these series is supposed convergent. 


Example 5.—Prove that either of the series 
(a) a, cos t+ ag cos 3t-+ a,cos S¢+ .. . 
and (b) b,sin t+ d,sin 3t-+ b, sin 5t-+ . . . represents a Class IV 
function of period 27. 
For the first series represents a function belonging to Classes II and 
III and the second a function belonging to Classes I and ITI. 


If the origin in the second series be shifted to the point ¢= z by 
2 


“oF Tt : 
writing => + 7? the series becomes 


b, sin (E+ ) + bg sin ("+ ')+ b, sin (E+ ) ae 
which 
= b, cos t’ — bg cos 3t’ + b, cos §¢7— . . . 
showing that both these series represent the same kind of functions 
teferred to different origins. 


Example 6.—Sketch the curves given by 


(a) f(t) = cos t + ‘2 cos 2¢. 

(b) f(é) = sin t+ °2 sin 2¢. 

(c) f() = sin t+ *1 sin 3¢. 

(2) f(t) = sin t — ‘I sin 32. 

(e) f(t) = sin t+ ‘I cos 32. 
and verify in each case that your sketch belongs to the correct class 
of function. 


CHAPTER II 
FOURIER’S THEOREM 


Statement of Theorem. Fourier’s "Theorem 
states that any periodic function, whether “ natural” 
or “ artificial,” of period 2 2 can be expressed by the 
serles— 


f(8) = 2" + a, c08 t+ ay C08 2¢ + a, co8 34 + He | 


+ 6, sint + b,sin 2t + b, sin 3¢ + , 
(1) 


where 


—— 1 9 COS #% 
iN 40 ra ‘ 
and by = t \ 70 eed) 


There are a few unimportant restrictions on the 
behaviour of f(t) which we will notice later in order 
that the theorem should be true. 

We cannot at this stage prove Fourier’s theorem, 
for this requires the proof that the representation of 
the function by the series given in (1) is possible ; 
but we can easily prove that if f(z) can be represented 
by the given series the coefficients must have the values 


given by (2). The two sides of (1), being supposed 


sie (a . 
* The reason for writing —° instead of a, is so that the value of a, 
2 


will be correctly given by putting » =o in the expression for ay. 
21 
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identical for all values of ¢, must remain identical 
when the same operation is performed on both sides ; 
unless this renders the series divergent. Let us 
multiply both sides of (1) by dz and integrate from 
t=0o to t= 2x. By equations (A), p. 7, all the 
terms on the right-hand side vanish except the first, 
so we are left with 


2 2 


J Kod 2h dem na, 


which is the value obtained by putting = 0 in (2). 

Again multiplying both sides of (1) by cos n¢ dt and 
integrating from o to 2m we see that all the terms on 
the right-hand side vanish by (B) and (D), p. 8, 


except the a, cos nt term; thus we get 


2m 2m 
\ T (2) cos nt dt = an cos? nt dt 
ie} fe} 
2n 


an 


= a (1 a COs 2nt)dt = Wan, 
2 ° 


which is the value of ay given by (2). Similarly, we 
have 
27 2 
\ JC) SIN aide \ sin? nt dt 
10) ie) 
be 27 
= \ (1 — cos 2nt)\dt = aby, 


arg 


also in agreement with (2). Hence the theorem is 
proved on the assumption that the expansion is 
possible. 

In the expansion (1) the various terms are called the 
“harmonics” of f(#), the cosine terms being called 
the “cosine harmonics ” and the sine terms the “ sine 
harmonics.” ‘Thus a, cos t is called the “ first cosine 
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harmonic” and d,sin 2t is called the “second sine 


harmonic.” Both the uth sine harmonic and the nth 
cosine harmonic, of course, go through » complete 
cycles in the period of the functions. The two nth 
harmonic terms a» cos nt + by, sin nt can be combined 
together in the form 


r/an® + bn? cos (nt — $), 
where tan ¢ = dy/an. 
This expression is called the resultant or complete 
mth harmonic; »/an? + 5,2 is its amplitude and — ¢ 
its initial phase. ‘The “ first harmonic ” is often called 


instead, the “‘ fundamental harmonic.” ‘The a’s and 
b’s in (1) are often called the ‘‘ Fourier constants” of 


S(t). 


‘Two or three proofs of Fourier’s Theorem will be given later. The 
student familiar with the theory of linear differential equations with 
constant coefficients will be interested in the following. If f(z) has a 
period 27 we have f(t+ 2x)= (f(t). By the symbolical form of Tay- 


lor’s theorem f(t + 27) = e?*P f(z) where D= qb hence (e2™P—1) . f(t) 
= 0, which is a differential equation for f(t). We know that if 


m is a root of e2™P — 1 =o, the solution consists of the sum of a 
number of terms of the form C”*, For two imaginary roots, + im, 
the terms take the form A cos mt+ B sin mt, A, B and C being arbi- 
trary constants. Now e?*™) — 1 =o has obviously no real root save 
D = 0; to search for imaginary roots we put D = 2) and get e?7A = 1 
or cos 2mA= 1 and sin2mA=o simultaneously; from which it is 
obvious that 4 can be any integer. Hence we see that f(t) can be 
expressed in the form 
f) = C-+ A, cos t+ Ag cos 2#-+ Ag cos 3¢ + 
+ B, sint + B, sin 2t + B,sin3¢ 4+ . 

The constants must be chosen to satisfy the initial conditions, which 

in this case consist of the behaviour of the function from 0 to 27. 


Simplification of Fourier’s Theorem for the 
Four Special Classes of Functions. If /(t) is 
an odd function, so that f(t) = —f(— #4) for all 


values of ¢, we have 
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fo 4 g, cos t + a, cos 2t + e., +b, ite 
2 
= — {22+ a cos ¢+ a, cos 28+. ee 
D 
— 6, sin t— 6,sin 2?— . .; 
since cos(— #)=cos¢ and sin( — #)= — siné. 
This reduces to 


Gi, + 20) COS -j-420, COS 2) ene 
for all values of t; hence all the a’s are zero, so that any 
odd periodic function can be expanded in the series 
ff) =b,snti+ &sin2t+ hsinzgt+ . . (3) 
Similarly, if f(z) = + f(— 2) for all values of #, so 
that f(#) is an even function, we find that all the b’s 
are zero and that 


fd) = 52 + a, cos t + ay cos 2¢ + <a) 


Again, if f(z) satisfies the condition that f(t + 2) = 
— f(z) for all values of #, we have 


a 
sit COSA oi 5g COSTAE A 


a 


Dy sin? + 6, sin 2? -- 4. 
= — © — ay cos(t + m) — ay cos 2(¢ + 2) — 

@s COS'3(7.-|- ar) oe 
— bsin(t-+ 2) — b sin 2(¢ + 2) — bg sin3(t +2”)... 
= — “2 + ay cos t — ay cos 22 + a, COSTS to ea 


+ bsin ¢— bd, sin 2¢+ b,sin 3t— . 

and hence 
a, + 2a, cos 2t + 2a, cos 4¢ + 2a, cos 6¢ + 
+ 2b, sin 2¢ + 2, sin 4¢ + 2b, sin 6¢ + 


= 0 for all values of t, whence 
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aj—44—=4—4...=0andd,—b,=),..:=0, 
and the function can therefore be represented by 
f(t )= a,c08 t + ag cos 3t + 
+h sint+ bd, sinzgi+ . . . 5) 

It will be seen that these three cases correspond to 
Class I, Class II and Class III functions respectively, 
and the reason for the names “ sine-harmonic,” “ cosine 
harmonic ” and “ odd-harmonic ”’ will now be evident. 
We will now give simplified expressions for the 
Fourier constants in these cases. The expressions 
previously given for the Fourier constants, viz. : 


an an 
peo ‘| HQ) cos nt dt and by = ~\ (2) sin nt dt 
ce fe) me ie) 
can be written 
4 wT 
an = \ f(t) cos nt dt and by = i F(Z) sin nt dt ; 


for, since the expressions under the integral sign have 
2x as a period, the result will be the same so long as the 
range of integration is 2m, independently of the posi- 
tion from which the range commences. 

For a Class I function with the origin suitably 
chosen, f(z) is an odd function and since cos nt is even, 
the product f(#) cos mt is an odd function; hence 
the integral of it from — a to +a is zero. This also 
shows that in this case all the cosine terms vanish. 
Again, since sin mt is an odd function and the pro- 
duct f(#) sin nt is an even function, its integral from 
— x to + a is twice its integral from 0 toa. 

Hence the result that a Class I function of period 
2n may be represented by 


f@) =b,sint+,sinzt+. . . (6) 
where = “| CE) PUL HE CE ee a ks (7) 
Ho 
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Similarly for a Class II function with the origin 
chosen so that it is also an even function; 4» vanishes 
since f(z) sin mt is an odd function, while a» is given 
by double the integral over half the range. 

Hence the result that a Class II function of period 
2a may be represented by 


f(t) = 2 + ay cos t + a, cos 2¢ + Be seh Cod 


where an = a) ji) COS Rial = a) see) 


= oO 

Consider next a Class III function and break the 
integrals broken into two taken over the ranges — a 
to o and o to a respectively. If we consider corre- 
sponding elements in these two sections at a distance 
of x apart we see that f(#) has equal numerical values 
but opposite signs in the two cases, since f(t + =) = 
— f(z). Further, since a is equal to an integral 
number of periods of the even harmonic terms, the 
product f(#) cos mt or f(#) sin nt will, when 1 is even, 
have equal numerical values but opposite signs at all 
corresponding elements and hence the integral from 
— zx to o will cancel that from 0 toz. If, however, 
n is odd cos mt and sin nt will have an odd number 
of half wave-lengths and simply change sign in the 
range of a, so that /(#) cos nt and (2) sin nt will have 
the same values at all corresponding elements and 
we may therefore take the integrals from o to a and 
double them. 

Hence the result that a Class III function of period 
2a may be represented by 


f(t) = a, cost + aycos3t+ a,cosst+. . . 
+ by sin’? > b; sin 38 asic ee eee (10) 


where 
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an = 2\ f(2) cos nt dt 
a le) 
(11) 
and Dn = 2 f(4) sin nt dt 


fe) 


Lastly, let f(#) be a Class IV function with the 
origin chosen so that it is an odd function. Then 
since it is also a Class III function we have 


SQi=Opsin fe by sini. (12) 


where 
ee *\ (2) sin nt dt. 
2 Oo 


If f(2) is referred to a new origin at ¢ = > it becomes 


an even function with respect to the new origin; so 
too does sin nt when n is odd, for if m= 2m +1 it 


becomes sin (2m + (2 = r) which = cos {ma + 


(2m + 1)t’} = (— 1)” cos (2m + 1)?’, hence the given 
integral from o to a is equal to twice the integral 


It 
from © to —. 
2 


Hence the result that an odd Class IV function of 
period 2x can be represented by 


(iserOnsi | 0, si 3h 2. |. (12) 
where 
7/2 
ote 4 fi) sin ntdt =... (13) 
A 1@) 


In an exactly similar manner it can be proved that 
an even Class IV function of period 2a can be repre- 
sented by 

f(t) = a, cos t + a, cos 3t + a, cos 5t+. ~ (14) 


where 
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7/2 
an = A ft) COS That ie eS) 
aE O 


The student should guard against the mistake of 
working out these latter integrals when n 1s even and 
then because the result is not zero putting in even 
harmonic terms in the expansion of f(t). The only 
reason why the range of integration is halved is because 
the values of a, and by are zero when n 15 even and can 
be obtained by doubling the integral over half the 
range when 2 is odd. 

In functions of Classes I, I] or IV, only one kind 
of integral has to be evaluated, viz. [f@ cos nt dt or 
[fa sin mt dt; but in the Class III functions, and in 
the case of the general function, both types are required. 
In these cases it is most convenient to multiply the 
expression for b, by z and add it to the expression 
for an, so that 

Tw 


ayn + 1b, = \ f(A[cos nt + 72 sin nt]dt 


= 


~ 


iia 
= | foe ‘di: 


and then after evaluating this integral to separate it 
into its real and imaginary parts which are equated 
to ay and by respectively. 

If the period of f(z) is not 2% but an aliquot part of 


: 2n : 4 
it, say me the formule we have given for the coeffici- 
ents still hold, since 22 is still a period of the function. 
In this case all the harmonics vanish except those 


which are multiples of the mth,* so that 


* Actually, of course, what appears here as the mth harmonic is 
physically the first or fundamental harmonic. 
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a= 2 + a, cos mt + a, cos 2mt + 


+b, sin mt+ bsin2mt+ . . (16) 
where 
I 2m 
an =— \ f(s cos mnt dt; 
ue 1@) 
but since f(#) and cos mnt both have a period of 
2G es : 
a this may be written 


an 
mM m 
aps =\ GCOS nid iene es A177) 
22 ie) 


and similarly for dy. 


Example.—Prove that if the period of f(t) be T, Fourier’s expansion 
becomes 


ft) = “2+ a, cos + a, cos #7 + ee 


a 
+ 4, sin 7 + by sin se +- 5, sin et ee ees (LS) 
AR 
where ag = a f(t) cos a at Pee a ee i 00) 
eo) 
fi 
and ba = 7 | fC) sin ~_ Lie epee (20) 
fe} 


We may conveniently close this chapter with 


Some APPLICATIONS TO ELEcTRICAL ENGINEERING. 
Isolation of Particular Classes of Har- 


monics in any Periodic Function. If /(#) = “2 


Bera pecoge? a. sa Oy Sin 2 + ie .gr and if we 
displace it by a distance , the displaced function, 


f(t +2), is given by 
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flt+m)=2 + a,cos(t +2) + ele b, sin ea) 


Sees 5 C08 F --d,-COS 27 — GCOS 3i a ae 
2 
— b, sint + bd, sin 2t — bg sin 3¢ + 
If we add this to the original function we get 


f(0) + ft +2) = 2] 42 + ag 08 2f + ay cos 4t + sun 


+ bsin2t + sin get. } 


that is, twice the sum of all the even harmonics. More 
generally, if we displace the function by the distances 
males oe ey: se) in turn and add these 
Cite ae ee P 

~ —1 displaced functions to the original, we get, 
since the sum of all the sine and cosine terms which 


are not multiples of the pth vanish by (G), p. 9, 
the Sum = p | % + ap COS pi + a,pcos2 pit’... 


+ bpsin pt + bp sin 2pt +. . | 
ioc es eee 200) 


that is, we get p times the sum of all the harmonics 
which are multiples of the pth. Again, if we subtract 


t(: + =) from f(t) it will be seen that all harmonics 


which are multiples of the pth will cancel out. As an 
example, let us find 

(1) The Terminal Voltage of a Three- 
Phase Generator. Let f(t) be the voltage in- 
duced in one phase of a three-phase star connected 
alternator. Then the voltages induced in the other 
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phases are given by f(1 - = aad f(: + #), the 
latter of which may be equally well represented by 
f(r-2 

The aint between the second and third phases is 
f(: Fa =) es f(: £ =). Since f(#) is necessarily a 


Class III function (p. 19) we have 


F(t) = 4, COs t 4+ -ag-cos 3¢ + a,cos5¢+ .. 
+ 6, sint+ db, sin 3t¢+ .. 


and hence 

+8) 1") 
= a,}cos (¢ + = cos (1 — =} Ee 
++ by {sin (¢ + =) — sin (: 9) ee 


which easily reduces to 
R40 ae sin te a;sin §t — a,sin 7¢ + a, sin IIt. yf; 
+ »/3{b, cos t — b, cos 5¢ + b, cos 7t — by, cos 11 t + ..}. 


5 A . 4 d 
If in this we write sin ¢ = — cos a an 


. TU 
cos ¢ = sin (: Lie =), etc we get 


/3{a cos(¢+)— a, COS (142) +a, cos 7(1+2)— : 
+5, sin (: +2) 4, sin (142) +4, sin 7(1+7) = 
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If we compare this with f(#) itself we see that all 
the harmonics which are multiples of 3 have dis- 
appeared while the amplitudes of all the others are 
multiplied by 4/3. We note also that the coefficients 
of the sth, 11th, 17th, etc., harmonics have their 
signs changed, which may be interpreted as a change 
in phase of the harmonic in question amounting to 
half its period. 

(2) Voltage Wave Form of an Alternator of 
given Field Form. By the “field form” of an 
alternator we mean the voltage wave form induced 
in a single conductor. This is not in general the same 
as the voltage wave form at the terminals of the 
machine, since the armature generally contains a large 
number of conductors connected in series which 
occupy a displaced position relatively to one another 
on the armature. We will suppose that the windings 
are uniformly spread over an angle of 2« radians, or 
of 2a ‘electrical radians ”’ if the machine has more 


than two field poles. If 

F(t) = a, cost + ag cos 3t+ . 

+. b, sinz+ db, sin 3+... 
represents the voltage induced in the central con- 
ductor, then 
f(t + 9) = a, cos (¢ + 0) + as cos 3(¢ + 0) 4+ . 

+ b, sin (¢ + 6) + bs sin 3(¢ + 6) + . 
will represent the voltage induced in a conductor at 
an angular distance @ in advance of the central one, 


and the total induced voltage in all the conductors 
will be proportional to 


\ + 6)d0, 


Late is. eto 
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a, | cos (t + Odd + ay | cos 3(¢+ @)d0@+ ... 
= —o 


Ls {sia Cee pee 
This is equal to 
a,{sin (¢ + «) — sin (¢ — «)} + Ffsin A fete) 
sin 3(¢ — «)} + d{cos(¢— «)—cos(¢+a)}+..., 


or to 
: 2Ger- 
2a, sin « cos ¢ + = sin 3% cos 34 .. 


+ 25, sin & sin t+ EP sin 3x sin 36 + re 


If we divide this expression by 2 sin a, so as to make 
the amplitude of the first harmonic equal to its 
amplitude in the field form, we see that the amplitude 
of the mth harmonic in the voltage wave form is 

. sin na : : 
reduced (relatively) to ———— of its value in the field 
n sin a 


form. 
If, for instance, a = 2/3 or 60° all harmonics that 


are multiples of 3 will vanish and the field form 


iii a, GOs i=. 2, COS 3% 2), 
+b, sint+ 5, sin 3t+ ... 


will produce the voltage form 
a a 
a cos f— —2 cos 47 -+- —L cos 7? —... 
1 F 5 7 7 
; bar’ eee 
Pg EE lee ea a sin 7?— . 


D 
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Example.—lf «= = or 72° so that the spread of the windings is 


144°, prove that the voltage wave form of a star connected three-phase 
alternator will be proportional to 
a, cos t+ *08Qa_ cos 7t + “0914, Cos I1t — *0484)3 Cos 13¢ Se 
+ b, sin t-++ ‘089, sin 7t-++ *ogtd,, sin 11¢ — *048),, sin 137 +. 
where 
a, cos t-++ ag cos 3t-+ 4, cos St-+ .. . 
+ b, sin t-+ bg sin 3¢-+ bs sin St+ .. . 
represents the field form. 
[There is a difference of phase between these two expressions as 
written, but this does not affect the amplitudes of the various 
harmonics. | 


(3) The Root Mean Square (R.M.S.) Value 
of a Periodic Curve. We frequently, especially 
in electrical engineering, require to know the square 
root of the mean square of the value of a periodic 
function; f(#) being a periodic function of period 2z, 
this quantity is clearly by definition, 

ae bye 
hae SY) at 
Let us multiply both sides of the equation 
f= 2+ a, cos t + dy cos 2t + 


+ db, sin ¢+ b, sin 2t + 
by f(z) dt and integrate from 0 to 2z. 


27 
Since \ i (2) cos nt dt = nan 
[e} 
ane \ Ovni ch 2k, 
ie} 
we get 


J {08 dt = S(na,) + aylnas) + alan) + 
Lh; Geb) a baler) Cee 
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and hence 


27 
=i Oa _ jo RSOUe ibe melee 
ie} 


2 2 


If we denote the constant term . by c, and the 


amplitudes of the rst, 2nd, 3rd . . . harmonics by 
Lena bai ete 0 sR WESES she the R.M.S. value of a 
periodic function is given by 


Ae ape n ed A> Ga" eee ee 


2 

(4) The Form Factor. Another characteristic 
of a periodic function which is required in electrical 
engineering is the “ form factor” of a voltage wave. 
This is defined as the ratio of the R.M.S. value to 
the mean value of the rectified wave. 

Since the voltage wave of an ordinary alternator 
can only involve odd harmonics, we limit ourselves 
to that case. Let the voltage be given by 

Vee psin t4- b, 31 3t -—- Cs 
Gy COR. t= Ge COS Sf 4-4). 
so that its R.M. S. value is 


(GCOS EERE 


We will suppose that ‘V =O whens 90 atid that 
b, is positive so that the positive half of the wave 
extends from t=otot=a. This of course neces- 
sitates the relation a4, + 4,+4,+ ...==0O among 
the cosine coefficients. 

The mean value of the positive half of the wave 


Tw 
is a) V dt, which is readily seen to be equal to 
vt 1@) 


Z by. Os ) 
at atet ..-), 
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So that the form factor is 


iy NG ere a eee 
2/2 Pee Uhe 3 oe oo 
The numerical factor is equal to I-11I approx. 
Whenever the equation of a periodic function is 
known, both its R.M.S. value and its Form Factor 
can be deduced directly from the integrals defining 
them : it is quite unnecessary then to use the expres- 
sions for these quantities given above in terms of the 
Fourier constants. 


CHAPTER Ii 


ON THE REPRESENTATION OF ARTIFICIAL FUNCTIONS 
BY FOURIER’S SERIES 


We shall in this chapter obtain the Fourier Series 
representing a variety of fairly simple artificial func- 
tions. The Fourier series representing a function is 
called briefly the ‘‘ Harmonic Analysis” of the 
function. 


Problem I.—Obtain the harmonic analysis of the 
eveniunction represented by /(7) = 4 from #= 0 
PORE tartan DY (tf) Orlronny? — a tori i* 

Since the function is even, we have 


f(t) = “2 + a, cost + ay cos 2t + oe tes eG) 


where He a) S(4) cos nt dt. 
a 1e) 


If we substitute the given expressions for f(t) over 
their appropriate ranges we get 


2 (04 
ans A 4 cos nt dt, 
Ag 1) 


; : 2Das 
which gives Gn = sin na. 
TH 
* It will be assumed in all cases that the period is 27 unless other- 


wise stated. 
37 
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This readily gives all the a’s save a, which takes an 
. . Oo ° ° 
indeterminate form of = If, however, we imagine 


n to be very small but not quite zero, we see that the 


2h na 2ha 
expression tends to — =", when # is zero. 
UH qt 


Fic. 9. 


Hence substituting the values of the a’s in (1) we 
get 


4 ; 
fa = 2a + 0 = C08 ¢ =p “= cos 2¢ 


sin 3% 


aC cosstt TAS) 


In particular, if = : we get 
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ye Cosi oe ee | 
Ae) = 2A + 4 }. @) 
A sketch of this function is shown in Fig. 9(a) 
where the dotted line indicates the value of the first 
term. The way in which the successive terms 
approximate to the portion ABCD is shown enlarged 
in (6) where the curves marked 1, 2, 3 and 4 indicate 
the sum of the first 2, 3, 4 and 5 terms of (3) 
respectively. 


Fig. <0. 


paras 
We may note here that the constant term me is 


always the mean height of the function ; for the mean 


27 
height is z\ f(t) dt, which is the value of 72, 
2s 2 


Problem II.—Analyse an even Class IV function 
which is equal to 4 from ¢ = 0 to t = a and to zero 


from t=ator—=". 
Here 
id) —"aec08 F 4, C08, Zt | 
where 
n/2 0 
4h . 
an = 40 cos nt dt = 4 cos nt dt = van Slt Mee 
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Hence 
ben 2 Paar si 3G eve 
Ol teeRe creer 


SEES reece 4. 2 2 lhe emt) 
The graph of the function is shown in Fig. to. 


The part from = to z is constructed from the part 


Fic. II. 


from ae to o by means of the relation f(t + 2) 
= — f(t), which holds for all Class IV functions. 


Problem III.—Analyse an odd function represented 
bt a —t 

bya7 (7 = from (f= OVfOrt = Gana by. b= =>) 
from ¢=a to t=a. ‘This function is represented 
by Fig. 11. 

Unlike the two previous functions, this has no dis- 
continuities in magnitude but only in slope. 

Since the function is odd, we have 


f(t) = b sint + b, sin 2¢ + bg sin 3t+ .. . 


where basse 2(" f(t) sin nt dt. 
Uo 
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If we integrate this expression by parts we get 


Ca : [Ease “| pn : \ fo cos nt dt. 


Tt TN Y 6 


This step would obviously be justifiable if f(t) was 
asingle function. As it is not, we must break the range 
of integration up into two, one from o to a and the 
other from a toz. If we do this we see that the value 
of the first term at the upper limit of the first range 
is the same as its value at the lower limit of the second 
range and hence these cancel, leaving only the values 
at o and m as written. We might have said simply, 
that the first term only takes cognisance of the mag- 
nitude of f(t); and since f(z) has no abrupt changes in 
magnitude it behaves, as far as the first term is con- 
cerned, as an ordinary analytic function. Now /(2) is 
zero both when ¢ =o and when ¢=2a. Hence the 
first term entirely vanishes. ‘The second term, the 


integral, breaks up into two; for f(z) = : from 0 toa 


and = — from « to z. 


Hence we have 


hee 21° cos nt dt — ae {cos nt at| 
tT Oo 


AAtAR : 


2h(sin na sin no ! 2 A mo 
nl no ne — a) a(% — «)n? 


sina . sin 2% sin 30 


2h : 
it) = ual a sint +. = in 2t + Fo sin 3¢ 


fe eC) 
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Some special cases of this result are worth noticing. 


If « = = we get 


8h{. sin3¢ , sin5z sin 7¢ 
1G) Se sin t— - a + 7 ec: ee s) 


so that the function is a Class IV function, which is 
geometrically obvious from Fig. 11 when « = 2/2. 

Again, let us put a = 0 so that the function takes 
the form shown in Fig. 12. 


BIG. 12. 


The expression (5) requires evaluating as an indeter- 
minate form, since both numerators and denominators 


- «SIN HX 
are zero. However, on replacing ——— by m when 
a 


& = O, we get at once 


fe) = loins 4 Si2# , singe , sin ge 


| 
oes) 
satay) 


Here the question arises, What is the value of this 
series when t= 0? If we put ¢t = 0 we see that the 
sum of the series is formally zero; but it really takes 
an indeterminate form, as is seen by imagining ¢ 
infinitesimally small when we get 
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f= Pie rtetetet 1 ba BE 


which is of the form 0 X o. 

The value of the function at t = 0 is also indeter- 
minate: it is obvious from Fig. 12 that when ¢ = 0 
the function possesses all values from — 4 to + bh. 
We shall prove subsequently (Chapter VI) that the 
sertes is also indeterminate when ¢ = o and that the 
range of indeterminateness of the series is 1-179 times 
that of the function. 


Problem IV.—Wet us analyse an even function 
which is represented by the same equations from 
o to x as the function in Problem III. 


Eleven ( (7) — 2 + a,cost+ 4,cos2t+.. 
where ae \ S(t) cos nt dt. 
Ae @) 


If this be integrated by parts as in the last case, 
we get 


=2(@s sin | -2\r (z) sin nt dt. 


As before, the first term vanishes at both limits, 
and the second term gives 


Tw 
I , 
a, = — a {EA sin me dt — 1, sin ne ar} 
O boy eG 


7H a 


2h {1 —cosmu , cos mm — cos nu 
—— ae : 
na a ™% — a 
This takes different forms according as to whether 
m is odd or even owing to the cos mm term. If 1 is 
4h sin? no /2 
a(n — a) ie 


even it reduces to 4d, = — 
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and if 2 is odd, to 

By fe = lL + eat 
ae = a) | n* 

The value of a, which is twice the mean height of 
the function is obviously 5. The shape of the function 
is represented by Fig. 7(¢), p. 17. 


an = 


Problem V.—From the two previous problems we 
can at once write down the analysis of the Class III 


: ¢ a t 
function which is equal to bt from t = 0 to a and to 
a 


ba — ; 

Aa — 4) from t=atoz. For it must be represented 
te = OF 

by 


Jit) = Gy Cost 1-0dy COs Fe 
+ 6,sin t+ 6, sin3jt+... 


and the values of these coefficients have all been 
found in the last two Problems, giving the result : 


(Oe [Foor sin ¢ + SRS sin 30 + ae 


a(za)L 1? 2 
2a 
—-— I+ cosa ate cossa 
by EA 
+ -=T $3 cos ¢ -- 
I 33 
COS tA Ei) ae |. 


The shape of the function is represented by 
Fig. 7(d), p. 17. If in this series we put a =~ 
2 


> 


all the cosine harmonics vanish and the function 


becomes, as is geometrically obvious, a Class IV 
function. 
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Problem VI,—Analyse an odd Class IV function 


which is equal to bt from ¢ = 0 toa, and to j from 
a 


t =a to ae 
D 
Here 
yf 0, ait t= basin 37 -E 


where 


7/2 
bn = 4 A(2) sin nt dt. 


Fre. 13; 


_ Integrating by parts we get 


7/2 t/2 
ios gp) eek \ FC) cos nt dt. 
ie) 


ae ne Ma 


This is justified, since f/(#) has no discontinuities in 
magnitude. ‘The first term vanishes at both limits; 


b 
also f(t) = - from o to a and is zero afterwards ; 
a 


hence 
Peds, Lee na — 42 Sin na 


an an man? 
so that we have 


Oe ae each ee ae sin 3¢ 


sin aN 


a sin §¢-+ . as poet eS) 
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V3. 


a : 
If, for example, we take a = x we get sina = SS 


/3 


Sif30s-—O -aeciNh 0 ee Sa etc., and so 


Ow Bee Si SP NGID, 77 es Ee 
Ors il p{sint — i eis oe ke Pilar } 
or f(t) = b{1-056 sin t — -042 sin §¢ + -022 sin 7¢ 
— 009 sin IIf = 224s 


This function together with its first harmonic is 
shown in Fig. 13. 
Problem VII.—Analyse an even Class IV function 
ancora) ag B(x — ) ELOTU ies Omt ON Stes = 
IU 


Here f(t) = a; cos # + a, cos 32 -+ ==. 


7/2 
where an = a JG) cos’ nt di. 
ve ie} 


Since f(z) has no discontinuities in magnitude we 
have, on integrating by parts, 


ae BIA * en/2 
is t ‘fr : 
Pees 2 pe i IF mat \/ (2) Siniwe dit al} 


aD 
The first term vanishes at both limits; /(#) vanishing 
7 : 

at © and sin mt at zero. Since f(#) has no discon- 


tinuities in magnitude, we may integrate (9) again 
by parts, giving 


zea | usc eb Ey 


7/2 
F'(2) cos nit dt. 


= 2 
tn nN TH oO 


Again the first term vanishes at both limits. The 


integral can easily be evaluated, since f@O=- a 
wt 
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Lh 
from 0 to z Hence we get 


Ps 
32h (7/2 22 ee TT 
an = - 5, \ cos nt di = 2_ sin — 
"82 i n3_3 1 3° 
Whe “sine term i¢° 4- 1 and. => 1 alternately for 


fi—=ai, 35 5. etc, and so 


(o= 324 \ cos ve coe op SES OE | (10) 


os 7 
This function is shown by the heavy line in Fig. 14. 


Differentiation of a Fourier Series. If we 
differentiate (10) we get 

7h ease ee sin 3¢ sin 5¢ 

T@ =| sin ¢ + 2 Fe eisai 

But /(¢) is represented by the thin line in Fig. 14. 
The maximum numerical value of this function is 

| PAGES 9301 9 : 
easily seen to be 7, 3. Since Tt 3 from 0 to > 
If we call this maximum value H, the expression 
for f(t) becomes 

: Site. si0esi «cin St ; 

fij=- sin t — 2 _ oye aa \ 
which is in agreement with (6), p. 42, for the same 
function. 
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If we differentiate this again we get 


” 8H COS Zhi 4. COSTST | 
fii) = — cos # — : + : aoe 
but f(z) is the function represented by the dotted 
lines in Fig. 14, consisting of lines parallel to the 


° e It 4 
axis and at distances H ~ : above and below it. If 


this distance be called H’, we have 


fO=- aH | cos — ee Ne a 
7 3 J 

which is in agreement with (3), p. 39, if the axis be 

shifted parallel to itself so as to make the constant 

term zero. 

These examples illustrate the theorem that we may 
differentiate a Fourier Series term by term and 
provided that the differentiated series 1s convergent the 
result will represent the expansion of the differential 
coefficient of the given function. This result hardly, 
we think, requires further proof than is afforded by 
this example. Conversely, we may in a similar manner, 
integrate a Fourier series and the result will represent 
the integral of the given function. No qualification 
about convergency is necessary here, since the inte- 
grated series is necessarily more convergent than the 
original one. 

All the functions we have analysed so far have 
belonged to one of the four special classes. We now 
give one example of the general type. 


Problem VIII.—Analyse the function given by 
f(t) = bt/x from t=0 to a and by f(t) =o from 


i 3 gh (OY Qi 


Since there is no symmetry about this function, we 
have 
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f(t) = 52 + a, cos t + ag cos 24 + ies 
+ bh snt+ &sin2zt+... 


27 
where fag tHe es \ f(pem dt 
7 e} 


Tw 
iy Ear ad. 
i" Je 


Integrating by parts, we get 


int YT wT 
an + tb, = e E ‘| = b ZUras 
ie) 


m?2L in 21n 5 


int 
= b we h enn __ ve) 


min mn 
Rememberine that’ é@ =1 and 2" = cosimn = 
(— 1)”, we get, on equating real and imaginary parts, 
b } b 
so (11)? 1; and “b, = —{— 1)?=?) 
Care 2 ( ) n me ) 


a . . 
The constant term < is obviously b so we have 
4 


fa) = aft — 2 (4 cos 3f , Cos St q 


I 3s ss 
ra rant fest 27 Hn sin se oe) 
2 i 


Hitherto our artificial functions ae been com- 
posed of straight lines or parabolas; we will now 
take some cases where they are composed of segments 
of sine curves. In the next two problems it will 
be convenient to take the period a sub-multiple 
of 27. 

E 
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Problem IX.—Analyse the rectified sine curve. 
The function consists of a repetition of the positive 
halves of the sine curve y= Asint. See Fig. 15. 

Here the period is x and the function is an even 
function, so that 


fd) = E+ ay cos 2t + ay cos 4t + : 


Tv 
2 : 
where an = — \ h sin t cos nt dt. 
MEA 
() 


Fic. 15. 


But the curve / sin ¢ is symmetrical about an axis 


of y through ¢ = = and so also is cos mt when 2 is 


: 14 
even; hence we may take the integral from o to ~, 


giving 
7/2 


ry) : 
an = a sin ¢ cos nt dt 
Tt Oo 


T/2 
rd.) : : 
= SS | {sin (n + 1)¢ — sin(n — 1)t\ de 
mer) lee (n— 1)t  cos(n+ wal ia 
1% n— I n+ tI 


fe) 


Since m is even, m — 1 and m+ 1 are both odd, and 
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therefore both terms vanish at the upper limit so that 
we have 


res 2h ie! } poe 40 
oa eT eee T a(n? — 1) 
This assumes no indeterminate form when »=0:; 
hence 
ft Ne _ COS 2t cos 4t cos 6t _ ie: He 
2 ie) en's Gay 
(11) 


The student of electrical engineering will recognise 
the rectified sine curve as the voltage wave form of a 


(© 
~“F% m7) Toe OP se w 5% 
Fic. 16. 


simple dynamo with a single coil in the armature and 
a two-piece commutator. The next problem will 
give the voltage wave form of an armature with three 
coils at 120° apart and a three sector commutator.* 


Problem X.—The positive halves of the sine wave 
y = Asing¢ are drawn at intervals of 22/3 and the 
ordinates added where they overlap. Analyse the 
resulting curve. 

This curve is shown in Fig. 16. 

In order to add the overlapping curves, we will 
draw the axis of y through a point of intersection as 
shown. The equation of the arc AB is then y = 


sin ( +) and of CD is y = Asin (" — ‘). The 


* The “field form” is supposed to be sinusoidal in both these 
problems. 
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sum of these two is y = 2b sin 7 cos t = h cos t, and 


and is the equation of 


. nt Sa 
this holds from — é {OL ; 


the arc EF. But the equation of the arc FG would 
be y = hcos t if the axis of y passed through its mid 
point; so that both arcs are similar and hence the 
required periodic function has a period of 7/3 or 27/6. 
The function is also an even function, so that, using 
the formule for the period an aliquot part of 2m 
on p. 29, we have 


ft) == + @, COS. Of -- Gy(COS 122 — ene 
12 7/6 
where an = — \ Td) cos nt dt 
vt [@) 


12h (7/6 
as \ cos ¢t cos 6nt dt 


fe) 


6h 7/6 
ie \ {cos (6m — 1)t + cos (6m + 1)thdt 
Oo 


A 
: 7/6 
e o|2 (6n — 1)é is sin (6” + zd 
qr 6n — I én +1 


oO 


_ 6h{sin (nx — 2/6) i sin (mz + a) 
Gi | kOe aay Gee 


Both sine terms are numerically equal to sin 2/6 
Ora ty 2. oe both are of opposite sign, so they give 

m(36n — 1) 
nately + and — starting with a, + ; 4, is also obtained 
by putting m =o in this expression. 

Hence we have 


an numerically, but with terms alter- 
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f) = eae 4. 608 Orem cochl 25 me cos of ae ! 
BOY HW ee £7 19 
= h{1-955 + 0:055 cos 6f — 0-013 cos 12¢ 
-- @006 cos 18f — . . iA 


which shows that the amplitude of the lowest harmonic 
is less than 6° of the amplitude of the constant term. 


We will now take two mechanical problems leading 
to artificial functions composed of arcs of sine curves. 


Fic. 17. 


Problem XI.—A body performs simple harmonic 
motion in a fluid in which the resistance is proportional 
to the square of the velocity. Analyse the periodic 
function representing the resistance to the motion. 

If the velocity is represented by & sin z, the numerical 
value of the resistance will be given by Asin? z; but 
since the resistance always opposes the motion, it will 
change sign whenever the velocity reverses, so the 
resistance will be a Class IV function represented 
graphically by Fig. 17 where also the thin line shows 
the first harmonic. 

Here f(z) is represented by / sin* ¢ from t = 0 toa; 
hence 


Ha) = 6, (sinit 4 Oy si. 38a. 
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/ 
AE : 
where bn, = ve sin? ¢ sin nt dt 
7 
[e) 


7/2 
= 2A (1 — cos 2¢) sin nt dt. 


La: 


The first term of this integral gives 


7/2 
2h E a ; 
Eo. i, e) 
aL—nad, 
which, since m is necessarily odd, is equal to 2h/xn. 
The second term gives 


h n/2 . 
ae | {sin (n + 2)¢ + sin (n — 2)t}dt 


1/2 


= ae (n + 2)t 4. 608 (n 4 


A n-+ 2 n— 2 


, : Yeh h I 
Remembering that is odd, this gives — = (- rs 


ss +) 2 Soe) Cres hiae heme ecae amen 
n— 2 a(n® — 4) 
_ — 8h 
first term we get by = pa Ph 


fl) = ie ¢_ (Sin 32s), sitgt Sesin Fie } 
Le ES, 1-3-5 353 /p anor? 
The loss of symmetry in this expression due to all 
the terms being negative save the first is in appearance 


only: really the factors in the first denominator are 
— I, I and 3. 


and so 


Problem XII.—A reservoir is divided into two by a 
vertical partition with a small hole in it near the 
bottom. The height of liquid in one part is kept 
constant and may be taken as zero; the height in the 
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other part varies according to kcost: analyse the 
function representing the velocity of flow through 
the hole. 

The velocity of flow is proportional to the square 
root of the numerical difference of heights and changes 
sign whenever this changes. So we have to analyse 
an even Class IV function of period 2% which is equal 


to 4./cost from t= 0 to t= = This function is 


shown by the heavy line in Fig. 18,* the thin line 
representing the first harmonic. 


Fic. 18. 


lero) (i), COS. f y- a5 C08 32 


4h m/2 
where an = 2) cos 4t cos nt dt. 
u% Jo 


A Reduction Formula may be found for this integral 
which will make a, depend upon ay-2, and thus finally 


? d : : 
upon 4,. Consider 7008 3/2¢ sin nt), where m is even. 


It is equal to 


— 3 cos # sin ¢ sin nt + n cos 3/2¢ cos nt 
2 


* By an oversight it is represented as an odd function in the 
diagram. 


56 A PRACTICAL TREATISE ON 
3 


= COS in| cos t cos nt — 5 sin ¢ sin nt 


MC Ose) 


Sam 


{en + 3) cos(n + 1)t + (2n— 3) cos(m — I)éy. 


Now integrate both sides of this from t= 0 tot = 


Nia 


7/2 
The left-hand side becomes [cos 3/2¢ sin nt] which 
(@) 


is zero since it vanishes at both limits. Hence we 
have 


m/2 \ ] —~3 n/2 ; 1 
t t = — — 
\. cos 4¢ cos(m+ I)tt ont 3 cos #¢ cos (n—1)¢ dt. 
Putting »=2 we have a; = — a1 ; putting 
1h == bt i = — ae = Z re a 
4 we obtain 4, - as he ane a,, and 
so on. 


H = { aed 1 3 
ence: (Zia, cocee ae: 3t+4 : F008 SE 


oe eee 
rates = | 
7/2 
where 2+ { cos 3/2¢ dt. 
T ie} 


This integral can easily be evaluated by Simpson’s 
rule, taking about 5 intermediate ordinates. It will 
be found that Gy =e Ise 


; The student familiar with the Gamma Function will know that 
iy Vr 15/4) a 

if cos3/? @d@ = 2° BG/4) (See Williamson’s Integral Calculus, 
Checter VI and the Table of log I'(p) at the end of the chapter.) We 


may point out that when 1 is large an approximates to ++ ) a 
77H 
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If we integrate the expression for f(#) we shall 
obtain a result proportional to the total flow through 
the hole up to the time ¢, viz. 


1-113h{sin ¢ — sin Sire Sigs ces Sl 7e our \ 
a: oy 39) 

A General Formula for the Fourier Constants 
in terms of the Discontinuities. We shall now 
obtain a very general formula for the Fourier con- 
stants from which the expansions in nearly all the 
previous examples can be written down straight 
away. Let f(t) be a periodic function of period 27 
represented in different portions of its range by 
straight lines or by arcs of parabolas of any finite 
deeree,™ 

Let the function, or any of its differential coefficients, 
have discontinuities of magnitudes Ta, Is, Ty . 
V’., I’g lic elacue ee, at ste, DOINTS 0-8.) Vn.*. 
using the notation of Chapter I, p. 14. We will 
suppose, to save extra words in the course of the 
proof, that there are no discontinuities at ¢= 0; 
although, since the origin is arbitrary, it is obvious 
that such merely require treating like any of the 
other discontinuities. 


We have f(s) = “2+ a cost+...+,sint+... 


where 
I 2 
ayn + thn = —\ fee dt. 
wt Oo 


Integrating by parts we get 


I ie |e z\ fem dt... (12) 


aL in 1 sain 


* A parabola of the pth degree is the curve given by the equation 


y= 6, ott cgtt . . . + epe?. 


58 A PRACTICAL TREATISE ON 


Let now f(z) be represented by one analytic function 
from o to a and by another from a to 2a. The first 
term when taken between limits gives 


AL fa — ope — fone + flameitar — fla + oe) 


By supposition there is no discontinuity at o and so 
fax) = flo); also ¢2*” = ¢i™°, so we are left with 


oS etna ‘ 
| fla + 0) — fla — 0) 

that is, with - Hees 
Nn 


If there are also other discontinuities at f, y, etc., 
it is clear, by breaking up the range of integration 
into ranges from 0 to a, a to , etc., that the first 


. 1 . 
term of (12) will become —2I,e"*, where the 
7M 


summation covers all the discontinuities of f(#) in 
magnitude. 

If we similarly integrate the second term of (12) 
by parts we get 


= J sae I Py . 
ce ae ae Sa ge t eint dt . . I 
Nn [ 1n ze 7121? i ( ) ( 3) 
In an exactly similar manner the value of the first 
term of this expression is 


ee ES [tie 
ane * 


where I’, is the discontinuity of f(#) at t= a and 
the summation extends to all such discontinuities. 
Similarly, if we integrate the second term of (13) by 


parts we get 
a0 int 
: E (Ye |+ : | Pein dt, 


mn? in min? 
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The value of the first term of this is similarly 


1 ’ 
— a gles 
un 


Proceeding in this manner, we get the final result 
that 


n(an aE 1by) pe —~L1,e" — —Z1'aei" roe “1” geina 
nN nN nN 


+ GET“ acite pee eld) 


This series terminates; for if p is the highest 
degree of any of the parabolic arcs, f(¢+ (2) is every- 
where zero, so there are no discontinuities, 1?** or 
any higher ones. 


Separating real and imaginary parts we get 


I : I I ? 
Nan = — wale sin na — wae le cos ma + ele sin na 


I ttt 
a acl a cos 7a — ° ° . (15) 
and 


I ; I 
th, = SAL eos 7a — => 14 sin 7a. —21". cos na 
n n* nN 


ee ein ma Fe (16) 
nN 


As an example, we will apply this last equation to 
the function in Problem I. Here we have a discon- 
tinuity + 4 at —a and one of —4/ at +a; and 
there are no discontinuities in any of the differential 
coefficients, these all being zero throughout. Hence 


: ‘ Be, 
Fat) sin) Ze h) sin na = 2? sin na, 
n n n 


the same as previously found. 
Again, let us take the function in Problem III. 
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Here only f(#) has any discontinuities; and has a 
; bot 


7 — a a(z — a) 


discontinuity of magnitude aie 


at +a. 


at — a and one of magnitude : = 


ha 
(a — 4) 
Hence by (16), 


ion & 

n? a(x — a) 
ee 2h sin na 
ae a(t — a) "42 


I ha 


sin na 
n? a(t — a) 


sin (— na) + 


ntby = — 


or 


again in agreement with the value previously found. 
Formule (15) and (16) can be simplified in the case 
of functions belonging to the four special classes by 
consideration of the relation of the discontinuity at 
—atothat at + a: we leave this to the student. 

The formule (15) and (16) may be employed even 
when they do not terminate, provided that they are 
convergent. As an example, we will consider the 
rectified sine curve dealt with in Problem IX, save 
that we will now make the period 22 instead of a. 
The function is then represented by A sin ¢/2 from 
t=0O to t= 2n. 

The function itself has no discontinuities in mag- 
nitude, but its differential coefficient has one of mag- 


, 2h ’ 
nitude oy atz= 0; so that |= A. * Its second order 


differential coefficient has no discontinuities in mag- 


nitude, its third has one of magnitude — x at 
~=0; Nie) that le = =o similarly, ib aoe Ly 5 he Nase 


16” 
= iy etc., whence we get by (15) 
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nde =) cia Ge 
HOLE) Se 
aye n®\16 


= I 
a ~ 20: Soe 
b 4h 


=) I n= — 1 
4n 


, 2 
Since the constant term is 2h we get 
Ls 


- lh Cost cos 27 | cos 3? | 
joss = = ie 
KY) am \2 Tens 3.5 Say. 


also in agreement with the previous result. 


EXAMPLES. 

1. Prove that the odd function of period 2m which is equal to 4 

TOMI — 10) tO f— 70 1S 
) , gin 37 Sl Ste Si 77 \ ; 
i — sin ¢ = SuaeA 
f= 2 (sine SE 4 Sy 

2. Prove that the even function of period 2x which is equal to ¢* 

from f= 0 to 7 is 


cos 2# cos 3¢ , Cos 4t \ 


A) = 4{% — cont hee 4 aes fa. 


3. (2) By putting ¢= 7 in Example 2 obtain the relation r= 


I I I 
7G tae 
(b) By putting t= im (6), p. 42, obtain the result = 


I 
Dy boats 


(c) By putting to in (10), p. 47, show that == I— zat 
an 

(d) By putting t= : in the result of Example 4 (below) show that 
ae 
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4. Prove that the odd Class IV function represented by /(t) = 
3 
o(¥— 2 from t= 0 to ~ is 
ue) 2 
96h. sin 3¢ , sin 5¢__ sin 7¢ Ae 
fo = 9B {sin 2 — ial Fc 7A 4 ty 
5. Four circles are drawn on the sides of a square as diameter. 
Prove that the polar equation of the quatrefoil defined by the four 
external semicircles is 


Pes ae cos 49 cos 80 , cos 120 \ 
jt Zeer AON Ait Lae ee 
where a is the side of the square. 


Tv 


6. Prove that the periodic functions of period 2x represented by 
sin mt and cos mt from — 7 to 7 are 

sin ¢ 2 sin 2t 3 sin 3t \ 
Pane Ae ae ea 


. PF hs 
sin mt = — sin m{ 
Tw 


and 
m cos t m COs 2t 1 


= see 


9 


Zan [fu 
cos mt = — sin mn — —_—- 
a Ea 12— m2 22 


respectively. 


7. A long uniform chain hangs over a number of small smooth pegs 
at a distance 2m apart in a horizontal line. Each segment. hangs in 
a catenary whose equation is y = A cosh x/A. 

Prove that the complete form of the chain is given by 


a gael ie 
Ale rf I+ 422 ererve ma 


Obtain this result, (2) by evaluating the integrals in the usual way 


aie mf cos cos 2¥ Cos 3x 


y= = sinh 


with the help of the integral fee cos nx dx; and (b) by (15), p. 59, 
from the discontinuities of the function at 7. 


8. Prove, by considering its discontinuities, that the odd Class IV 
function of period 2m which is represented from 0 to « by the line 
joining the points (0, 0) and (a, 4); is represented from « to B 


by the line joining the points («, 4) and (8, &), and from 8 to . by a 
line parallel to the axis, is given by 


S(t) = 4, sin t+ bg sin 3¢ + bg sin st+ ... 


t= ACfED ane =) se} 


where 
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g. Obtain the values of the coefficients in Ex. 8 when « = A g=*, 
. 3 
I 
Dis a and) kie= “4 so that the points lie on a sine wave. 


to. Prove that the even Class IV function of period 2x which is 
equal to cos ?¢ from t = 0 to t= 7/2 is given by 


= Shas oo ead de PEM 
() = C{cos # AIS Le ESE See ok 
pe ho eae 
Ps ppc oe Pp 
m/2 
where c=t{" cos Pt tp dz,* 
wT 


oO 


cos 7#-+ . 


——— 


pani ‘ : aul 
[Consider 7, (cos? + ¥y sin nt) in a similar manner to that adopted 


in Problem XII.] 

Note that when ¢ is odd this series terminates; also that the given 
function is then represented by cos ?¢ throughout its whole range so 
that we have expressed cos ?¢ when is odd in a Fourier Series. 


11. A periodic function of period 27 has no discontinuities of any 
kind except at t== 7; that is, it is represented by one and the same 
analytic function from — 7 to-+ 7. Prove that its Fourier constants 
are given by 

ioe pa ii 
ioe CoS as r+(- ea Fea Die er ye 


and 


I I" wv 
thn = (— Bier (—1)*— oa (— ate A eae 
where 10") = f— n) — fe). 


Obtain the expansion when f(t) = ce*/A. 


* The value of this integral in terms of Gamma Functions is 


ao 
a, es (Williamson’s Integral Calculus, Chap. V1.) 
ae: 
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12. Prove that the odd Class IV function of period 2x represented 


Tv. 
by 4 sin ct from t= 0 to t= 18 


ch crf sin t sin 3¢ sin St } 1 
fi) = = cos ae a eye 


2li— 4. og 25 


By making c tend to zero while he remains finite, obtain the expan- 


Tt 
sion when f(t) = ¢t from t = 0 to 2 


13. Prove that the harmonic analysis of the even function which 


Tv 715 
is equal to hcost from to to ia and to zero from i= to 


a= IS 
cos 2t cos 4¢ , cos 6t \ 


Aplin . ss bE £ 
fi) = p+ 7 cos tp jae pe 


all odd harmonics being absent except the first. 


14. The positive halves of the curve / cos ¢ are repeated at intervals 
ofee Prove, by the result of the previous example, that the harmonic 
n 


analysis of the sum of these curves is 


zoe cos 2mt cos 4nt 
ales (2n—1)(2n+ 1) (4n— 1)(4n-+ 1) 
cos 6nt 7 
7 (n= 1\(oncd- 1) nea eae 
if 2 is odd and is 
2nh| 1 n[z—1f cos nt __ Cos ant 
aE eee: \@—D@tD @nx—Dean+n 
cos 3nt = i 
pcre eens 
if m is even. 


Hence note that the voltage from a commutator with 27+ 1 


sectors is exactly the same as that from one with twice the number 
of sectors. 


<r : : ™ 
15. Two sets of equidistant ordinates at distances of — apart are 
m 
drawn, the first set commencing at =o and the second set at t-= 
Tw : 
om °° that the members of the latter are midway between those of 


the former. Through all the points of intersection of the first set 
with the curve f(?) = sint short lines are drawn, parallel to the 
axis of t, terminating at the adjacent ordinates of the second set. 
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Prove that the harmonic analysis of the step-like sine function 
defined in this manner is 


od 2 ae ee sin (2m —1)t sin (2m-+ 1)t 
ee 
sin (4m— 1)t , sin (4m- 1)t 

bee aN oe eg 

16. A body performs oscillations in a medium in which the resist- 


ance is proportional to the square of the velocity so that the equation 
of motion is 


at aC 2) + ly =o, 


the sign changing every half oscillation. 

Prove, from the result of Problem XI, that, if the damping is so 
small that its effect on the period may be neglected, the motion is 
given by 

gah sin (nt + y) . 
ee m j 


where A, and y are arbitrary constants. How do you explain the 
fact that when ¢ is large this result is approximately independent of 
the initial amplitude A,? 


17. Show that the result of Ex. 16 may be written 
ee eee : 
le aaa TS 


where — ?, is the time at which the amplitude would have been infinite 
if the result was valid for large amplitudes. 


CHAPTER IV 


ON THE REPRESENTATION OF PERIODIC ANALYTIC 
FUNCTIONS BY FOURIER’S SERIES 


Ir is clear that any algebraic function of sines or 
cosines of @ and/or integral multiples of 6 must, since 
each of the sines and cosines is periodic in the interval 
2a, likewise be periodic in that interval, and therefore 
must be expressible in a Fourier’s series. ‘The Fourier 
constants will, of course, be given by the integrals in 
the usual manner; but they may be obtained in a 
much simpler way in these cases, for the Fourier 
expansion is nothing more than another representation 
of the same function which is already defined for all 
values of @ by the single given expression. It must 
be possible, therefore, to transform the given expres- 
sion into the Fourier series by mere algebraical manipu- 
lation. ‘Thus instead of obtaining the expansion by 
evaluating the integrals we may use the expansion to 
obtain the values of the integrals. 

To effect this algebraical transformation, we replace 
the sines and cosines by their exponential values ; 
then, for brevity, denote e’® by a single letter, say w, 
and then arrange the expression in positive and 
negative powers of u, grouping u” with u-”. We 
finally replace u” + u-” by 2 cos 10 and u” — u-™ by 
21 sin 78. A few examples will make the method clear. 


Problem 'I.—Express, in the form of a Fourier’s 
series, the inverse square of the distance between the 
two points given by (7,6) and (4,0) in Polar Co-ordinates. 

66 
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The distance squared is a? — 2arcos@ + 72, so 
that we require to represent 
I ue I 
a* —2arcos8+ 7% a? —aru+tu-)+r 
e u 
(a — ur)(au — r) 
in ascending and descending powers of u. By the 
method of Partial Fractions, the expression can be 
written in the form 
gear 


2 Awe 


where P and Q are independent of u. 
Adding these fractions and equating to the given 
expression we get 
Qa— Pr=o 
and Page O7 2s 1: 
Solving these for P and Q we get for the required 


expression 
I ( ae ! Tr ) 
2 Nae a 


If r is less than a this may be written 


I I r I 
2 ( 1 ; Jk 
a@—r\1—urla au 1—r/au 


the expansion of which is 


I r NAS (1) 
aot + Get G)* 18 eg Bs deka 
2 
+ (C)u- + () Uae sale Ry Ue a 
a a 
== a p-{1 + 2(Z) cos@ + 2(7)° cos 20 
Aree -5T a 
“N32 
+ 2(7) cos 30+... fe, 
a 
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which is the Fourier series required. It is convergent, 
since r/a is less than unity. Ifr/a > 1, we can obtain 
in a similar manner the expansion 


{x +2(%) cos 6 + 2(2) cos 20 + pa: \. 


Problem II.—To obtain the Fourier series repre- 
senting cos” @ and sin” 6. 

These expansions take a different form according 
to whether 7 is even or odd. 


We have 
ann ¢ ae yg 

Z 

ae u2” +. 2ny2"-2 2n(2n — 1) 2n-4 

a, nu + —~—___‘4 +. 
_— 

4- \20 Tie + | u° 20 SS) a 
ites || n |e 4-0 (ait 


= Rae 


‘|2n {3 . cos 20 + mn — 1) cos 40 


~ 28"-I lal n n+I (n + 1)(n + 2) 
betta 1)( 2) | 
Cie NCEE BD COs 00.2. le 
In a similar manner we have 
2n-- I § — j2u +1 Ac ee 
cos 2" Inn +1 au Gog 30 


n(n — 1) deg aie nn — 1)(m — 2) 


(n+ 2)(n-+ 3) (w+ 2)(n+ 3)\m+4) °°” 


Basal 
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and likewise 
| 21 {2 n n(n — 1) 
8, 0 ais 
Zoe 2 Aiea ne Gs) 


cos 49 — ile Iie 2) Cos 
(n+ 1)(n + 2)(n + 3) ar -} 


sin2” §@ — 


and 
Soeeriee |2n-+ I {si es 
In ier. ane in]a ht sin @ Sorry sin 30 
oS mn — 1) Ene ss Des in 2) 
(+ 2) +3) > Get 2\(n + 3) + 4)" 
sin 70+... } 
In an exactly similar manner it may be proved that 


n(n 


Os" 0 COs 710 == mae I + ncos 20 + me coe 40 


a: ee 2) COS, Ob me: 21, } 
whether be even or odd; and in a similar manner 
like expressions can be obtained for cos” @ sin 10, 
sin” 8 cos 76 and sin” @ sin 78, these expressions taking 
different forms in the latter two cases according as 
m is odd or even. 


Problem III.—Let us represent the displacement of 
an engine piston in terms of the stroke 2a, the length 
of the connecting rod /, and the crank angle 6. The 
expression giving the piston displacement, x, is 

#=acos6+ //? — a*sin?g—1 . . (1) 

Let the ratio a/l be denoted by e 

Then we have, by the Binomial ‘Theorem, 


x = acosé + tr = I A2 sin? @ — T 4 sint 0 
0 Z 8 


— Lo sinde — 5 o8sinta — ... — 1}. 
icon 128° ° 
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If we replace the powers of the sines by their expan- 
sions as found in Problem II, namely, 
2 sin? @ = I — cos 20, 
8 sint@ = 3 — 4.cos 20 + cos 48, 
32 sin’ @ = 10 — 15 cos 20 + 6 cos 49 — cos 60, etc., 
and arrange the result in the form of a Fourier’s series, 
we get 


Zane ( eS Cane Se ) 6 
ee Comer ee = COS 


a (“ at OS 150° it .) cos 26 
Asma © mies hz 


oe (70 se ) 10 (es ) 60 
(6, + 85 + ---Jeosaa +(5, + - --)oos 


This expression is chiefly required in practice to 
obtain the acceleration of the piston. If @ = nt, so 
that mis the angular velocity of the crank, we have 


3 5 
a BAN cos @ — (c a g = Ee: se cos 20 
| Fi Pres: 


o° , 30° ) = (2 ) 
+( + + ...)cos 40 Sp EE oy 


The numerical maximum of this expression is when 
8 =o and will be found to be m?a(1 + e@) which can, 
of course, be obtained by differentiating (1), or the 
Binomial expansion, twice. 


Tue Direct SUMMATION OF FourIER’s SERIES. 


Closely allied with the preceding representation of 
an analytic periodic function as a Fourier series is the 
direct summing of a given Fourier series, which is 
simply the inverse process. How to tell whether a 
given Fourier series represents an analytic function 
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or an artificial one will be dealt with in Chapter V1) 
A single example will suffice to illustrate the method 
adopted in either case. 


Example I—Find the sum of the two series 


c? c3 
Rp wagin ya s-as tor CO8 30, ee, 


‘ ce. Gens 
and be at 


where ay ees 

If we denote the sums of these series by C and § 
respectively, we have, on multiplying the second one 
by z and adding it to the first, 

2 3 
C+ 18 = Lente aera oh Oe 

which is easily identified as — log(1 — ce’), 

Taking exponentials of both sides after changing 
their signs, we have 

e~“cos § — zsin S$) = 1 — ce’, 


whence 
e-S&cos§ = 1—ccosé@ and e-©sinS = csin# 
c sin 6 
and hence tan’ S =( 
I — ¢cos@ 
and 


e-2((cos?S + sin?S) = ¢~2© = 1 — 2¢ cos 8 + ¢? 
and so 


C = — #log(1 — 2¢ cos @ + c*) 
and eo ae 
I — ccosé 
If c= 1, these simplify to 
C = log(4 cosec 4/2) (from 6 = 0 tom) 
and Ses ys (from 6 = 0 to a), 
oe se2 


the latter agreeing with the result in (7) p. 42. 
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EXAMPLES. 


1. The length of the focal radius vector of an ellipse of eccentricity 
¢ and semi-latus-rectum / is given by 


1 

T= ————_.- 
I — e cos 

Prove, by comparing this expression with that in Problem I, that 
the Fourier’s series for r is 
r =Isec B{1 + 2 tan B/2 cos @ + 2 tan? B/2 cos 20 

+ 2 tan? 6/2 cos 30+ . 

where sin =e. Note that if a and d are the semi major and minor 
axes respectively, we have b= a cos B and, since / = 5?/a, the factor 


lsec B is simply equal to 5. 


2. Prove that the sums of the series 


c c 3 
1+ — cos 6-+ -— cos 20+ — cos30-+ ... 


E 2 [3 
and 
2 
cohen Oe sin 26 ++ Gar 30+... 
[2 2 EB 
are e° £089 cog (c sin 8) 
and e€ £089 sin (c sin 8) 
respectively. 


3. Show that the sums of the series 


cos 6 + cos 30 + ee : 
and sin 9 + # sin 30 + oe Ns sin 50+ . 
are cy 
2" cos" 8 cos (n + 1)0 
and 2” cos" 8 sin (n + 1)0 
respectively. 


4. Show that the sums of the series 
1-+ ¢cos0-+ ¢? cos 20+ c?cos304 ... 
and esinO+ c*sin 20+ c3sin 30+ . 
where | c| is less than unity are 
I — ccos@ c sin 0 
I — 2¢cos0 + c? I — 2¢ cos6-+ c? 
respectively. 
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5. Obtain the results : 


cos 94 SE 4 SEY ~. « = log | +4 cosec 6/2 | 
cos 20 , cos 30 
BOS et oe. . = — log | $ sec 6/2 | 
cos + SEI 4 oe .« « = Flog | cot 6/2 | 
sin § — = + eS . . . =F log | tan (0/2 + 2/4) | 


= flog|sec@-+ tan 0| 
These results should be compared with the series 
sin 20 , sin 30 


sin 6 +- —- ; eee 
and 

: sin 20 _ sin 30 

sin § — . + 3 


which represent the function in Fig. 12, with the origins at the points 
marked o and x respectively ; and also with the series 


pee ee ee 
3 5 
and 
cos 30 , cos 50 
cos § — eee 
5 3 5 


which represent the function in Fig. 9(a), where the origin is taken 
at the point A, and midway between A and D, respectively. 


CHAPTER V 


ON THE REPRESENTATION OF EMPIRICAL FUNCTIONS 
BY FOURIER’S SERIES 


Tue empirical function which we desire to represent 
by a Fourier’s series is necessarily given us either: 
(a) by some experimentally obtained graph; or (2) by 
some experimentally obtained set of points. In either 
case, the ideal results that were aimed at have only 
been reached to a certain degree of approximation. 
This limited accuracy of experimental results should 
always be remembered when subjecting them to mathe- 
matical operations. 

When a graph is given, two courses are open to us: 
we may either apply some form of mechanical con- 
trivance called a Harmonic Analyser to the graph to 
evaluate the Fourier Constants, or we may select a 
finite number of points on the graph and proceed as 
if these points were all that we knew about it. 

Several mechanical contrivances for the harmonic 
analysis of graphs have been invented which possess 
very various limitations with respect to the number of 
harmonics they are capable of dealing with; some- 
times, indeed, too severe to make the apparatus of 
much practical value. Space prohibits us from dealing 
with these machines here, since nothing but a com- 
plete explanatory description with diagrams and illus- 
trations would be of any value. A useful account of 
several such machines will be found in the Article on 
Calculating Machines in the 11th Edition of the 
Encyclopedia Britannica. References to the original 

74 


FOURIER’S THEOREM as 


description of several of these machines will be found 
in this Article. 

These machines are all both complicated and expen- 
sive, and, apart from the mechanical errors of the 
machine (and sometimes theoretical errors too), the 
accuracy of the results depends upon the accuracy 
with which one can move the pointer of the machine 
along the given graph. We do not consider these 
machines essential for practical use ; and a new method 
of harmonic analysis which is described below, we > 
think, renders them all the more dispensable except 
perhaps where a very large number of similar graphs 
have to be analysed for the same harmonic. 

We will suppose, therefore, in the remainder of this 
chapter, that we are given a set of points on a periodic 
function, which we will suppose to be situated on 
equidistant ordinates, and that we require to obtain 
a Fourier’s series which will pass through the given 
points. ‘There are two distinct ways in which this 
problem may be solved. In the first method we find 
the Fourier’s series with the least number of the 
lowest harmonics which will give a function passing 
through the given points. This is the usual method 
of harmonic analysis and will be explained in detail. 

The new method referred to, which has not, as 
far as we know, been published before, is quite distinct 
and consists in joining the given points by straight 
lines (or by arcs of parabolas of the second or third 
degree) and then obtaining the Fourier’s series repre- 
senting this built-up curve from the discontinuities 
of the function defined by this series of arcs. ‘This 
method is somewhat more laborious than the first 
one, but it has the advantage that we know exactly 
what the resulting series represents. ‘There is, how- 
ever, no utility in applying this method when all 
that is known is a set of ordinates of the given 
function; for we could not then tell which of two 
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different curves passing through the given points 
represented the function most accurately between 
the points. But the method is valuable when a graph 
of the function is given if we take ordinates at such 
a distance apart that it is obvious to the eye that the 
form of the curve between two ordinates is very 
approximately parabolic. In this way we may obtain 
a Fourier’s series which represents the given graph 
with much greater accuracy than does the simplest 
Fourier series passing through the given points. Since 
points can be measured on a graph much more accur- 
ately than the pointer of an instrument can be moved 
over the graph by hand; and since, also, with a reason- 
able number of ordinates, the various parabolic arcs 
will agree more accurately (through their entire length) 
with the given graph than a pointer moved over it 
by hand however carefully one does it; we consider 
this method of harmonic analysis as being more 
accurate than any Harmonic Analyser. 


The Simplest Fourter’s Sertes passing through 2p points 
per period. 


We will suppose that the period is 2% as usual and 
that we require to find the curve with the lowest 
possible harmonics which passes through the 29 points 


(5,4), (Ba). Ae. 


We are at liberty to assume any series with 2 
unknown coefficients. If we assume 


f(t) = “2 + a, cos t + ag cos 2¢ + cin Tr iap COsmys 


+ b, sin¢g + bygsin2i+ ... + bpsin 9%, 


we have 2p + 1 unknowns; more than can be deter- 


mined, but we note that when ¢ = Ne pt = Na and 
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30 the last sine term is seen to be zero at each ordinate : 
it therefore makes no difference to the curve passing 
through the given points, so we omit it. We will 
assume, then, the expression 


(i) = “Say cost +... + aps cos (p— 1)t-+ *P cos pt 
+ b,sint+ ... + bp-r1sin(p—1)t (1) 


a ae 
= and - instead of a, and ap 
is, as will be seen later, so that all the a’s can be obtained 
from the same expression. 


the reason for writing 


If in this expression we put t=~, rik 
successively the left-hand side becomes A,, Ay, .. . 
Aop respectively by hypothesis. If we then add the 2p 
equations so obtained all terms cancel save the first 


by (F) and (G), p. 9, and so we have 
Pa, od ZAg. 


If we again write t=~, 2 Ran successively 
in (1); but before adding the equations multiply them 


2nT 


NIC ZPnmn . 
ByaCOce-.4COs =, 5... 4 COS ade respectively, where 


n> p and tee add, we can see by equations (H) 
to (J), p- 10, that all the cosine terms will vanish 
save an, and that all the sine terms will vanish also ; 
and hence we have 


an& cos? 177 — Aq cos ky 
p ig 


i] 
By the use of the formula cos? 6 = ime a) we 


easily see that the value of the sum on the left-hand 
side is p, since there are 2p terms, except in the case 


7 
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when n = p or n = 0, when the sum is 29, but since 
ae An 

we have allowed for this by writing ie and = instead 

of a, and ag, we see that 


I qu 
tn = =a NG 00s 2—— ee 
ee a 


in all cases. 
In an exactly similar manner, we find 


bye A sine ee 
ee 
We have therefore solved the given problem. 


We will now prove that, as the number of ordinates 
tends to infinity, the constants a, and by, tend to the 
Fourier constants a, and Dy. 


Aj is the ordinate at t = 2: we may denote it by 
p 
; n 
A; instead. In this case the factor cos ”” must be 
denoted by cos nt. The factor . is very small, when 


p is large, and is = of the distance between successive 
ordinates; but this distance, in the limit, is dz, and 
hence + becomes Z in the limit when the summation 
becomes an integration. The summation extends 
from t=” to t= 2a which, in the limit, is from 


o to 2x. Hence 
27 
I no I 
Beg COs v= becomes 2 A: cos nt dt, 
Pp P uM Jo 
which is the value of the Fourier constant a,. Simi- 
larly, b, becomes dy. 
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The expressions in (2) and (3) which consist of 2p 
terms can be reduced to one-quarter of the number of 
terms in the case of the four special classes of periodic 
functions. In what follows we shall denote the 
distance 2/p between two consecutive ordinates by a 
and by Ay we shall always mean the ordinate at 
t= qa. 

For an odd function, let us divide the half-wave 
into p parts and let the ordinates at 0, a, 24, . . . be 
o, B,, B,. . . By-1,0. The complete set of ordinates 
over a whole period is thus 


OME ebe re Dp as) ON — Does, 2 06. 1) By,,0. 
If these be inserted respectively for the A’s in the 


fommula Dy == ; & Ag sin gna, and if we note that the 


sine factors for B, and — B, are numerically equal 
but of opposite sign ; and that the sine factors for B, 
and By_; are also numerically equal, and of the same 
sign when x is odd and opposite when n is even, we see 
that 


b, = iB + (— 1)"-? Bp-r) sin na 
+ (3, + (— 1)"-1 Bp-2) sin 27a 


4- (3, + (— 1)" Bp-3) sin 3na + .. | 
(4) 


the series continuing till all the B’s have been taken 
once only. If there are an odd number of B’s, 1.2., 
if p is even, the last term in this expression will be 


Bpj2 sin”. The assumed series for the function is, 
2 
of course, 


AO OE ee as a 
Orr aan O 
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For an even function, let the ordinates be C,, CG, 


Cres tee Cpa at) 210,60, eer thatentie 
ordinates over a whole period may be denoted by 
CC. Criteesn Gham, Conn kG ree Inserting 


these for the A’s in ay = poe cos ngx, and grouping 


together, as before, the four terms which have the 
cosine factor of the same numerical magnitude, we get 


pee ; {fot = 1)"Cp ae (c + (- 1)"Cp-1) COS N& 


in (c, tyes 1"Cp-2) SEER es ‘ aera) 


the series continuing till all the C’s have been taken 
once only. ‘The series in this case is 


f(t) = “2 + ay cost + ay cos 2t + nee 


+ ap-r cos (p — 1) ¢ + *? cos pt oe sc ae 


For an odd harmonic function, let us take the origin 
at one of the points where the curve cuts the axis, 
which it necessarily does at intervals of z. Let the 
ordinates over a half period be denoted by 0, D,, D,, 
. . . Dp-1, 0, so that the complete set of ordinates 
over a period are given by 


op AB}, ID) ak. 5 Dp-1,0, — D,, — IB sas — ‘Dp-1. 


Substituting these for the A’s in (2) and (3) we 
easily get, on remembering that 7 is odd, 


an = 5 (Ps — Dp-1) cos ne + (D, — Dp-2) COs 2n% 


iy (D, x Dp-s) coe | 2 PERS) 
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and 


Da =7{(2 i Dp-r) sin na + (Dy + Dp-2) sin ana 


“ (D, i Dp-s) sin gnu +. . | Sea) 


the series being continued till all the D’s have been 

taken once only. It must be noted that these expres- 

sions are not to be worked out for even values of n. 
The series representing the function is 


Ji) = a, cost + a, cos 34+ a,cosst+ ... 
PD, sine bysin 37-4 (0) 


If p is even the cosines stop with the term 
ap-1 cos (p — 1)¢ and if pis odd with . cos pz, and the 


sine terms stop with bp-1 sin(p—1)t or bp-2sin(p—z)t 
respectively. It can easily be seen that the total 
number of terms in both cases is p. It may be noted 
that since f(0) = o the a’s satisfy the condition that 


aqacipeda teapot OS 


For a Class IV function, we need only measure 
ordinates over a quarter of a period. Let us divide 
the quarter period into p/2 parts, p being assumed 


. Ut 
even, by ordinates at — = a apart, as before. 
For an odd Class IV function, let the ordinates be 


n 
Pete pio at.O,- 0,20, <% * >: Then from 


(9) we have 
bs = aa sin na + 2E,sin2ma +... 


+ Epis sin | Pee ake (ll) 
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the factor 2 inside the brackets being associated with 
all the terms except the last. The series is 


fat) = b, sint + by sin 3¢ + . . . + bp-1sin(p — 1)f 
(12) 
For an even Class IV function, let p again be even 


and let the ordinates over a quarter period be F,, Fi, 
. . . Epjz-1, 0: ‘Then from (6), since now, 


C= = GG, = — F,; Cp-1 = — = — Bete. 


we have 
ae a + 2F, cos ma + 2F,cos2ma+ .. 
(13) 


the series being 

f= 4, cos t+ 4, cos 34+ 2. 2 4 apy cos(p ye 
(14) 

The student must remember that formule (11) and 

(13) are only to be worked out for odd values of n. 


There is still another great simplification which 
may be made in the practical calculation of these 
coefficients. Let us change ” into » — m in any of 
the expressions for ay or by. Since gna i.e. gna/p 
then becomes g(p—)x/p or ga — qnx we see that 
a cosine term will be unaltered if ¢ is even, and merely 
have its sign changed if g is odd, and vice versa with 
a sine term. So that the series for ap-n or bp-—n is 
always the same as that for ay or b» with the signs 
of alternate terms changed. Hence we add up the 
odd and even terms separately ; and the sum of these 
two sums gives an or by and the difference gives 
ap-n Or bp-n. In order to avail of this simplification 
in the case of odd harmonic functions, it is obvious 
that p, the number of ordinates per half-period, must 
be even. 
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Schedules for Harmonic Analysis. 


We are now in a position to construct Schedules 
for Harmonic Analysis. Such a Schedule for a Class I 
function when using 10 ordinates per half-period is 
shown on page 84. ‘The given ordinates are denoted 
by the y’s; both yo and 410 being necessarily zero. 
The sums and differences of these when arranged as 
shown are denoted by the U’s and V’s respectively. 
The first column contains the numbers from o to 5, 
and the second contains the value of the factor ~ sin x 
when p = 10, for the values of m from o to 5. Nine 
further columns for the nine harmonics are taken, 
arranged and headed as shown. We then multiply 
the U’s in turn by the successive sine factors and 
place the products, which are denoted by the w’s, in 
the columns for the Ist and gth harmonic; placing 
them alternately in the third and fourth columns as 
shown. If the sums of these columns by denoted by 
PpacdO, we have PO; = b,-and P, — 0, = b,. 

In filling in the w’s for the 3rd and 7th harmonics 
we place uw, in column § in line 3; wu, should go in 
column 6, line 6, but the sine factor of this line is the 
same as for line 4, so we place it in this line; “g should 
go in column §, line 9, but line 1 gives the correct 
sine factor, so it is placed there; «, which should be 
in line 12, has the same sine factor as line 2, with its 
_ sign changed, so it is placed as shown. ‘The directions 
for placing these w’s may be stated as follows: starting 
with u, in line 3 go up and down the lines from 0 to § 
in steps of 3 at a time placing the u’s alternately in 
columns 5 and 6, and changing the sign of the term 
each time the sine factor passes through zero. For 
the 5th harmonic, ™, v3, and u, occupy line 5, while 
u, and u, occupy line 0; but since the sine factor of 
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SCHEDULE III.—FOR ANALYSIS OF ODD HARMONIC OR CLASS III 
FUNCTIONS. (14 ORDINATES PER HALF-PERIOD.) 


Ordinates © Y,; Yo Ys Ya Ys Vo V2 
OS Viysu ae 11). 30) 9S 

Sumi anes CG We WA WRU, We Wa We 

Difference | O) \Vigi Via Wan) \eal Visue au Nie 


Harmonic. 
@| Sine | Cosine 
5 factor. | factor. 
Ist (& 13th). 3rd (& 11th). 5th (& gth). 7th. 
fo) fo) -1429 +, 
—v 
I] +0318 | -1392 |, v,|—4U; —v,;|—Us —U3 5 
2| :0620 | -1287 U5 | Ug Us| tg +U,|+U, 
3| 0891 | -L117 | us Us| My v,|—4; —vs, 
Zui ami? || AOR) U,| Ug —U,|—U, Us| ~ te 
5| °1287 | -0620 |u, v;| Us Us| My V4 
6| -1392 | -0318 U6 | Ue Us| Ug —U,|—U%, 
7| 1429 fe) Uy o|—u, o| U, o| uU,—tU, 
+U,—Uy, 
Sum of u’s Py (OF, 1B (O) PE 0; P, 
Sum of v’s Ry Si RG Ss R,; Sr R, 
P,+Q,=b, Ps; Q; = bg P; + Q; = bs P, = b, 
P,—Q:= bys} Ps — Qs = Dir Ps; — Qs = by 
R,+5,=a, | R,; +S; = 4s; R; + S; = a5 R, =a, 
R,—S,= 443} Rs; —S3 = ay; Rs; — Ss = ay 


The w’s are the U’s multiplied by the sine factor, and the v’s are the V’s 
multiplied by the cosine factor, of the line on which they are situated. 
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SCHEDULE III.—FOR ANALYSIS OF ODD HARMONIC FUNCTIONS. 
(14 ORDINATES PER HALF-PERIOD.) 


Ordinates o 104 192 264 320 360 384 392 
O 104 192 264 320 360 384 


SUM eee. 0 208 384 528 640 720 768 392 
Difference o o fo) fo) fe) fe) ow) 392 
Harmonic. 
3 Sine | Cosine 
= | factor. | factor. | 
4 ist (& 13th).| 3rd (& 11th). 5th (& gth). qth. 
fo) fo) "1429 
I | -0318 | :1392 6:6 ' | —22°9 —16:8 
2 °0020. || *1287 23°8 39°7 47:6 
BiOsOL | “biL7 | 47-0 18°5 —64°1 
4| 1117 | -o891 71°5 —85'8 42°9 
5 | +1287 | -0620 92°8 68:0 26:8 
6 | +1392 | :0318 106:9 53°4 — 89:2 
7 | *1429 fe) 56-0 — 56:0 56:0 29°7—70°4 
102:9—5'6 
Sum of “ sums.”’ | 202-4 | 202-2 7:6 as 19 1608} 
Sum of “ diffs.’’ 
by, = 4046 | bs = 14:9 by = 3-2 b, = 1:2 
by = i= Dir 33 Dyt=* 6 
a fo) a = to) a = 6) ar —10 
= fo) a fo) yy == © 
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TABLE OF SINE OR COSINE FACTORS FOR HARMONIC 
ANALYSIS. 


p or Number of Ordinates per Half-period. 


4 6 8 fe) 12 14 16 18 20 22 24 


fo) fo) fo) ° fo) fe) fe) 

353 | 167 | -096 |-0618 |:0431 |-0318 |-0244 |-0193 |-0156 |-0129 |-0109 
500 | -289 | +177 |:1176 |-0833 |-0620 |-0478 |-0380 |-0309 |-0256 |:0216 
333 | °231 [1618 |:1179 |-0891 |:0095 |-0555 |:0454 |-0378 |-0319 


*250 |*1902 |*1443 |*II17 |-0884 |:0714 |-0588 |-0492 |-0417 
*2000 |*1O10 |°1287 |*1039 |:085I |:0707 |"0595 |*0507 
*1667 |*1392 |*1155 |0962 |-0809 |:0687 |-0589 
+1429 |*1226 |*1044 |-089QI1 |:0765 |-0661 


+1250 |*1094 |-0951 |:0827 |-0722 
*LIII |-0988 |-0872 |-0770 
+1000 }-0900 |:0805 
-0909 |:0826 
0833 


this line is zero, they are omitted. ‘The even harmonics 
are calculated from the V’s, but the arrangement of 
the v’s, which are the V’s multiplied by the sine 
factor of the line on which they are situated, being 
exactly similar to that for the w’s, need not be further 
explained. 

In an exactly similar manner, a Schedule for the 
analysis of any odd function using any number of 
ordinates per half-period can easily be constructed. 
A table of the sine factors required for a Schedule for 
any even number of ordinates from 4 to 24 per half- 
period is given above. 

Page 85 shows Schedule I worked out for a par- 
ticular example. The equation of the resulting 
Fourier’s series there obtained is 


H(4) = 19°43 sin ¢+ 3:08 sin 22 — 1-43 sin 3¢ — 1-38 
sin 4¢ + +73 sin 6¢+ +37 sin 72 — +33 sin 8t — 49 sin 9¢. 

The given points will be found to lie on two straight 
lines and the Fourier analysis of this function is 
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2,200 | 2m sin't |, ... 4%. sin 2? | 
I(t) ae {sin er + ees ) + bay 


Or 


F(t) = 19°27 sin ¢+ 2-98 sin 2f— 1-32 sin 3¢— 1-20 sin 4t 
+°53 sin 6¢+-24 sin 7# —-18 sin 8—-24 sin 9f+ ... 

with which the Harmonic Analysis is in very good 

agreement especially for the lower harmonics. 

Page 86 gives the Schedule required for a Class II 
function, using 10 ordinates per half-period. It is 
constructed similar to Schedule I and is an embodiment 
of formula (6), just as Schedule I is an embodiment 
of formula (4). Page 87 shows the working out of a 
concrete example. It will be seen that the given 
points lie on two straight lines. ‘The Fourier analysis 
of this function 1s easily obtained from its discontin- 
uities. It has four discontinuities, of slope only, at 


o,a and + cA From this we easily find that 
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an = ree {- I — cos mm + 2 cos e 

from which we obtain 

f(t) = 17 + 6:26 cos t — 9:16 cos 2t — 1°82 cos 3¢ 
— :87 cos 4t + 81 cos §t — +39 cos 6t — +33 cos 7¢ 
— 57 cos 8 — -o8 cos gt +... 

as the Fourier analysis, in contrast with 

f() = 17 + 631 cos £ — 9:47 cos 2t — 1°96 cos 3¢ 
— I-99 cos 4¢ + I-00 cos §¢ —-53 cos 6t — +51 cos 7t 
—I-oI cos 8¢ + -16 cos gt 

as the simplest function passing through the given 

oints. 
Schedule III, on page 88, for the analysis of Class ITI 
functions, is also constructed in a very similar manner. 
This Schedule, as far as the sine harmonics go, is exactly 
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like Schedule I for the odd harmonics. A separate 


Schedule can, if desired, be constructed for the cosine 
harmonics, but these two Schedules may be put 
together, as it will be found that the terms in one 
Schedule occupy the spaces left vacant in the other. 
For this purpose, it is necessary to have a column of 
the numerical cosine factors adjacent to that of sine 
factors. Further, it is necessary that the v’s should 
be distinguished from the w’s in the same column so 
that both sets may be added up separately. It is best 
if all the v’s, and the cosine factors, are written in red 
ink; or, alternatively, the w’s kept to the right of each 
column and the v’s to the left. Care must be taken 
that the w’s are not multiplied by the sine factors 
by mistake, which is not likely if different coloured 
inks are used. 

Class IV functions present no difficulty whatever. 
If arranged as an odd function, they may be equally 
well worked on Schedule I or III and only half fill 
the Schedule as the “ differences” vanish. If arranged 
as an even function they may be equally well worked 
on Schedules II or ITI, and in this case the ‘‘ sums ” 
of the ordinates vanish. Page 8g illustrates an odd 
Class IV function worked on Schedule III. The 


; : : Sx 14" 
given ordinates lie on the curve f(#) = ——_—+ in—2), 
7 


from which it will easily be found that the Fourier 
analysis of the function is 


y a ; 
fy) = ees a sin ¢ ++ oe ie ae daa i) 


or 

F(4) = 404°6 sin ¢ + 15-0 sin 3¢ + 3-24 sin 5¢ 
+ 1-18 sin 7# +°56 sin of ++30 sin 11¢ +-18 sin 13¢ 
Sik: Mens 

with which the result 


FOURIER’S THEOREM 93 


f(4) = 404-6 sin ¢ + 149 sin 3f + 3-2 sin st 

+12 sin 7t + -6sin gf +-3 sin 11¢ + -2 sin 13¢ 
obtained by Harmonic Analysis isin excellent agreement. 
The agreement with the Fourier analysis is seen to be 
very much better than in the case of the examples 
worked on Schedules I and II. This is partly due to 
the fact that we have used 14 ordinates in place of 
10; but is chiefly because, with this function, the 
discontinuities are in its second differential coefficient 
and not its first. ‘The same function analysed by the 
use of 10 ordinates per half-period, but with the 
maximum ordinate-reduced to 250, gives 

f(4) = 258-0 sin ¢ + 95 sin 3¢ + 2-0 sin St 

++ +7 sin 7¢ + -2 sin Of 

against the Fourier analysis of 

F(t) = 258-0 sin ¢ + 9-60 sin 32 + 2°07 sin 52 

+75 sin7f +-35sing¢+ ... 

‘There remains now only the general periodic func- 
tion to be dealt with. This is easily analysed by 
resolving into its odd and even components. We 
measure ordinates over a range from —a to +2. 
If yg and y=, are any ordinates at ¢ and — # then 


Nt Meee eT is the corresponding ordinate of the even 
3 


component and J4 = Y-4 that of the odd component. 


These two functions are then analysed on Schedules 


II and I respectively. 
We frequently, in Applied Mathematics, when we 


desire to calculate a curve to pass among a number of 
experimentally obtained points, have a much larger 
number of points than disposable constants in the 
equation of the curve selected. In this case we cannot, 
in general, make the curve pass through all the points : 
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instead, we determine the constants so that the sum 
of the squares of the differences between the calculated 
and observed values of y,is a minimum for all possible 
variations in the constants. If the curve selected be 
represented by a Fourier’s series containing all har- 
monics up to the wth, when x is insufficient to make it 
pass through all the points, it may be proved that the 
values obtained for the coefficient of any harmonic in 
this way is the same as if we had assumed an equation 
with sufficient harmonics as in the first part of this 
chapter. The proof of this will present no difficulty 
to the student familiar with the Method of Least 
Squares. Hence in dealing with periodic functions 
the difficulty of determining the constants of an 
equation (when we have less constants than points) 
so as to represent the observations as closely as possible 
does not arise: we simply adopt the ordinary method 
of harmonic analysis for the lower harmonics and just 
neglect those harmonics which are higher than we 
want to consider. 


A Graphical Method of Harmonic Analysis. 


Harmonic analysis may be performed graphically 
by the aid of diagrams which save the necessity of 
constructing Schedules. ‘The accuracy in this case, 
though rather less, is often sufficient for practical 
purposes. Fig. 19 shows the method applied to 12 
ordinates per period. A large sheet of millimetre 
paper is taken and 6 lines inclined at 30° to one another 
are drawn as shown. 

These lines must be drawn as accurately as possible, 
their inclinations being set off by the tangents or 
cotangents of the angles, using as large a base as the 
paper permits, and not by a protractor. ‘Three 
concentric circles are then also drawn as shown. ‘The 
intersections of the radii with the smallest circle are 
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then numbered 0, 1, . . . 11 consecutively. The 
intersections with the intermediate circle are alter- 
nately numbered 0, 1, up to 11 and left blank. Finally, 
every third radius is numbered on the largest circle 
from o to 11 with the result that each numbered 
radius in this case bears three numbers. 


FIG. 19. 


The harmonic analysis is then performed as follows : 
We measure the ordinate A,, given either by figures 
or graph, along Ow and read its length in millimetres ; 
then we measure off A, along radius marked 1 on the 
small circle and read off its projections both on Ox 
and Oy; we likewise measure A, along the second 
radius to the same circle and read its projections, and 
so on up to A,,. The sum of all the projections 
along Ox divided by 6 gives a,, while the sum of all 
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the projections along Oy divided by 6 gives b,. If 
instead of adding the sums of the odd terms to those 
of the even ones, as we have just done to get a, and b,, 
we subtract them, we get a, and b,. 4 

If we measure A,, A,, Ap... along the radii 
marked 0, I, 2, . . . on the intermediate circle, and 
add up the projections and divide by 6 as before, 
we get a, and bg, while if we subtract the sum of the 
even terms from those of the odd terms we get a, 
and b,. Likewise a, and bg are obtained by using the 
radii to the largest circle, while a,/2 is one-twelfth 
the sum of all the ordinates. 


The Error of Harmonic Analysts. 


It is important to have some idea of the deviations 
of the a’s and b’s found by harmonic analysis from the 
true Fourier Constants the a’s and the b’s. We shall 
now show that the a’s and b’s can be expressed in terms 
of the a’s and b’s by means of the following 


Theorem. If 2 ordinates per period be taken at 
o, 7, ™ . . . and if each ordinate be multiplied by 
the sake of cos mt at that ordinate, the sum of the 
products (symbolically denoted a tN cos gnn|p) 

q=1 


is equal to 


plan + aop—n + Gopin t+ ayp—n t+ Qaptnt see } 

weer mcs LS) 
and if each ordinate be multiplied by the value of sin nt 
at that ordinate, the corresponding sum is 


Pion — bop—n + boptn — bypin + bypen— .. 2} 
an (16) 


FOURIER’S THEOREM 97 


For we have 
Ag cos grn|p— | 22 + a, cos qnu/p+ .. 


=> Dy sin gr} p+ ..3.. ! cos qux/p, 


and if the trigonometrical products be replaced by the 
sum or difference of two terms in the usual way, we see 
that the coefficient of by in the given sum is always 
zero, and the coefficient of ay is zero unless r — n or 
r + mis zero or a multiple of 2p, when it is equal to 
%.2p3; from which it easily follows that the required 
sum is equal to p(an + dap—n+ dapin+ ...). 

The second result follows in the same manner, 
but since the product of two sines is equal to the differ- 
ence of two cosines, the coefficients of the b’s are 
alternately + and — giving (bn — bop-n + bapin 

Since the left-hand side of (15) gives pan, we see at 
once that, in particular, 


a 

<O = 2+ aap + aap + aop + poi 

a; = 4,-+ A2p-1 + A2p+1 + Ayp-1 “ts Asyp+t a oe 

dg = dg + A2p-2 i A2p+2 a 

; Ss Ree re oe 1: f in 

Ap-1 = Ap-1+ Aptit 43p-1+ 43p+1 4 45p-1 1 --- 

ap = Ap + 4p + 43p + 43p + asp 4 
this last can be written 


“P= ap + asp + asp + arp + Jee 


Similarly, 
b, = b, — bap-it bop+i1— byp-1+ Dipti ie 
by = b, — bop-2+ bopt2— . lige 


bp-1 = Dp-1 — Dp + O3p-1 — O3p41 + AG, A 
H 
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If the a’s and the J’s diminish regularly and pro- 


gressively, we sce from these equations, that the 
accuracy of the a’s and b’s regularly diminish from 
a, and b, till ap-1 and bp-1; which latter two cannot 
be very accurate, since we have, roughly, 


aAp-1 = 4p-1 + ap+i and bp-x = bp-1 _ bp+t- 
The remarkable accuracy with which = represents dp 


should be noted. 

It should also be noted that in the equations (17) 
each a occurs once, and only once, and similarly with 
the J’s in (18). Hence each harmonic in the function 
analysed affects one, and one only, of the coefficients 
of harmonic analysis. For instance, if we are using 
10 ordinates per half-period, and if there is a large 
thirteenth harmonic while all the other harmonics 
beyond the ninth are small; all the a’s and b’s will 
be nearly correct except a, and b,; for we have in 
this case ay = ad, + 413 + dar + Gag + Gye + « < « 
and by = b,— dyg + Day — Ogg + Og — ss 


or, approximately, a, = a, + ad, and b, = db, — dys. 


Determination of the Amplitude of any Selected Harmonic. 


We sometimes desire to find the amplitude of some 
particular harmonic, say the mth, which we are interested 
in, without troubling at all about any other harmonic. 
If the harmonic is completely unknown, three ordinates 
per period of the harmonic is the least which will 
suffice to determine it, although if we only want either 
the sine or the cosine component two ordinates per 
period will suffice. With a larger number of ordinates 
per period, the harmonic can naturally be determined 
more exactly. We shall give formule for the cases 
of 3, 4, or 6 ordinates in the period of the harmonic 
respectively, which will at the same time indicate the 
error of determining a harmonic in this manner. 
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If in the formula ied) EAE 
2p-1 
~ Ag cos is = Pian + Gop-n + Gopin+.. (19) 


q=0 
we put 2p = 3n, so that the ordinates are at intervals 
of 2/3, we get 

ApeosA es 


JN ca its Bae _ 4 _ AAs ae 
cake a + As AG teal 


= BN ay an + dan gn + an + some Js 


and if in the formula 

SpE 

aa Ag sin as SS Pion = bop-n <p bop+n ee AOU Ay (20) 
q=0 


we make the same substitution, we have 


NSA, — Ay + Ag Ast Ay As + ¥ | 


= 3 1b — Ban + Ban — Dyn + bon — vee | 


If now the 2mth, 4mth, 5th ... harmonics are 
negligible compared with the mth, these formule give 
us approximately ay and dy, using three ordinates per 
period of harmonic. 

If in the formule (19) and (20) we put 2p = 4a, 
that is, if we use four ordinates in the period of the 
harmonic, we get 
A,—AgtA,—Agt - . - =20 {ant asntagnt+an+ ...} 
and 
A,—A,+A,—A,+ .. . =22{bn—Dgntbgn—dant ...} 
which give a, and b, provided the 3th and higher 
harmonics are negligible. 

If we put 2p = 6m, that is, use six ordinates per 
period of the mth harmonic, we get similarly 
2 a ae 
2 Z 2 2 


= 3n{an + dgn + gn + Ayn + Asn t+ - - ‘ 
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ane VSIA, + Ay An Ae Ae oe 


= 3n{bn Den ban Dak ae Dis i GF o * 


which formule give us a, and by, provided the 5th 
and higher harmonics are negligible. Similar formule 
for eight or any other number of ordinates per period 
can easily be written down. The simplification of 
these formule for the four special classes of periodic 
functions may be left to the student. 

As an example, suppose we wish to find the amplitude 
of the third harmonic in the odd Class IV function 


defined by f(t) = 18! fom t= 0 to t=" If we 
4 


use four ordinates per period of the third harmonic, 
that is twelve per period of the function, the ordinates 
over the first quarter of a period are 0, 3, 6, 9, so we 


easily find 
—6= 6{b, —b, +, ...4 or db; = —1 approx. 
Using six ordinates per period of the harmonic we 


find that the ordinates over half a period of the 
function are 0, 2, 4, 6, 8, 8, 6, 4, 2, o, from which we 


find 

a ( 8) = 9[0s — dig + bs mam | 
or by = — 9 V3 = — "70. 
Actually, 63 for this function is — ore SII. 


12 


A Simple Method of Harmonic Analysis. 


There is also a modification of this method of har- 
monic analysis, which we do not consider very satis- 
factory, but which may be employed when we are 
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given a graph of the function to be analysed. In 
this modification we draw different sets of equi- 
distant ordinates, and use one set of ordinates for 
determining one coefficient and another set for 
determining another, choosing the ordinates in each 
case so that we have a definite (small) number of 
ordinates to the period of each harmonic. The work 
involved in placing and carefully measuring such a 
large number of ordinates renders, we consider, this 
method more laborious than the one described, while 
the results are generally less accurate. The only 
saving is in the simplification of the sine and cosine 
factors—or their abolition in the case when only two 
ordinates per harmonic are chosen and placed at the 
positions of the maximum and minimum values of 
such harmonic. 

This method may conveniently be developed from 
the formule connecting the ordinates with the co- 
efficients. If, for instance, we have a general function 
and neglect all harmonics beyond the third and take 
12 ordinates per period it can easily be seen that we 
have the following equations for the coefficients: 


Wo = qy(A, + Ay + Apt «+ Aro + An) 
Ce ees (A, —A,+ A, Ps A, si Ag ae Ao) 
bs = ¢(Ay —A,+ A,—A,+A4,- 1) 
ag = HA, A; It Ag ae As) 
b, = HArs — Ags + Avs — Aros) 
a, = a3 >= Uae pee 6) 
b, adi bs ey 3(A; a As) 


Here Ay;-; means the ordinate half-way between 
A, and A, and similarly for the other terms in the same 
line. To use these formule, we thus have to place 
and measure sixteen different ordinates, while the same 
number of equally spaced ordinates would, with the 
ordinary method of harmonic analysis, enable us to 
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determine the first seven harmonics. All that has 
been saved is the trouble of multiplying the ordinates 
by the sine factors—a matter of very little trouble 
with a slide rule. 

For the general periodic function, a sufficiently 
accurate result for rough work can sometimes be 
obtained by assuming that all the harmonics beyond 
the third are negligible. Eight ordinates are then 
sufficient and it will be found that the ordinary 
formule for harmonic analysis give: 


8.22 = AFA + A+ A+ A + Ap+ Ast Ay 


4a; —= Al Al 707 Ry a see 1) 
4a, = A, + A,— A,— A, 
433°) As A707 At Ags tA) 


t= (A, + Aa + Ay + Ae) — (Ar + Ag + As + Ay) 


4b, = A, — A, + -707(A, + A, — A, — A,) 
4b, = A, + A, — A,— A, 
Abs = — Gia — Ay) JOT Ay Be he ee) 


A New Method of Harmonic Analysis. 


We will now explain a totally different method of 
harmonic analysis which, though it is rather more 
laborious than the ordinary method, possesses consider- 
able advantage from the point of view of accuracy 
and is quite easy to apply when it has become 
familiar. As mentioned above, the method consists 
in joining the tops of the given ordinates with straight 
lines or by parabolas of different degrees. This can 
be done in quite a variety of manners, but the three 
varieties of the method given probably include the 
most useful that can be obtained. 

Let us, as usual, take 2p equally spaced ordinates 
Yoo Nin » + + per period at o,2/p, . . . and for brevity 
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let x/p be again denoted by a. Let us join the tops 
of the ordinates by straight lines; then the slope 


between y, and ¥y, is 1 2e and that between y, 


and y, is”? = Yt; so that there isa discontinuity of slope 


at y, of 22— “v1 ms Yo. that is, of He ae Az, 
7 


is that second difference of the y’s which is vertically 
below y,. Substituting this value for I’, in the 
formule (15) and (16) on p. 59* we get 


m2? 


i ——— P sAry, cos mA 


(21) 
and by 


ayy? SA*y, sin 4 


ee 
B I 2 FYE 0 ; 
ut —2A*y, cos— is the an of ordinary harmonic 


analysis for the function whose ordinates are A*y,, 


el eae 
Hence we have 
pe 
Qi — aa oa | 
< eee See 22) 
and Dy — = P* on | 
mn” 


The quantities an and by can be calculated in the 
usual manner on one of the given Schedules so that, in 
labour, this method of harmonic analysis exceeds the 
ordinary method only by that involved in forming the 
second differences of the ordinates. "The Schedules, 


* The student must not confuse the present use of « with its use 
in these formule. In (15) and (16) « was used to specify the position 
of one of the discontinuities; now it is being used to indicate the 
width between successive ordinates; but no confusion need occur 
through this double use of this very convenient letter. 
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of course, only give an and by when n > p, but if we 
write p + in the expression for an we get 


I gp + n)r 
a n = XA? cos poy: TE fee 
Picts ? vq A 


which is the same as if we had written p — n, hence 
we have ap+n = ap—nand similarly bp4+n = — bp—n, 
and we see also that d2p+n = an and bapin = dn, 
hence we are able to write down the complete analysis 
of the function defined by the straight lines joining 
the tops of the ordinates. In actual practice, of course, 
the function to be analysed does not consist of these 
straight lines and so the harmonics beyond the pth 
are too inaccurate to trouble about. 

As a very simple example, let us take the function 
analysed in Schedule I, page 85. If we form the 
second differences of the ordinates we find them all 


27 27 
zero save at —, where we have — Io, and at — —, 


5 
where we have + Io. 
Hence we get 


10 . 200 ; 2nn 
bn = — —} — 10 sin —— + Io sin {| — ——}}, 
Mal 5 5 


ZOOS Re ING 
or bn SS any oa 


men? oy 
agreeing exactly with the harmonic analysis given on 
page 91, which it obviously should do. 

In a second manner of applying this method we 
draw a parabolic arc through the tops of the three 
ordinates ., , ¥23 another one through the tops of 
ey V3, Vy and so on; we thus get a function with dis- 
continuities at the even ordinates y,, yo, yy. . -, and 
at these ordinates there will, in general, be a discon- 
tinuity both of slope and of curvature. We must 
find an expression for the magnitude of these dis- 
continuities. 


FOURIER’S THEOREM 105 


If we take the origin at y, it will easily be found that 
the equation of the parabola through the tops of y,, y, 
and ‘yp. is 


Naas +44 29 Roe meee ata a, t?, 


2a2 


where a is, as usual, the distance between the ordinates. 
Similarly, ‘the equation of the parabola joining ys, Je 
and y, referred to an origin at y, is 


Pern vik) Neos bak er 
y Ys 3 2a i i 2a" 8 


Differentiating and putting t = — a, we see that this 
parabola has a slope of 


= Maer 49s = "342 
2a 


at Vo. 
Similarly, it will be found that the first parabola 
has a slope of 
Nest 4) Vo 
2a 


at the same ordinate; and so the discontinuity of 
slope at Yq is 


fS Sento Oy, + 4V1 — Yo 
2a ‘ 


4 
or — es where A4y, is the fourth difference of the 
a 


y’s which occurs vertically below yp. 


, > ee eed? 
We must now find the discontinuity in For 


dt? 
the second parabola we see that - as he ou ye. 


. ‘ ‘ A’ De d* 
or in our difference notation, eos Similarly, oi 
a 


106 A PRACTICAL TREATISE ON 

for the first parabola is Ae, and hence the discon- 
a 

Sill. Cee , p ATVs =e : 

tinuity in os at the ordinate ‘yz 18 ere that is, 

the difference between the second differences which 


lie each side of that one below y,. Replacing a by 
x/p and substituting in the formule (15) and (16) of 


page 59, we get 


an = ee BaCOS ee 
2 F . 2nT 
=i Fd AY, — A?y,) sin ay. a eee 
and 
n= fly sin = 
2 2nm 
= FB A*ys — A?y,) cos ana eee 


The summation extends to all the ordinates at 
which there are discontinuities. "The expressions actually 
written under the %’s represent the discontinuities 
at y, only. 

If a’, and b’» are the ordinary harmonic analysis 
coefficients of the function with ordinates A*y, at yo, 
etc., and a”, and b”, are those of the function with 


ordinates A*y, — A®y, * at Ye, etc., we may write these 
coeflicients 


ae eS n a 7373 n . . . (25) 
and 


eee ae ene ea ne UE Ses 


With the obvious notation of using the suffix 1°5 to indicate a 
term standing mid-way between terms of suffixes 1 and 2 we have 


A?ys — A*y, = A®ys.s + ASyy.5 while Aty, = Ay,., — ASy1.5. 
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As an example of this, we will take the case of an 
even Class IV function whose ordinates are 2-000, 
mu 
> 6” 3 
to the symmetry, we need only obtain the differences 
over a quarter of a period, but to do so we must write 
down two ordinates before and after the quarter 


period chosen. We thus get the Table 


14 . . 
1-732, I-000 ando ato and respectively. Owing 
2 


bd aa ee ° es i = es oe 
3 6 6 3 2 3 6 

yi 12000" “3'732) | 2"000' 1°732 11000 .@  —1000 =1'732 
Ay "732 ©=-268 268 — +732 —1-:000 —1:000 —+732 
A*y —-464 —+536 —-464 -—-268 0 +268 
A’y —:072 :072 196 +268 “268 
Aty *T44 "124 «+072 
A*y,—A’y, o ‘268 © +464 +536 


Thus, in a quarter of a period, we have two discon- 
° Ged . . MA . 
tinuities in A*y, viz.,atoand~—. Ina complete period, 
there are two discontinuities similar to the first, 


and four similar to the second. In A?y, — A?y, we 
only have one discontinuity in the quarter period, 


viz., at ~. Substituting in (23) and remembering 

that 2 is odd and p = 6, we get 

Oe 6 te x3 cos o + x 072 cos 
"ont? 44 4 3 


as el4 x +464 sin =| . 
Working out this formula, we find 
f(2) = 19975 cos t — -0097 cos 5¢ + :0081 cos 7t 
— +00033 cos IIt + -00163 cos 13¢ 
+ -00007 cos 17t + :00064 cos I9t+ ... 
which agrees remarkably closely with the simple curve 


f(t) = 2 cos ¢ which the given points evidently lie on. 
As an improvement in this method we may also 
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draw parabolic arcs joining the ordinates ,, Ye, Vs and 
Yay Var No» etc. and then add the analysis of this case 
to that of the former case and divide by two. This 
will be allowed for if we sum the expressions in (23) 
and (24) for every ordinate, instead of only for the even 
ordinates, and divide the result by two. We leave 
it as an exercise for the student to apply this method 
to the above example and find that the agreement with 
2 cos ¢ is made considerably better. 

A third manner of applying this method of harmonic 
analysis is to draw a curve joining the tops of ordinates 
y, and y, such that: at y, the curve is parallel to the 
straight line joining the tops of y, and 42, and at y, 
it is parallel to the straight line joining the tops of y, 
and 4s. Since this is virtually drawing a curve through 
four points, the curve must be a parabola of the third 
degree. ‘The tops of all the other consecutive ordinates 
are joined in a like manner. This function then has 
no discontinuities either in magnitude or slope; but 
it has discontinuities in its second and third differential 
coefficients. 

We leave it to the student to prove that the equation 
of the curve joining y, and ys is 


t 2 
Yigact Oni Yi) rate ate eee 592 + 2) 
3 


+ (%4 — 398 + 35a — My) 
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the origin being taken on y,. The equation of the arc 
joining y, and y, can be written down from this by 
writing y, for 3, y, for y, and changing ¢ into — ¢; 
and so is 


t 2 
same aeiG perth snl betintnar Whee She 2s) 


WES RENTED  S 


referred to the same origin. 
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The third differential coefficient of the first curve is 
constant and = 3(y, — 373 + 372 — 4,)/a3, while for 
the second it is 3(y, — 3. + 3, — y,)/a8. Hence 
the discontinuity in this coefficient at yp is 3(y4 — 4/3 + 
62 — 4y, + y,)/a3; that is, 344y,/a? in our previous 
difference notation. 

The second differential coefficient of the first curve 
at t= 0 is (— ¥,+ 43 — 52 + 2y,)/a2, while for 
the second curve at the same point it is (2ys — 57, + 
4, — y.)/o?. Hence the discontinuity in this co- 
efficient at yp is (— y, + 273 — 2y, + y,)/a?; that is, 
it is (— A®y, + A*y,)/a?. If we reduce these suffixes 
by unity so as to get the discontinuities at the ordinate 
y, and then substitute in the formule (15) and (16) 
of page 59 we get 


: Ast 
ay, = — F249, — A.) sin oe 


+ SEZ At, cos Peon 27) 


and 


i 42? ae 
3 
= 3P EE Aty, sin zs hope a Pep 


Or if a’, and b’» are the ordinary coefficients of 
harmonic analysis for the function with ordinates 
Aty, at y,, etc., and a”, and b”, similarly for the 
function with ordinates A?y, — A®y, at y,, etc., we 
have 


3 4 
an = — Fa b'nt Fan. . . (29) 
and 


3 4 
PS we. ain + Ff ioe eg A ey 


men 
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It is noteworthy that we are concerned here with 
exactly the same differences as we were when dealing 
with ordinary parabolic arcs joining three ordinates. 

Applying this case to the same example as we took 
previously (page 107), we get 

Be eel by teeta gin re: 
an = — 3 14 X °200 sin 6 + 4 X °404 sin 3 


n3n3 | 
+ 2X +536 sin =) 


E104 ON eae det 
Poe 144 cos O + 4 X *124 COs 6 


+ 4 X%**072,c0s =I 


Working this out, which is not the trouble it appears 
to be, for after having worked out a, there is very little 
to do save divide the same figures by n° and 4, we get 

f(2) = 1°99814 cos t + -00002 cos 5¢ + 00138 cos 7t 

-+ ‘OOOOOI cos IIf + *00003I cos IIt 

+ +000000 cos 17#+ ... 2 foamy 
a remarkably good approximation to 2 cos ¢t. It may 
be objected that the ordinary method of harmonic 
analysis would have given an even better approximation 
as it would have given 2cos ¢ exactly. Quite so, but 
that is accidental and due to the fact that the selected 
points lie exactly on the curve 2 cos ¢ which is one of 
the terms of the assumed harmonic expansion. Had 
we, on the other hand, used the very similar function 
8z(a — t)/x® and taken 12 ordinates per period and 
joined successive even ordinates by a parabolic arc 
passing through the intermediate odd ones, our 
present method would, in this case, have given exact 
agreement with the Fourier analysis of the function, 
which is 

f(t) = 2:0641 sin ¢ + -07627 sin 3 + -01651 sin St 

+ -oo602 sin 7¢ + 00283 sin gi+ . : 


<5 S 
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whereas the ordinary method of harmonic analysis 
with the same number of ordinates would have given 
us 


f(t) = 2-0577 sin t + -07407 sin 3¢ + -01630 sin 52, 


the higher harmonics being quite undeterminable. 
The accuracy of this is “miles behind” the accuracy 
with which (31) represents 2 cos 7. 

In the last manner of applying the present method 
the labour is about twice that of ordinary harmonic 
analysis apart from the labour of finding the fourth 
differences, since we have two different discontinuities 
at each ordinate and each discontinuity has to be 
treated just as if it were an ordinate. The present 
method is not put forward for use where any rough 
analysis is good enough, but only when the highest 
accuracy is required, and especially in those cases 
where we require the amplitudes of one or more of the 
lower harmonics very accurately, together with as 
high an accuracy as is obtainable, on the higher 
harmonics. ‘The last variety of the method described 
should then always be employed. 

The student who ever has any practical harmonic 
analysis to do will soon find he has had ample exercise 
in the subject. We would recommend all students, 
however, to construct one or two Schedules for a 
different number of ordinates per half-wave-length 
than those we have given. 


Interpolation Formule for Periodic Functions. 


When the value of a periodic function is known for 
non-equidistant ordinates the harmonic analysis can 
be performed by first writing down an interpolation 
formula and then transforming it by the method 
described in Chapter IV. 

If we have points, (a,A), (,B) .. . (&,K) (,L) 
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the equation of the parabola of the m( — 1)th degree 
passing through the points is evidently 
a (8 De) ea) 
I) (ara By ey ee 
(4 — ax— 6) 22 5 = B®) 
12 (I a) ey ee 
If now we imagine all these points lying between 
o and 2 and imagine the whole set repeated an infinite 
number of times at intervals of 22 to both — o and 
++ 00, we find that in the limit (32) becomes 


_sn#{a—b)...snk«x—J 
NGI sin #(a — db)... sin #a—J) aa 


sin H« — a)... sin #(x — k) 
in HT tal ee ee ae 
when 7 is odd and 


f(x)=cos 3(”—a) 


sin #(a—b) .. . sin #(x—l) 

ana oh 2 Sade) 
xv — p) SH He — 4)... sin (x — b) 

108 Ke 7) ae oe ea 


(34) 


when x is even. 

The proof of these formule, though rather tedious 
to write out, is quite easy if the following results be 
borne in mind: 


= eT I ao x 
sin —*¥ = —v{ I — —— I — —~— Aye 
2 2 4n? 1672 
I 


$ cot w=! ie : 
z Os Xv -+ 27 x — 22 X-b 4m 

I vy I 
a x+ 27 x — 272 x + 4m 
The expressions in (33) and (34) obviously satisfy 
the conditions of passing through the given points 


a ae 


I 
-COSCG == 
2 


Ht Hee 
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and, since each of the factors in the numerators 
merely changes sign when » is increased by 2x and 
there is in both cases an even number of such factors, 
the period is evidently 27 in each case. 

We know that if in (32) each given ordinate be 
increased by the same amount, f(x) will be increased 
by this amount for all values of x; that is, the formula 
(32) still represents the same curve : we conclude there- 
fore that formule (33) and (34) likewise satisfy this 
condition. 


Ux 
If in (33) and (34) we write ¢2 =,/u and arrange 
the result in powers of u it will be seen that no 
fractional powers occur and that further, if the result 
is expressed in sines and cosines of multiples of x, 


the highest harmonic when n is odd is the i es lth 


- ° n 4 
and when 1 is even is the 5th. But these are just the 


harmonics which we found it necessary to assume present 
for our harmonic analysis when taking equidistant 
ordinates: we infer, therefore, that for equidistant 
ordinates (33) and (34) give the same result as is 
attained by ordinary harmonic analysis. 

When the ordinates are equidistant the expressions 
(33) and (34) may be greatly simplified. Let the 
points be (a,Y;), (2a,Y2) . . . (2,Yn), where a = 2n/n. 
Then it may be shown, by a well-known factorisation 
result in trigonometry that (33) becomes 


HO? ee “sin saa — Y, cosec #(” — a) + 


Ye cosec Hx — 2a) Sh =e We, cosec (x = 2m) (35) 


which holds only when is odd; while similarly (34) 
becomes 
I 
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Kx) = = sin na Y, cot #4” —a)+ 


Va COt 4 (4 — 20) oe + Yn cot (a — 27)| aaa) 


which holds only when a is even. ‘That these expres- 
sions satisfy the conditions of passing through the 
given points and of having a period of 2m can be very 
easily verified. These formule, then, are mathematic- 
ally the same as the results arrived at by ordinary 
harmonic analysis: they suffer from the results not 
being in so useful a form from a physical point of 
view, but they possess the advantage that they can 
be immediately written down without any calculation. 

Formule (35) may be usefully employed for inter- 
polation with non-periodic functions when these are 
given for a fairly large number of equidistant ordinates. 
They are much more easily written down and calcu- 
lated from than (32) and give practically the same 
except near or beyond the end-points. 

We will now simplify (34) in the cases of functions 
belonging to one of our four Special Classes. 

First let f(~) be an odd periodic function and let the 
m given points (a,A) . . . (/,L) be considered as lying 
between o and x. If we make (34) in addition to 
passing through these points also pass through the 
nm points (— a, — A), etc., it will be found to become 


_sin x (cosx—cosb) .. . (cos x—cos /) 
Kx) sin a (cos a—cosb) . . . (cos a—cos/) bee? 


sin x (cos ¥ — cosa) . . . (cos x — cos k) 


sini ‘(cos 1 — cosa) . . . (cos 1 — cos &) 


(36) 
Similarly, if f(«) is an even function and if we make 
(34) in addition to passing through the » given points 


also pass through the points (— a,A), etc., it will 
become 
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fi) = (cos * — cosb) . . . (cos x — cos J) 
(cos a — cosh) . . . (cos a — cos /) 


A+... 


(37) 
Again, if f(~) is an odd harmonic function, and if 
we make it pass in addition through the » points 
(a — a, — A), etc., it will become 
po sin (4 —")) sin (4 — )) 
eae ay an Ge 
sin(#— a)... sin(« — k) 
ee eae 
if 2 is even and 


Fa) = c0s (w — a) FO At... 
Bee 5 (39) 
if 2 is odd. 


Finally, if the f(*) is a Class IV function and if the 
nm given points be considered as lying between 0 and 2/2 
and if we make (34) pass through similar points in the 
other three quadrants we get 


fle) = sin x (cos 2x — cos 2b) . . . (cos 2¥ — cos 2/) 
sin a (cos 24 — cos 2b) . . . (cos 24 — cos 21) 
ee ae AL LAD) 
for an odd Class IV function and 
cos % (cos 2% — cos 2b) .. . (cos 2% — cos 2)) 
ie) = cos @ (cos 24 — cos 2b) . . . (cos 24 — cos 2) 
See ee LenS 


for an even Class IV function. 

These results are most readily obtained from (36) and 
(37) respectively. All the formule (36) to (41) 
evidently satisfy all the conditions imposed upon 


them. 
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The formule (36) to (39) can be further simplified 


to half the number of factors in each term if the 
ordinates are chosen symmetrically so that =a — a; 
k=—2a— 5, etc. In this case, we will suppose the 
ordinates at a and a — a to be A, and A, respectively 
and similarly for all the other ordinates. 


We then find that for an odd function (36) becomes 
ai= ‘fe + A, sin x it A, — A, sin 22 e 


2 sin @ 2 sin 24 
(cos 2% —‘cos 2b)).. = (cos 2x cos2)) 4 (42) 
(cos 24a — cos 2b) . . . (cos2a—cos2)) 


Likewise for an even function (37) becomes 


Aen = {os — A, cos x & A, + A, cos 221 ¥ 


2 cos 4 2 cos 24 
(cos 24% — cos 2b) . . . (cos 2% — cos 2/) rs (43) 
(cosi24a — cos 2b). . (Cos 2a/—€as 2)) 43 


While for an odd harmonic function both (38) and 
(39) give 
eg {a + A, sin x an A, — A, cos “| $s 


2 sin a 2 cos 4 
(cos 24% — cos 2b) . . . (cos 2% — cos 2l) (44) 
(cos 24: cos2b) .. . (cos 2a— cosa) iD 


In these three formule we have supposed that all 
the ordinates signified by the letters from a to 1 occur 
between o and 2/2. The term written down is the 
one for the ordinate at a and a similar term must be 
added for each of the other ordinates from 8 to 1. 

In any of the formule (33) to (44) we can take two 
of the ordinates undefinitely near together and so 
obtain the equation of a curve which has a given slope 
at a given point. By taking three ordinates indefinitely 
near we can obtain a curve with both a given slope and 
a given curvature at a given point. ‘Two examples 
will suffice to illustrate this point. 
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Example 1. Find the odd Class IV function of 
period 22 passing through the points (6,$6) and 


(4 H), when 6 is infinitesimal. 


Formula (40) gives us at once 
sin * cos 2x — coszx 
CO 
sin 0 cos 20 cosz 
sin * cos 2% — cos 26 
sin 7/2 cosa — cos 26 


or 


J(&) =S-sin x cos? + H sin3x 


which may be written 
S iH 
fle) = 2 cae 


Example 2. Find an even Sate curve of period 


sin 4 -- Dies Sin ax. 


2x through the points (0, 0) (Qe —— ) and (@, H), when 


6 is infinitesimal. 


Here (37) gives us at once 
i) = (cos « — cos o)(cos x — cos 2) Réd? 
(cos 6 — cos o)(cos 6 — cosa) 2 
(cos « — cos 0)(cos « — cos 9) 
(cos % — cos 0)(cos % — cos 0) 


or ieee sin’ a+ At — cos x)?. 


Problems of this kind, however, are often best solved 
by writing down a Fourier series, with the least number 
of the correct kind of harmonics which are seen to be 
necessary, and then determining the coefficients to 
satisfy the given conditions. 
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EXAMPLES. 
1. In an odd Class IV function, if 0, E,, E,, E, are the ordinates at 
o, = = = respectively, prove that 
Wi pin TENE: : 
by = 3(Ey + E,V3 + E,), 


bg = 4(2E, — E,), 


bs rz 3(E, aa E,V3+ F3). 


2. Show that b,, bs and b, in Ex. 1 are respectively, 
b, = 5, — by + 43 — O93 + F955 — - 
bs = 5 — 99 + Py, + Dig + Ogg — - 


3. Analyse the odd Class IV function defined by f(#) = 2? from 
™ 


™ A : : 
E=0 to— by using 3, 4, 5, 7 ordinates per quarter-period ; and show 
that in the first case the coefficients obtained satisfy the relations in 
Ex. 2 and satisfy similar relations in the other cases. 


CHARTER VI 
THE THEORY OF FOURIER’S SERIES 


Tue present chapter may be entirely neglected by 
those whose interest in the subject is purely practical. 
Nevertheless, we hope we have made the subject so 
far sufficiently interesting to induce many, who have 
no intention whatever of making mathematicians of 
themselves, to peruse this chapter out of an interested 
curiosity in the subject. 

For the possible benefit of a few readers we will 
first recall one or two theorems on the convergence of 
infinite series. First, any infinite series s, — 5, + 53 — 
5g +... in which the terms are alternately positive 
and negative, and in which each term is numerically 
less than the one before it and such that the terms 
decrease ultimately to zero, is convergent. For, 
writing the series in the form (5s, — 5) + (53 — 5,)+... 
we see that the sum is positive; and writing it in the 
form s, — (5, — 53) — (54 — 5g) — . . . we see that it 
is less than s,. Hence the series is convergent. 

Secondly, the series 1+4+4+4+... 1s 
divergent. For if the reader constructs a series of 
adjacent rectangles of unit width on the same base and 
of heights 1, 4, 4, etc., and also draws the curve 


Ve | through the vertices so as to lie inside the 
a 


rectangles he will see at once that the sum of all the 


oc 
rectangles is greater than \ a which equals log « 
I 


which is infinite. 
11g 
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Thirdly, the series 1 + = + A + Ssh at 


convergent when n> 1. For if we draw rectangles 
as in the last case, and a curve through the vertices 
lying outside them, we see that the sum of all the 
wae 
aes 


terms after the first is less than \ which equals 
I 
I 
n—1 
We will next prove that if 5, 5, 53 . . . be any 
series of regularly decreasing positive terms decreasing 
to zero ultimately, no matter how slowly, then the 


series 


5, cos B + 5, cos2B + Ss cos3B+ . . (1) 
and Soin fy Si 26 45 ee ee 


are convergent whenever f can be expressed in the 
form px/q, where p and gq are integers of which p must 
be odd. Whatever the value of f this can be done to 
any desired degree of accuracy. If, ¢.g., B = 47/7, the 
theorem will hold for B = 4,000,0017/7,000,002, so 
the two series can only be divergent for certain isolated 
values of f. 

Consider the first 2g terms of either of these series 
and divide these terms into two groups, the positive 
and the negative terms. Let the sum of the first be 
denoted by P, and of the second by —N,. If 
Sy cos rB, say, is any term in P, then 5,44 cos (r + 4)B, 
that is, — sy4 cos 7B is a term in N, which by hypo- 
thesis is numerically less than the corresponding term 
sy cos 7B in P,. 

Similarly, all the terms in N, are less than the 
corresponding terms in P, and hence N, <P. If 
now we consider the terms from 52941 to 54g and denote 
the sum of the positive and negative terms by P, and 
— N, respectively, we see, by a similar argument, that 
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P, <N, and that N, <P, and so on. Hence the 
series of terms, 

Pe eget Ng Po), 
are such that each is numerically less than the preceding 
and of alternate sign. It follows in this manner that 
both the series (1) and (2) are convergent. 

Writing 6 + ¢ for B in (1), we see that 


5, cos B cos t + sg cos 28 cos 2¢ + “| (3) 
— 5,sinfsint — s,sin2Bsin2t— . je 


oe 


is convergent as a whole, for all values of ¢. But the 
first line is an even function of ¢ and the second one an 
odd function. Now the infinities of an even function 
cannot everywhere neutralise the infinities of an odd 
function, for if they neutralised for positive values of ¢ 
they would re-inforce for negative values. It follows, 
then, that the series for both the odd and the even 
functions must be separately convergent; and also, 
that a simular result must hold for the two series, 

| 


5, cos B sin ¢t + 5, cos 28 sin 2t + | 4) 
and Sg sin B cost + 5,sin 2B cos 2#+ . Mess 
obtained by putting 6 = f + #1in (2). 

One obvious exceptional case must now be noted, 
that is, when 8 = oora multiple of 27 in (1). All the 
terms are then positive and the series 5, + 52+ 53+. . 
may be either convergent or divergent, according to 
the rapidity with which the s’s decrease. We conclude 
from this that either of the series in (3) may be infinite 
when f + ¢ is zero or a multiple of 22, but must be 
finite for all other values of 1.* 

Each term of the series (2) is identically zero when 
B =o, but the series may then behave in a very 


* If the reader draws, from an origin, radii of lengths 5,, 5... ., 
in the directions B, 2B . . ., respectively, it will be nearly (if not 
quite) obvious to him that the sums of (1) and (2) are necessarily finite 
for all values of ® save zero or a multiple of 27. 
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peculiar manner which must now be investigated. 
Let us suppose that £ is indefinitely small so that we 
can replace sin mB by nf, the series then becomes 


5,8 + 25,8 + 35,8 +... 
or Bly A 25g easy Free <9: 

Now the series in brackets may be convergent or 
divergent: if convergent there is no doubt about the 
sum of (2) being zero when £ = 0, but if divergent 
the result takes the form o X o when Bf = Oo s0 that 
in this case we cannot say anything about the value 
of the series 

fSin Pct sisin 2p eae 


when f =o except that it is formally zero on the 
strength of the fact that each term is zero. 


I et 
Let us suppose that sn = —, so that the series is 
n 


spat pee mae ES 


We may write this 
Sift Pay SSID 2p ee sin ap sin NB | 
Ble tape ap tt meet |: 


The student must now recollect that the limit of 
the sum 


BL/(B) + f(28) +... + f(NA)] 


when f = 0 is 
Np 
" S(x)dx. 


by the definition of a definite integral. 
Hence the limit of the sum of the above series up to 


the Nth term is 
N6.: 
\ sin #5 
a eS 
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If now we make B = 0 and N = «, we do not get 
a unique value for the result, but a value which is a 
function of the product Nf. We conclude that the 
given series, when f= 0, may have any value that 
the integral 


can take. ‘This integral clearly increases with X from 
X = 0 to X =a, when it assumes the value 1-8519; 
after then it decreases till X = 2%, owing to the 
additional elements being negative, when it takes the 
value 1°4181; it then goes on oscillating with a 
diminishing amplitude till X = o, when it has, as we 


saw on page II, the value of < or 1°5708. 


This is the behaviour when # has an infinitely 
small positive value. Since the series in question is an 
odd function of #, it behaves in an inverted manner 
for negative values. 

Iience the series (5), swhen 6 = 0, ican” assume 
any value between — 1°8519 and + 1:8519. We 


shall call the value 2 the principal value of the series 


when f is + 0, and — the principal value when 


fB = — 0, and zero the formal value of the series when 
B=o. The principal value corresponds to the value 
when Nf = «, and is thus the value of the series 
when f is an exceedingly small finite quantity, or 
7 
ze 

We can similarly investigate the behaviour, when 
B = 0, of the series 


gs, sin B+ 5,sin2B+... 


as we may express it, (+ 0) = . and ¢(— 0) = — 
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when fn, = 3 and p is any positive quantity. The 
n 
series becomes in this case 
sin B , sin 2B sin NB ere | 
pr pP a 2? pP rm Os ay N?p? 


Fic. 20. 


which, when £ tends to zero, becomes 


NB .: 
aaa sin #7 etted Nis (6) 


xP 


The student familiar with definite integrals will 
recognise that this integral is finite and convergent 
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whenever i021 OOw if I <p <2 the result of (6) 
is zero on account of the factor B?—1. It can be 
proved that the result of (6) is zero even if Die 2, 
since the zero of the factor f?—' is of higher order 
than the infinity of the integral owing to the factor 


Fic, 21. 


xP in the denominator. If, on the other hand, lies 
between o and I the result of (6) is infinite when f = o, 
unless the range of integration Nf is infinitely small, 
owing to the B?—* factor which is then infinite. In 
this case the series can take any value from 0 to « 
when f = + 0, and any value from o to — « when 
B = — 0, since it is an odd function. 
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The appended diagram (Fig. 20) shows qualitatively 


the behaviour of the series 
sin sin 2 
p , sin ap 
1? 2P 


for the values of » of 4, 1, 3/2, 2 and 3, in the neigh- 
bourhood of B = o. 


Fig. 21, similarly, illustrates qualitatively the 
behaviour of the series 


++ 


COS 6 9 casi2p 
1? oy: 2P 


for the same values of p in the neighbourhood of 
PieOe 

In both these diagrams the value of the series is 
infinite when p <1 and also when p = 1 in Fig. 21. 
The length of the heavy vertical line in Fig. 20 shows 
the extreme range of undeterminateness of the series 
(5) when 6 = 0. Both arcs for the: value p 13) as 
in both these figures, have an infinite slope but a 
finite radius of curvature at 8 = 0; in both diagrams 
the slope at 6 = 0 is Minite if @-2eand initemias 
p = 2 in Figeer; the curvature at f — 0 1s Hniterin 
both cases when p> 3 and for p = 3 in Fig. 20 if 
the arcs on the two sides of the origin are considered 
as separate curves.* 


a 


The Relation between the Fourier Coefficients and the 
Properties of a Function. 


It is not necessary that a periodic function f(/) 
should be everywhere finite for it to be expanded in 


* The arcs on the two sides of the origin are necessarily represented 
by different analytic functions, and each arc has a finite radius of 
curvature at the origin. ‘This must not be confused with the point 
of inflection of an analytic curve for which the radius of curvature 
increases to infinity as the point of inflection is approached. 


FOURIER’S THEOREM 17, 


a Fourier Series; but it is necessary for all the Fourier 
constants to be finite, and for the higher ones to 
decrease indefinitely, if the series is to be convergent. 
The first condition will obviously be satisfied if 


ou bere aC. eet ( 7) 


an 
is finite; for the integral | f(z) cos nt dt consists of 
ie) 


elements none of which are larger than corresponding 
elements in the first integral; while some are taken 
positively and some negatively in the second case, and 
all positively in the first case. 

We will next prove that if f(#) is everywhere finite 
the Fourier constants a, and by cannot be greater, 
when 1s large, than c/n, where c is some finite constant. 
In other words, the Fourier constants, 4,, a2, dg . . +5 
must decrease, ultimately, at least as fast as the series 
Piet O) 230/325 +. 

The function f(#), in the interval from 0 to 2m, can 
be divided into a finite number of ranges in each of 
which range it is of one sign and either constantly 
increasing or constantly decreasing (or stationary). 


Let us consider the value of the integral \ f(#) cos nt dt 
over any of these ranges when 2 is fairly large. Over 
a short range of ~ the integrand will be positive and 

n 


will contribute s,, say, to the integral. Over the next 
range of z/m the integrand will be negative and give 
— 5g, say, and so on up to 5, where s, corresponds to 


the last complete half-wave in the range considered. 
Now the terms 


Sa Sg217 Sgt — 28 @ ae hy 


are alternately positive and negative and either steadily 
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increasing or steadily decreasing. In either case we 

can see by a similar argument to that used on page 119 

that the sum cannot exceed either 5, or 5, whichever 

is the greater. Also, it is clear, neither of these terms 
2 2M 

can exceed “.M.“ = ~~, where M is the greatest 
7 n n 

value of f(t) in the interval considered. To this result 

there may have to be added the contribution from any 

fractional half-waves at the beginning and end of the 


: : 2 
interval, but neither of these can exceed ——, so that 
n 


the value of \ 40 cos nt dt over the range in question 


cannot possibly exceed —-. Summing for other ranges, 
n 


we see that the whole integral from 0 to 2m cannot 
exceed c/n where c is some finite constant. The same 
result obviously holds for the sine integral. 

From .this important result we can easily prove 
that if f(z) has no discontinuities in magnitude, and if 
f'(2) is everywhere finite,* then the Fourier constants 
an and by cannot exceed c/n? where c is some constant. 


For 
' el ee 
an + thh = — \ f(tem dt; 
10) 


It 


integrating this by parts, we get 
an + tby = [AF] 


If f(z) has no discontinuities in magnitude the first 
term vanishes, since it is periodic. ‘The integral gives 
the Fourier constants for the function f’(#) and these, 


2m 2 
ea \ f'(Deim dt. 


Gu Stic 


* The first condition has been stated to prevent misconception : 
the second condition includes the first. 
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we have seen, are <c’/n where c’ is some constant, 
since f’(#) is everywhere finite. Hence an + iby, and 
therefore ay and by, are < c/n? where ¢ is a constant. 

In a similar manner, we can prove that if there are 
no discontinuities in f(#) and /’(z), and if f’(z) be 
everywhere finite, then the Fourier constants a, and 
bn must be <c/n® where c is some constant, and so 
on for discontinuities in the higher differential 
coefficients. 

We will now investigate the order of magnitude of 
the higher Fourier constants when f(t) becomes 
infinite at one or more points. Ifit becomes « at more 
than one point in the range 27 we can resolve it into a 
number of functions each of which becomes infinite 
at only one point. Further, in considering any of 
these functions, we may take the origin at its point of 
infinity. We will accordingly suppose that any one 
of these functions can be represented by 


fli) = 4 + $i) or by f(t) = kloge + $()), 


where ¢(#) is under no restrictions save that it 1s every- 
where finite; so that its Fourier constants, when n is 
large, cannot exceed c/n. ‘The Fourier constants for 
the first term are given by the real and imaginary 
parts of 


27 
k eint 
ao, tf 


Let us write nt = x so that the limits are now o and 
2an; the expression then becomes 
2TH 
k e* dx 


spe aes gee 


If p lies between o and 1 the integral is finite when 
m is made oc in the upper limit, so we see that the 
K 
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Fourier constants are of the order of magnitude not 
greater than c/nt~? when 1 is large. Compared with 
this, the constants due to the finite function ¢(#) are 
ultimately negligible. 

Hence we see that the order of magnitude of the 
higher Fourier constants depends on the nature of 
the infinities of f(2) if it has any, but if not, it depends 
upon the nature of its discontinuities in magnitude ; 
if it has neither, they depend upon the infinities or 
discontinuities of f’(¢) and so on. 

It should be noted that if p<1 the integral 


27 
\ | f(2)|dz is finite; for then 
e) 


2m a 
\ pence, 
oO 


We thus see that whenever the condition that 
Bine 
\ | f(2)|dz is finite is satisfied, the Fourier constants 
fe) 


decrease indefinitely with 7, since both these conditions 
are satisfied at the same time, viz., when 7 lies between. 
Oand 1. 


The Convergency of a Fourier Series. 
The convergency of a Fourier series, representing 
any function f(z), say the series 
a, cos t + a, cos 2¢# + <5 
only needs to be considered when the series 
[4 | + l@e| + |as| + 
is divergent, for if this series is convergent the first 
series is obviously convergent. 
It follows from this, and from what we have said 


above, that it is only the effect of the infinities and 
discontinuities of magnitude of f(#) on the Fourier 
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constants that need to be considered here; for 
discontinuities of slope or curvature only contribute 
a term of the order ¢/n? or ¢/n® to the Fourier constants 
which could not affect the convergency in any way. 
Let the function /(#) be represented as the sum of 
several functions each of which has only one point 
of infinity or of discontinuity of magnitude. When 
any of these functions are being considered let this 
point be chosen as the origin. 

Now any periodic function of period 22, which 
has no infinities, and only one discontinuity, and 
that at the origin, can only differ from a multiple 


of the odd function defined by (2) <7! from 


= 0 to t= 2m, by a function whose Fourier series 
is absolutely convergent and therefore need not be 
considered. But the former function is represented 
by the series 

sin 7 Ae sin 2 , sin 3¢ at 
I 2 3 

(see page 42), and this serzes, we have seen above 
(page 120), is convergent for all values except for t= 0 
(or a multiple of 22) when it is indeterminate and can 
assume all values from — 1-8519 to + 1°8519. Since 
the function at t = o is indeterminate within the range 


from —™” to +, the range of indeterminateness of 
2 2 


the series is 1:179 times the magnitude of the dis- 
continuity of the function. ‘The formal sum of the 
series is zero, which is the point midway between the 
extreme values of the function at its discontinuity, 
and this property will obviously be preserved if this 
series is added to another series which represents any 
other function with no discontinuity or infinity at t=. 

Similarly, if the function has a single infinity at 
t = o, and in the neighbourhood of t = 0 approximates 
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to k/t?, where p lies between o and 1, we have seen 
that the series representing it approximates, in its 
higher terms, to 


COSEE COSeZE cos 32 
vl ae ey 
VES 


| 
21-P g1-p pret 


or to 


Sy t =P a 31-P OY ene i 

where ¢ and c’ are constants, according to whether 
the function is odd or even: if it is neither, it can be 
represented by the sum of both these series. It is, 
of course, only the behaviour of the higher terms of a 
series which affects its convergency. Now both these 
series, we have seen, are convergent for all values of 
t save t = O (page 120). 

When ¢ = 0, the first series, we saw, was infinite ; 
while the second was indeterminate over an infinite 
range, which is exactly how the graph representing 
the function behaves in the two cases. 

We have thus seen that for any function, with any 
finite number of discontinuities, or infinities which 
are such that the integral 


J ixalas 


is finite, the Fourier series zs convergent and can thus 
be used to calculate the function. 


We may note here that it will easily be found that 
the series 


Ae Pe sin 2h sil eae 


COSaz cos 2t COSEAE 
ef pees he oe 


: ; pt es 


which we have not mentioned represents the artificial 
periodic function defined by f(#) = log |fcosecdy|. 
When fis small this becomes — log | ¢| approximately ; 
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which is infinite when ¢ = o, but is a lower order of 
infinity than is represented by 1/2?. 

Since a periodic function (whether even or odd) 
which has a single infinity of the order k/t? at t = 0, 
where p lies between o and 1, has its higher Fourier 
constants approximating to c/n'~? we see, by integrat- 
ing both the function and the series (s — 1) times, 
that a series in which the higher Fourier constants 
approximate to c/n‘? represents a function whose 
(s — 1)th differential coefficient has an infinity of the 
order k/t#? at t= 0. Also from the foregoing results, 
we see that if the higher Fourier constants approximate 
to c/n’ where s is an integer, then the (s — 1)th 
differential coefficient will have an ordinary discon- 
tinulty in magnitude if it is an odd function and will 
have an infinity of the nature of log|#| if it is an 
even function. 

We may remark that the functions in the last 
paragraph are not necessarily artificial functions 
though they usually are. For instance the function 
defined by y = sin(#+ ¥) is an analytic function. 
It represents an odd periodic function of period 2z ; 
a distorted sine wave which has an infinite slope at 
the origin and attains its maximum value of unity 
when t= 2/2 —1. Near the origin we have t = y3/6 
approximately, or y= 61/3z1/3 s from which i = kt~?!3, 
so that the higher Fourier constants of the differential 
coefficient approximate to c/n"/3; and so the higher 
Fourier constants of the function itself approximate to 
c [n'l3, 

The Fourier constants may, of course, decrease 
much more rapidly than is represented by the law 
c/n’; an and b, may, for instance, when n is large 
approximate to ck” where k <1, as in Problem I 
of Chapter IV. In this case, if we differentiate the 
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series s times, a» and by for the differentiated series 
approximate to ck”n* ; and the ratio of one coefficient 


3 : nN Tn\s : ° £6 
to the previous one 1s e( x ) , which, in the limit, 


n 

when 1 is infinite, is k, so that the series is still con- 
vergent. Hence such a series represents a function 
in which none of the differential coefficients are ever 
infinite or have any discontinuities. Such a function 
is necessarily an analytic function.* 


Proof of Fourier’s Theorem by Direct Summation of the 
Series. 
We have proved in Chapter II that if any periodic 


function of period 2x can be expanded in the series 
40 4. a, cos t+ iia) A eee 
2 


then the coefficients must have the values there found 
for them. We also gave later a proof that the expan- 
sion was possible (page 23). Just how far the reader 
will regard these two proofs taken together as constitu- 
ting a rigorous proof of Fourier’s Theorem will depend 
upon what we may call his mathematical temperament. 

Another proof of Fourier’s ‘Theorem has virtually 
been obtained in Chapter V. We there obtained a 


* ‘The student familiar with the elements of the Theory of Functions 
will recognise that the Fourier series in this case will be convergent 
for complex values of ¢ as well as real ones. If we write t= u-+ iv 
ny Ukr 608 (n+ 1)t 

Ay, COS nt 
ceding, we find that when mis large it rapidly approaches a,4 el” an. 
Hence the series will be convergent so long as el?! < h, i.e. so long 
as |v | < log 1/k. 

If the series of Fourier coefficients is still more rapidly convergent, 
it is possible for the Fourier series to be convergent for all complex 
values of t. Such will be found to be the case when a, or b, approxi- 


» which is the ratio of one harmonic to the pre- 


Re” F 
mate to ie even when & is greater than unity. 
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sine and cosine series which passed through the tops 
of a number of equidistant ordinates and we showed 
(page 78) that when the number of ordinates increased 
indefinitely, the coefficient of each sine and cosine 
term moved to a definite limit, and so the finite 
number of terms became in the limit the infinite 
Fourier series. ‘The fact that the coefficient of each 
term approximates to a definite finite limit as the 
number of ordinates is made infinite, taken in con- 
junction with the convergency of the series, which 
we have proved above, is a satisfactory proof of the 
validity of Fourier’s expansion. 

If we took any infinite series of different functions, 


$,(4), ¢.(#) . . . . we could of course determine the 
constants in the expression 
Cybi(2t) + CoPa(t) + . « « + endn(2) (8) 


so as to make it pass through the tops of any given 
ordinates; but it would not in general happen that 
the ¢’s tended to any definite limit as m was increased 
to infinity. Nor would it generally happen that the 
graph represented by (8) would, when the ordinates 
were sufficiently close together, approximate to the 
straight lines joining the tops of successive ordinates. 
In general, the curve would make a half oscillation 
in running from the top of one ordinate to the next ; 
a fact which is consistent with, and explains, the fact 
that the expression (8) has no definite limit when n 
tends to infinity. The fact that the sine-and-cosine 
series on the contrary approaches a limit, means that 
the graph represented by it approaches a limit, and 
this limit cannot be anything else than the limit of 
the straight lines joining the tops of consecutive 
ordinates. 

We will now sum directly the whole of the first 
harmonics of Fourier’s expansion including the 
constant term. We have 
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fd) = 2+ acostt+... 
+b, sint+ . 


2m 
where ae \ f(4) cos nt dt 
a [e} 
1 (eee 
and Dipeos if f(2) sin nt dt. 
a oO 


Let us use the letter x for the variable in these 
integrals and, since the starting-point of the range is 
immaterial so long as the range is 2x, make the limits 
t—zaandt+a; so that 


I tr 
an = = \ (*) cos nxdx. 
Jim 


The apparent dependence of a» on ¢# in this result is 
illusory, for, as we have just said, the value of ¢ does not 
affect the integral. 

Replacing by, by a similar expression and adding 
together the two terms constituting the mth harmonic, 
we have 


I t-+1 
ay cos nt + by sin nt = = \ F(x) cos n(w — t)dx. 
t—t 
On writing « = u + ¢ this becomes 
wv 
= 1 { f(t + u) cos nu du. 
ve th 


Hence the sum, Sy, of the first m harmonics plus 
the constant term is 


i{ K(e+ Oe +cosu+cos2u+ ... + cos nud 
Replacing each cosine term by its exponential value, 


we easily find that the value of the series is sua (ee 
2 sin $u 
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Hence we have 


Sn, = a\ ne ae u) sin(m + 4) 9 


sin $u 


Now change the variable again by writing (w+ 4)u=v, 
and for brevity write n’ for n+ 4. We then get 


n't : 
iL 7) sin Uv 
S=2\s(+3). Sw... 
7 WL ros VU 
eye 2h sin 
ZH 


Tig7sis large; f(z + : , for moderate values of v 


say for v<V, is sensibly equal to f(z) unless f’(z) be 
very great ; also 2n’ sin = is then sensibly equal to v, 
hence, taken over for the range from — V to V, the 
integral is sensibly equal to 


VV: 
pet Ciaauey hy Aaa (10) 
Tt v 

Now this converges when V = , the contribution 
from large values of v being really less than seems to be 
the case, since each positive half-wave is very nearly 
neutralised by the following negative half, when v is 
large. If V is moderately large, then, the value of 
this integral is sensibly that of 


“sin v 

s 

\ === av; 
Soe 


which we have seen (page 11) is x. ‘There is clearly 
no difficulty in choosing a large finite value of V which 
will make the integral in (10) as near to x as we choose 
and then in taking n’ much larger than V so that the 
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approximations (of 2m’ sin ie v and f(t+ ) S616) 


in (g) are as accurate as we wish them. 
We conclude, then, that the limit of S,, when 1 is 
infinite, is f(#) provided f’(#) is finite. This furnishes 
a rigorous proof of Fourier’s expansion. 
If the function f(z) has a discontinuity in magnitude 
at the particular value of ¢ considered, the result (9) 
breaks up into two, giving 


ne : 
I VU sin VU 
oe = = \ f(: — )) dv 
Tt nN nes U 
fo) 2n sin —; 
2n 
fe) ‘ 
J Vv sin Uv 
+= f(2 — dv. 
Tt ; n pee ee 
—NT 2n es 7 


As previously, the first integral becomes in the limit 
co 


i as ) \ oe dv =4f(t+ 0), 


and the second integral becomes likewise 
@) 
f(t — 0) \ sin v ed ae 
ar aan Mis dv=4tf(t— 0); 


so that the sum to infinity of the Fourier series is 


FH f(t + 0) + f(t — o)}, 


at a point of discontinuity, and f(#) at other points.* 


* The student with a little mathematical courage will probably 
obtain this result by boldly putting n’ = cc in (9) when it becomes 


Sa =i | (6) 2% aoe™ fe, 


and this is not so reprehensible as some authors would have us believe. 
It is necessary, however, before writing oc, or any other particular 
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At a point of discontinuity this investigation clearly 
obtains what we have called the formal sum of the 
series; and we have seen before that at such points 
this sum of the series is midway between f(z + 0) 


and f(t — o).* 


value for a variable, to make sure that the expression assumes a single 
definite value when this is done. A case in point is the series 


sin 1/e-+ 4 sin 2/e4+- 4 sin 3/f+ .... 


As ¢ increases towards ++ oC, this series, as we have seen, approaches 


Chant TT . : . 
the limit of + —; while when ¢ is negative and moves towards — oc 
2 


the series approaches the limit of = while actually at t= +c the 


series is indeterminate and can assume any value from — 1°8519 to 
+ 1°8519. Such behaviour, however, is not difficult to detect, and 
in the case of infinite series, is always bound up with the fact that at 
such a point the series does not begin to converge till after an infinite 
number of terms. 

* It is clear also that the limit of (9) when n’ is oc will only be 


© fA 2 do ot fl) 
Tisch ae 
provided that f (: fe 5) is sensibly constant for all values of v below a 
n 


reasonably large quantity, say V, when n’ is very large. This is 
generally expressed by saying that f() must not have an infinite number 
of maxima and minima or discontinuities in a finite space in the 


neighbourhood of the value of ¢ considered; for if so, e435) 


could fluctuate in magnitude during a finite change in v when 1’ 
was infinite: and so we could tell nothing about the value of the 


integral 
v\ sin v 
= dv. 
iG ) v ‘ 


This is the reason why so many books dealing with Fourier’s series 
continually repeat the condition that the function must not have an 
infinite number of maxima or minima. We have generally omitted 
specifying this condition, since no practical function ever does behave 
insucha manner. Such behaviour is exclusively confined to functions 
invented by mathematicians for the sake of causing trouble. 
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EXAMPLES. 


By putting t=o in the formula (21) of page 30, and by 
Hose the a’s by their integral values, prove that 


af +42) +18) +... +afem = 2 | "poe 


=e 2/10 cos ntdt + *[ 0 cos 2ntdt-+ ... 


[The reason for writing 4/(0) + $/(2m) instead of either f(o) or f(2m) 
is in case the periodic function f(t) has a discontinuity at t= 0 which 
will of course be the case if f(2m) = _f(0).] 

De ane by repeated oe by parts, that 


Hones that all the differential coefficients are finite and the series 
is convergent. 
3. Prove, by substituting the result of Ex. 2 in Ex. 1, that 


[ foa= ® 4/00) Ee) dak a 


mips dices NO Paes 25 


pl OR se ae 


The numerical series in HN ee are tae to have the values 
os mon eal ears 0 
* 90" 945” 9450" 93555° 
mn By changing the range of integration in Ex. 3 from 2m to nh 
see that 


eee) 


. respectively. 


ie Sat = BAS) +O) + JOB) +--+ A fla) 
=a -FO} + ge MC'm—F")} 


Mh Vv (ann VII 
— glen —r'@} + 2 {perma —pmo)} 


= ae 


provided the series is convergent. The series can often be used when 
it is not convergent, provided the terms decrease to a small value before 


9 
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they finally increase, for it may be proved that in many cases the 
error made by stopping at any term is less than that term. 

5. By taking the limits of integration as w and oc prove that the 
result of Ex. 4 becomes 


[ Fede = PBs) +FO4+ D+ fe tat. 4 
2 
A 


supposing the integral convergent. 


a 


4 6 
FO apf WtghO— Sire +. 


6. By taking f(7) = = in Ex. 4 and the limits of integration as m 
and n, which are integers, prove that 


Me ee I fF eer 


nm mn a-+ 2 m—I1- 2m 


ae tee a) - 2G - 3) 
12\n2 m2, 120\n* m4 252\n® m8 
iPr I 
ts) ee 
Calculate log,2 to seven decimal places by taking »= 10 and 
m= 20. [Ans. *6931472.] 
7. Prove, from the result of Ex. 6, that the limit when m is infinite 
of : + : + ; ae! Biss . + = — log m is a constant. Show, by 


taking n = 4, that this constant is 577215 approximately. 
8. By taking f(t) = log ¢ and the limits m and n as integers, show 
that 


m log m/e — n log n/e = 3 log n+ log (n+ 1)+... 
+ log (m — 1) + flogm— E17 Saag 


2\m n 


eis eee 


1260\m>5 nd 


m I I I I I 
I 


Regarding m as a variable and m as a constant, show that this result 
can be written 
I I 
log |m = Const. + (m + 2) log m m +- 


12m  360m3 


I I I 
+ y360m® — 1680m? T 1188m° 
By taking m= 10, show that the constant is “91893850... It 
can be shown theoretically that this constant is 4 log, 27. 
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Deduce from the last result, that when m is large, | m approximates 
to a ratio of equality with 
ee ON 
Jama) 


g. From Ex. 5, page 73, obtain the following results 


sin2x , sin4x , sin 6x 


x 
—| log sin x dv =x log 2+ 3{ 12 “is 22 5 32 laa 
[e) 


es 


wheno<x<7T; 


x sin 2% sin4x , sin 6x 
— [log cos x dx = a log 2 — 33 — t nr e 
1} 
Tw 715 
when — =< 47 <= — 
D 2 


and 
sin2x” . sin6x . sin 10x 


C4 
— [Jog tan x de = 3 + 733 t 32 


eae). 


™ 
when o<w< a 

10. Prove, by differentiation or by replacing the sines and cosines 
by their power series, that the curve represented by the series 

COSI EP ECOS:2 6, ae COSE3 (6 
Teree rm oe 

is finite when 8 = oif p > 1 and has a finite radius of curvature at 
8 =oif p> 3. Also prove that the series 


sinB , sin 28 
3 “1 a ge 


has a finite slope at the origin if > 2 and an ordinary point of 
inflection at the origin if p > 4. 
11. If @(x) = o when x = 0, prove that the series, 


(x) , (2) ae 9(3) re 


I 2 3 


can have any value when x = o which the integral 


[22ax 


o * 


can assume, where X is any positive quantity. 


CHAPTER VII 
FOURIER’S INTEGRAL THEOREM * 


Fourter’s theorem, for a periodic function f(z) of 
period 'T, is that f(¢) can be expanded in the form 


(see page 29) 


fd) = “2 + a, cos 7 + a, cos 1 + she 
+ by sin T+ by sin Fe ee er) 
where 
an = Slo co dt... (2) 
and bn = i Hos dt ee ane) 


If we change the variable under the sign of integra- 


tion from ¢ to 24, we may multiply (z) by cos ia and 


(3) by sin ie and put these terms under the 


sign of integration; since they are mere constants 
during the integration with respect to 4. The two 


2nn _ 2nrt wee 
terms dp Cos ie + dy sin “rp constituting the nth 


* This chapter is intended for the student of physics; it has no 
engineering applications; but we hope it will not on that account be 
found uninteresting to all engineering students. 
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harmonic then become 
, 2am a | t) 1. 
alae IQ) cos 


e C3 SIF e 
We will write a =o and now suppose T increases 


without limit so that « becomes indefinitely small ; 
the above expression then becomes 


“| 42 cos no(A — #)dd. 


—t/o 


Now replace all the harmonics in (1) by expressions 
of this pee ; we get 


t/s 
hs fs yy + { f(A) cos o(A— t)da 


—T/G fue 
we K) cos 20(A —f#)dA + . 
Now we know that if F(x) is any function of x, the 
limit of o{F(0) + F(o) + F(2o)+ ... + F(x)} when 


di sO nis \ F(x)dx; and, since each single term is 


indefinitely small in the limit, it makes no difference 
whether the first term is oF(o) or $ oF(o). Applying 
this to the above case, we see that if we write o = da * 
and replace the limits +a/o by + 0, which they 


then become, we have 


Wes ‘| a\ F(A) cos «(4 — t)ddA . . (4) 
0) eS} 
This result is known as Fourter’s I ntegral Theorem. 


* This is quite a different use for « than that in the earlier part of 
this book: we do it for the sake of uniformity with other writers. 
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The theorem, of course, only holds provided the 
integrals are convergent. The integral 


ioe) 
\ fd) cos (A— tdi. 2 2. 65) 
==0C 
will necessarily be convergent provided 
(o @ 
| #4) | da 
—c 


is convergent. We will accordingly suppose that this 
condition is satisfied. This requires that f(A) should 
decrease faster than 1/4 when 4 is large, for both 
positive and negative values. ‘The student will also see, 
in a similar manner to that adopted in the last chapter, 
that if f(A) is finite everywhere, but has discontinuities 
in magnitude, each integral with respect to 4 in the 
equation 


Oe a [cos at 1 cos aA dA 


+ sinat\ f@) sin af dil da, fuk) ah) 


in which form (4) may clearly be written, 
will tend to proportionality with 1/« when « is large ; 
and if f(A) has discontinuities in slope but not in 
magnitude, each integral will likewise tend to propor- 
tionality with 1/«?, and so on. ‘The integration with 


oC 
respect to « is always convergent when | | f(A) | dA is 
—co 


convergent owing to the cos «# or sin «t factors, which 
virtually make the integral a series whose terms are 
alternately positive and negative such that, when « 1s 
large, each term is less than the one before it. 

L 
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Like Fourier’s series, Fourier’s Integral Theorem 
applies both to analytic or artificial functions, subject 
to the integral 


\ 1 | da 


being convergent. 

We can have functions of two distinct types; either 
functions gradually increasing from zero at t= — © 
and decaying to zero again as ¢ tends ¢ + © ; or those 
which are zero from t= — © to t= 4, say, when they 
suddenly become finite and remain finite from t= 4 
to t = tg, say, when they suddenly become zero again 
and remain zero up tot = o. A function of either of 
these types is conveniently called a “ disturbance ” or 
a “Pulse.” It is important to bear in mind that it 
must necessarily be defined from t= — ~ tot= + 
or Fourier’s Integral Theorem cannot apply to it. If 
iwe take the limits of 2 as ¢, and #, in Fourier’s Integral 
Theorem, instead of — o and + , we are simply 
applying the theorem to the pulse which is equal 
to f(t) from t='4 to t= 7%, and Zero woutsidesthe.- 
limits. 

Fourier’s integral theorem takes a simpler form 
when f(t) is either an odd or an even function, as might 
have been expected from the behaviour of Fourier’s 
series. If, for instance, f(7) is an odd function of 2, 
the first integral with respect to 4 in (6) vanishes ; 
since f(A) is an odd and cosa an even function of 4 
and the product is odd: also f(A) sin a4 being an even 
function, the integral from — o to +c is equal to 
twice the integral from 0 to «. 

Hence, if f(#) is an odd function of t, we have 


fi) = 2 { sin wl \ 4 sin ad al da. (7) 
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Similarly, if f(z) is an even function of ¢ we have 


HO 2{ cose] | cos a ai da. (8) 


If we write 


| 10 sin ai di = 4(a), 


and 
\ f() cos ad dd = y(a), 
re) 


these results may be written in the following form : 


diate \ 4 sien ere | | ms 
fen Ko sin at da =" f(2) | 
and if \ © Corea) | 

then jx cos at da = 7 f(i | nase 


In the first of these f(A) is odd, and in the second it 
is even; and, since we see that the relation between f 


and ¢ and f and y is reciprocal, (save for the factor =) 


we see that, if we wish to interpret > and p for negative 
values of a, ¢ is to be taken as odd and yw as even. 

oran example, let f(t) = ¢-4 from t= 0 to t+ ox. 
Then we have by (9) 


\ a sin addi = 9(«), 


fe) 


- a A a 
Dut tis: inteoral i¢ ————_; 
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whence we infer that 


iS a It 
+——5 sinatda = —e-s 
cata 2 


or, changing symbols, that 


CS . 
\ %SiN Me pT yom, 
oO 


Similarly by (10), since 


C 


oe 
=6X ee eee, 
\ e~° cos ah di ree 


oO 


we have 


oo} 

(A It 

\ +—; cos at da = —e-t; 
Pe alg i 2: 


or, in other symbols, 


~ cos mx a 
peepee —_ - e-em, 
tae os 25 


Fourier’s integral theorem in the forms (9) and (10) 
thus becomes a useful theorem in the integral calculus. 
The way in which Fourier’s integral theorem per- 
forms the “‘ harmonic analysis” of a pulse can be 
seen by the following considerations. ‘The ratio of the 


period of the mth harmonic to the ( + 1)th is ate 
n 


> 


which tends to unity when is infinite. Also, when the 
period is made infinite, all harmonics of finite period are 
necessarily harmonics of an infinite order; they are 
thus infinitely close together, or in other words, there 
are an infinite number of harmonics whose periods lie 
between two given finite limits. Now the equation 
(6) can be written 
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I Cc Cc 
fee | [| Heo! a] sone 
1) Ses 
I CS co 
ae | [ \ F(A) sin ad au sin at da 
a (@) OS 


Dropping the integration with respect to a, we see that 
I 


se [| T(4 cos ad a cos at da 


7 


a's 1) sin ad a Sie. d done wee Te (11) 


is the part of f(¢) which consists of harmonics whose 
? ; 2 
periods lie between “* and ae or whose fre- 
a a+ da 
a + da 
aa 


quencies lie between = and If we denote 
if A 


the integrals within square brackets by P and Q 
respectively, (11) can be written 


C=) cos at + (ee) sin at. 
It Ut 


Here a is not quite constant but varies over the 
small range da from a to a + da, so we see that the 
amplitude of the cos at term at a may be regarded as P/ax 
per unit range of a or 2P per unit range of frequency. 
Similarly, the amplitude of the sin at term at a may 
be regarded as Q/z per unit range of a, or 2Q per unit 
range of frequency; and the amplitude of the com- 
plete harmonic as 2./P? + Q® per unit range of 
frequency. 

It will now, we think, be best if we adopt the 
language of that branch of physics in which this sub- 
ject has its chief applications. ‘The pulse is generally 
supposed to be some disturbance travelling through 
a medium ; it may be, say, an air pulse caused by an 
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explosion, or an electromagnetic pulse due to the 
stoppage of an electron by collision. In the first place 
we should take f(#) as proportional to the condensation 
in the medium at any given point; and in the second 
case to the electric (or magnetic) field at the point 
considered. In either case the energy per unit volume 


of the medium at time fis proportional to f()?, and so we 
(oa 


may take \ f(7)?dz as a measure of the energy of the 


whole pulse.* 

Now on page 34, in finding an expression for the 
mean square of the value of a periodic function in terms 
of the coefficients, we were really expressing the fact 
that the energy in the function was equal to the sum 
of the energies in each of its harmonics. An exactly 
similar result holds in the present case, which we will 
now obtain. ‘Taking the equation 


oc cc 
Or = af F(A (cos «A cos aé + sin aA sin at)dd 
te) —— OC 
or, with our previous notation, 
(i Bee 
bee “| dx (P cos af +_Q sin a?) ; 
oO 


let us multiply both sides by f(#)dz and integrate from 
—e« to + «. Remembering that P and Q are func- 
tions of a only, we have 


| farar 2 x| a [P\ 0 Sine aff ane at] 


af ox = Ae + Q%)de., 


* ‘There is no need to introduce the factor 4, which the student may 
be tempted to supply, as we have dropped other constant factors 
depending on the properties of the medium. 
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This shows that the total energy may be divided 
up among the different values of «, the energy lying 
2 


9 


a 


between « and « + d« being a ‘da; that is, the 


P? + Q: 


energy per unit range of « at « is —————; or, the 
es 
energy per unit range of frequency at the frequency 


a . . 
57 18 2(P* + Q*). Since the energy of an oscillation 


is proportional to the square of the amplitude, this 
result is consistent with the one arrived at on p. 149. 


When the different frequencies in a pulse are separ- 
ated in any manner the result is called the “spectrum” 
of the pulse. The light given out by a piece of red- 
hot metal or carbon, for instance, consists of the 
“‘ pulses” produced by the collision of the electrons 
with one another or with the atoms of the substance. 
There is not necessarily any periodic motion whatever 
associated with the origin of the light from an incan- 
descent solid body. A glass prism or a diffraction 
grating analyses this succession of pulses into its 
constituent harmonics in exactly the same manner as 
Fourier’s integral theorem analyses it. The fact that 
the spectrum, in this case, shows the presence of 
harmonics of all frequencies within the range of the 
instrument, shows not only that the origin of the 
light is in pulses, each containing all frequencies in 
itself, but that these succeed one another at irregular 
intervals. For a succession of exactly similar pulses 
at exactly equal intervals of time would obviously 
constitute a periodic disturbance, from which all 
harmonics would be absent except those whose periods 
were sub-multiples of the given interval; hence the 
spectrum would be a “bright line” spectrum and 
not a continuous spectrum. The distribution of 
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energy in a succession of similar pulses can only be 
the same as the distribution in each pulse if they 
follow one another at perfectly random intervals of 
times 

If we make a slit in the screen on which the spectrum 
is focussed we can let through any desired small range 
of frequency at any point, and by absorbing these rays 
we can measure the energy in this small region of the 
spectrum. In this manner we can get experimentally 
the distribution of energy in the spectrum.* 

We have seen above that if we knew the form of the 
pulses we could calculate this distribution of energy. 
We shall now show that, knowing the distribution of 
energy, we can find the shape of an odd pulse, and 
also of an even one, which would give such assigned 
distribution of energy in the spectrum. ft 

Let {F(«)}2d« be the energy in the spectrum be- 
tween aanda+ da. ‘Then we have xF(a«)? = P? + Q?. 


Accordingly, if the pulse 1s even so that Q = 0, we 
Cc 

have /aF(«) = P = 2\ F(A) cos «A di, where f(#) is 

the form of the pulse. : 


Equations (10) at once give us the result that 


| Fe cos at dt = 7 f(i 


or py val FO) cos at da. 


* It is, of course, necessary to correct the result for the light lost 
by reflection or absorption in the prism or grating. 

+ The student familiar with radiation from electrons will recognise 
that if an electron is projected against another and is brought to rest 
by the interaction and repelled back along the line of advance, the 
pulse will be an even one; while if it had passed on without being 
brought to rest the pulse would have been an odd one. This is easily 
seen, since the disturbance at a great distance is proportignal to the 
acceleration of the electron. 
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Similarly, for an odd pulse its form would be given 


by 
1 re 
oe Jal Fe) aoe de. 


Example. Lord Rayleigh proposed the expression C2a2e~2c«dw for 
ehes the 


constants C and ¢ depending upon the temperature of the body 
emitting the light. Here, F(«) = Cae-, so that an even pulse to 
give this distribution would be given by, 


: wu 0 
the energy in the spectrum between the frequencies — and 
2T 


G io) 
f= Jal on cos at dx 


and an odd pulse by, 
Cc ioe] 
G) = =| ae sin at da. 
AO = al. 


Multiplying the second by 1/ — 1 and adding it to the first we get 


| ri re ee oe ok CG 
FFA = Fal oe am ae re: 


T 
And so 
_ Ce — 79) 
Fo(2) = V7 (c2 ae 72)2 
and 
ZCzr 


A® = Vn (2 rs #2)2 


Similarly, we could find the shape of the pulses to produce any other 
law of distribution of energy in the spectrum, such as Planck’s, whose 
formula agrees much better with experimental determinations than 
the one used above; but this calculation has not any very great 
physical value, since white light does not consist of a succession of 
pulses of exactly the same shape, so what we have arrived at is only a 
a sort of mean shape of the pulses. 
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EXAMPLES. 


Todt | F(x) cos nxdx = @(n) 
eo} 


prove that 


[ cos nxdx = $¢(0) + 9(I) cos # + 


@(2) cos 2n-+ 9(3) cos 3n-+ . 
provided o <n <7. 


[For within this range each side is equal to =F (n).] 


2. Prove, similarly, that if 


| BEG} Gn aed 


then 


[ve sin nxdx = (1) sinn + (2) sin 2n-+ (3) sinju-+ ... 


provided thato <n <r. 
3. Verify the result in Example 1 when F(x) = x — x. 
4. Prove that the energy, lying between the frequencies «/27 


and («-+ da)/27, in the pulse given by f(t) = 
tOee 3 


5. Prove, similarly, that for the pulse f(#) = e~?*, the distribution 
of energy is proportional to e~2*/2?*, 
6. Prove, similarly, that for the pulse f(#) = sech pt, the distribu- 


: ; : Teo 
tion of energy is proportional to sech® —, 
2 


c : : 
=——— 18 proportional 
ao es 


7. Show that the distribution of energy in the pulse, which is 
equal to unity from ¢=-—1I to t=1 and to zero outside these 


Pere 
sete ; sin % 
limits, is proportional to (a) : 
om 


8. Show that the distribution of energy in the pulse, which is 
zero when ¢ is negative and equal to e~* sin pt when t¢ is positive 


(where & is very small compared with 9), is very approximately given 
by 


I 
BBE (a = p)* 

Hence note that a damped oscillation, when examined in a perfect 
spectroscope, will not show a perfectly sharp line, but will show a 
maximum of energy when « = 9, and the energy will be appreciable 
so long as « — p is not many times greater than &. 
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9. Show that the forms of the odd and even pulses, which will 
ae ‘jie anal ( supe of energy over the frequencies between 
«% — 0)/2m and (a + 6)/2m, are 
Pore ein a 
and to) = (= =) COS at 
respectively. 
to. The energy lying between the frequencies «/27m and («+ dx)/2x 


is proportional to (a2/-2e-22)dx ; prove that the forms of the odd and 
even pulses to give this distribution of energy are 


i ame g/D 
A ae ae 


and 
__ cos (p tan™! ¢/b) 
Si = (b? + 72)p/2 
respectively. 
11. If f() represents any convergent pulse, prove that the succes- 
sion of pulses 


2 Eft — 2m) + f+ f+ om) + fe+ant... 


is represented by the Fourier series 
fe + a, costo. ee 
2 
+ b,sint+... 


where 


rites 
Ae | Fd cos nt dt 
mJ —-e 


im Ree 
ya b,(/ooe | f(d) sin nt dt. 
Us I} anise 


CHAPTER VIII 


THE SEARCH FOR PERIODICITIES IN OBSERVED DATA NOT 
PERIODIC 


Tue physicist, astronomer, or meteorologist— 
especially the latter two—frequently obtain curves 
which are suspected of containing one or more periodic 
components. If the curve is found to be definitely 
periodic, and the observations have extended over a 
sufficient number of complete periods to establish 
the fact beyond doubt, recourse can be had to the 
methods of harmonic analysis previously given when 
as accurate an analysis of the curve can be obtained 
as its own accuracy permits of. More generally, how- 
ever, the curve contains several harmonics the periods 
of which are not in simple ratios to one another, so 
that the period within which the motion recurs is very 
large—if, in fact, it recurs at alli= Forgimstancess: 
curve composed of two harmonics of periods 12 and 
17 days has a period of 204 days, but if the ratio of the 


periods, instead of being 1°416666 . . . , was in error 
by about 1 part in 760, so that the ratio of the periods 
was actually 1-414213 ... or »/2, the “period,” 


instead of being 204 days, would be infinite; since 
there is no common multiple of two numbers bearing 
this ratio to one another less than infinity. Thus to 
calculate a period in cases like this, where the periods 
of some of the harmonics are only known to a limited 
degree of accuracy, is to tread on far too uncertain a 
ground, 

We have thus to determine the presence of a har- 
monic constituent in the curve, or set of observations, 

156 
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in the absence of any indications of a definite period. 
This can be done by the application of Fourier’s 
integral theorem explained in the previous chapter. 
Let the given curve or observations be represented 
by f(@) from t = — T to t = T, and let us apply the 
theorem to the pulse defined to be equal to f(#) from 
t= — T'to t= T and to be zero everywhere outside 
this range. Then we see from page 151 that the 
energy in this pulse lying between the frequencies 


a. a + da. : 
Bande" is proportional to 
27 27 


([\ @ cos a di a [| @ sin at di] [da ober (1) 


or, say, R*da. These integrals can be evaluated by 
Simpson’s rule or other suitable method of quad- 
ratures; such evaluation must be done for a fairly 
large number of values of a. 

If R? is plotted against a the result is called the 
* Periodogram ” of the original curve or set of obser- 
vations. We thus see that the periodogram is nothing 
other than the distribution of energy in the spectrum 
of the pulse defined as being equal to the observations 
over the range for which they extend and zero outside 
that range. Practically it is advantageous to split 
any given function into its odd and even components ; 
R for either of these functions is then either 


Ag 
a) Hi ICOS UN GN a (sa bs 3 00( TA) 
eh 
or 2 Fa AGh ce Bone wt ae DD) 


It is further advantageous to plot |R], instead of R, 
against a and to call these diagrams the periodograms 
of the functions in question. 
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We shall now show that the periodogram possesses 
pronounced maxima at each value of a which corre- 
sponds to a harmonic constituent in f(t) provided that 
the range of observations extend over several such 
periods. As an example, let us suppose f(#) consists 
of the harmonic c sin pt only and that the observations 

N. Na 
extend from — ” to + —" so that they cover N 
periods. In this case the first integral in (1) vanishes 
and we have 


Nr/p 
R= 2c\ sin pA sin ad da 
@) 


=— 5-/sin (p — a)Na/p = sin (p ac Bs (2) 
Ee pta 

The first term attains its numerical maximum value 
of 2cNz/p when a = p, while the second term neces- 
sarily lies within the limits + 2¢/(p + a), and lies within 
the limits + c/p in the neighbourhood of a = p; that is, 
it is only +1/2Na of the maximum of the first term. 
The first term vanishes when (p — a)Nz/p is a multiple 
of x other than zero: that is, it vanishes when 
a= p(1 +1/N), ora = p(1+2/N), etc. It will have 
subordinate numerical maxima at roughly half-way 
between these points; the value of the maxima when 


(a — p)/p = + ( + 4)/N being 


= 2¢ = 
(n + 4)p/N 
approximately, which is oat of the maximum 


when a = #. 

Fig. 22 shows a graph of |R| when N= 5; the 
extent to which the values of |R| are uncertain, 
owing to the neglect of the second term in (2), is 
shown by the distances between the two dotted curves 
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and the line midway between them. It is clear that 
such a relatively small uncertainty can do nothing to 
mask the subordinate maxima, let alone the very 
pronounced central maximum. 


If we had considered the harmonic ¢ cos pt instead of 
c sin pt, we should have had 


Nr/P 
Riz 2c\ cos pA cos ad di ; 
Oo 


which would have given the same as (2) with merely 
the sign of the second term changed, so that the graph 
of |R] would still lie within exactly the same limits as 
previously. It should also be noted that all these 
results hold whether N is an integer or not. 

Let us now suppose that f(t) contains two harmonics 
of different periods, say, csin pt+c’sin p’t. If the 
two periods p and p’ are widely separated there 
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will be no difficulty in recognising the two central 
maxima corresponding to a=p anda=p’. If p 
and p’ are nearly equal, the two maxima become 
confused ; but if the amplitudes of the two harmonics 
are at all comparable, separation of the periods is 
always quite easy when the two central peaks stand 
just clear of one another: that is when (p’ — p)/p = 
+2/N; or when the number of periods of p’ com- 
prised within the range of observations differs from 
the number of periods of p (which is N) by two. 
Hence to separate harmonics with nearly equal periods 
an extended series of observations is necessary. 

Separation is generally easy when the maximum 
of one central peak coincides with the zero of the other ; 
so that one of the harmonics has only one more period 
in the range of the observations than the other. In 
this case, if c and c’ are of the same sign and com- 
parable, the two central peaks will combine into a 
single peak whose width is three times that of the 
subordinate peaks. ‘The subordinate peaks for the 
two periods will, in this case, reinforce one another 
and so will be prominent. If c’ = — ¢ the peak for 
p would appear negative if we plotted R but since 
we plot | R | it will readily be seen that the combination 
of the two peaks will give two equal peaks each of 
a width 1-5 times the width of the subordinate peaks. 
In this case the two systems of subordinate peaks will 
largely neutralise one another and so may be indis- 
tinguishable ; nevertheless their width is known theor- 
etically (viz. a variation of 1/N in the ratio «/p), and 
the presence of two equal peaks each 1-5 times this 
theoretical width is conclusive proof of the presence 
of two harmonics of nearly equal periods and of 
amplitudes of opposite signs. 

In our search for periodicities by this method we 
are not interested in the exact shape of the graphs: 
we merely use them for the positions of their maxima : 
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it follows that we could use instead any other curve 
which possessed equally distinct maxima in the same 
positions. A much more easily obtained curve, and 
one which is quite as satisfactory, can be found by 
replacing the cosad and sin ad factors in the trouble- 
some integrals 


T ae 
\ f(A cosaddi and \ F(A) sin ad da 
° o) 


by another periodic curve of the same period which 
will make the integrations easier. Such a curve is that 


Fic. 23. 


in Fig. 23, which is alternately + 1 for half a period 
and then — 1 for half a period. When the origin is 


chosen so that it is an odd function, its equation is 


in 3a’ sin sad 
$,(aa) = 4 {sin ad + mie a: ce eye | 


and when arranged as an even function, it is given by 


A 7} | 
$,(at) = 4 {cox an — SSE 4 OES 


Substituting ¢,(a4) for sinad in the second of the 
given integrals gives us 


Meee fe 4) 1. °c 
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while the other one becomes 


n/2u p3n/2% e5n/2u it 
[ie \ Ree ee T Za ] oat. 
fo) n/2a J 3r/2% 


Each of these component integrals is very easily 
obtained, since we have the graph of f(A) (or can easily 
get it by plotting the observations) and from this we 
can read off the ordinates required in evaluating the 
integrals by Simpson’s rule. How the graphs of |R| 
are modified by using these functions ¢, and ¢, instead 
of the sine and cosine may easily be seen from their 
harmonic analysis. When using the function ¢,(a4) 
instead of sin a/, the first term in (2) becomes 

4 (sin(p—a)Na/p , rsin(p—3a)Na/p ) 
Wi oes 3B 


instead of 


sin (p — a)Na/p 
fa 
the second term of (2) being just as negligible in this 
case as in the former case. 

Thus the shape of the maximum in the neighbour- 
hood of a = p is unaltered, but we have introduced 
into the graph for |R| another maximum at a = /3 
which is 4rd of the height of the maximum at a = p 
and which need not correspond to any frequency in 
f() at all. Similarly, there is another maximum at 
a = p/5 which is 1th of the maximum at a = p and 
which also does not exist in f(#). It follows, therefore, 
that we can use this graph just as the correct one, | 
provided we neglect these maxima at sub-multiple 
frequencies. Should, however, the maximum at 
a = p/3 be much more than 4rd of the maximum 
at a = #, it is clear that a harmonic of the frequency 
corresponding to a = p/3 must exist in (2), since 
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there is more of this frequency present than we have 
introduced by using the functions ¢, and ¢,.* 

We may vary the functions ¢, and ¢,, so as to reduce 
still further the labour of evaluating the integrals 
involved, without affecting their effectiveness in 
detecting the various harmonic constituents. If we 
add unity to each of these functions, for instance, the 
result will only be to add to |R| the constant quantity 


ae 
\ f(4)d1; and will therefore make no difference to 
—T 


the peaks of |R|; yet it will have reduced the labour 
of finding |R| to one half. Hence we may use for 
¢, and ¢, functions which are zero for half a period 
and unity for the other half. It can readily be seen 
that we can go further in this direction and use 
functions which are zero for three quarters of a period 
and equal to unity for the other quarter. Going to 
the limit, we see that we can use functions which are 
zero throughout the whole period, save in the imme- 
diate neighbourhood of one point. In this case the 


trouble of integrating ¢,(«A) . f(4) and 4,(aA) . f(A) is 


* The student of optics will readily see the optical interpretation 
of this. The first periodogram corresponds to the analysis of a pulse 
by a grating which only produces a spectrum of one order—the first. 
The second one corresponds to a grating which not only produces a 
first order spectrum but a third order one of jth the energy of the 
first one and a fifth order of },th the energy and so on, the image 
in the third order underlying the image in the first order correspond- 
ing to }rd the frequency or three times the wave-length. No 
diffraction grating, only giving one order of spectrum, has ever been 
able to be constructed, although it is theoretically possible. We 
have to be content with one behaving more or less in the second 
manner. There is, similarly, no reason why we should insist on 
employing the analogue of the first case in constructing the periodo- 
gram. ‘The student may also note that if the intensity of the image 
in the first case is unity, the intensities of the different order images 
in the second case are 8/n?, 8/gn?, 8/25n.. . which add up to 
unity. 

M2 
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done away with, as the result simply reduces to the 
sum of those values of f(A) for which the function 
¢, or $, is not zero. This is made the basis of the 
practical method described in the next Section. 
From the above we deduce the following method 
for searching for periodicities in any given series of 
observations. We take the origin of time at the 
middle of the range and represent the observations as 
the sum of an odd and an even function by the formula 


fa) _ ft Fak St 4 FY) ns ek 


For each of these functions we plot |R| and find its 
central maximum peaks. Each of these corresponds 
to a definite harmonic in the series of observations. 
In this manner all cosine harmonics are found from 
the even function and all sine harmonics from the odd 
function. 

Having found from the two graphs for |R| the 
frequencies that exist in /(#), we must next find their 
amplitudes and phases. Let a, be the value of a 
corresponding to any frequency present in f(#). Let 
M, be the central numerical maximum at a, of the 
graph for |R| constructed with the use of the cos ad 
factor for the even component of f(/), and let M, 
similarly denote the maximum for the odd component 
of f(z) constructed in a similar manner with sin aA. 

Then if a, cos a,t + 6, sin a,t are the two correspond- 
ing harmonics in f(t), we have from (2), (page 158), 

Mis Ae 
\ ay 

But 2Nz/a, is the total range of time, T, over which 
the observations extend. Hence we have the simple 
relation 

CVs 
similarly, 4, = NET. 
If the graphs for | R| had been constructed with the 
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aid of the functions ¢,(a4) and 4,(a4) and if Mj and 


M, had been the values of their central maxima at 
a, respectively we should have had 


GNA Ts ts, 5) 
and Praia 5 (6) 

We must similarly write down the harmonic terms 
corresponding to any other plainly indicated central 
maxima flanked by their subordinate maxima and 
minima which occur in either of the |R| graphs, 
neglecting, as we have explained above, central 
maxima corresponding to sub-multiple frequencies if 
the functions ¢,(a4) and ¢,(a/) have been used unless 
such maxima are unduly prominent. 

Having found in this manner a number of periodici- 
ties in f(t) and the amplitudes of the sine and cosine 
components of each of them, we may subtract the 
sum of all these terms from f(t), when we shall be able 
to see whether the residue consists of periodic terms 
whose periods were outside: the, region tested, or 
whether it consists of a non-periodic function, or 
whether it appears to be nothing more than the errors 
of the given observations. 


A Simple Way of Applying the Method in Practice. 


The foregoing method may often, in practice, be 
greatly simplified—so much so that it hardly seems to 
depend on the above theory—by making use of the 
fact that if we add the values of a function for equi- 
distant intervals, all harmonics, save those for which 
the interval is a period, tend to cancel out.t 

Let us suppose that the given observations can be 


* Similar easily obtained results will hold when any other periodic 
functions are employed in place of @, and 9. 

+ If the accuracy of the given observations is poor they should be 
plotted, and smoothed values read off from the graph in the ordinary 
way. 
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represented by several harmonic terms together with a 
slowly varying residual (if any), so that we may write: 


f(t) = ¢, cos (ayt — 1) + Cy COS (gt — Ye) +... + (2). 


If we form the sum 
fat ser) ++ #)+. 2 @ 


when a = a, taking as many terms, say N, in the sum 
as the extent of the observations permit, the first 
harmonic will become Nc, cos (a,;t — y,). The other 
harmonics will be in different phases in different terms 
of this sum unless their frequency is a multiple of a, ; 
and so their resultant will in general be small com- 
pared with that of the harmonic a,. The effect of 
any harmonic, say, ag, will not be small if its extreme 
phases in the sum (7) do not differ by more than a 
from one another, for then every term has a positive 
component along the line of the mean phase. This 


270 


will be so whenever a, =) differs from a multiple of 
a 


1 

2n by less than z/N; that is, when a,/a, differs from 
an integer by less than 1/2N. On the other hand, if 
the total range of the observations—that is, 2%N/a,— 
is an integral multiple of the period of any other 
harmonic, such harmonic will absolutely cancel out 
in the sum unless its period is also a sub-multiple of 
a,. he effect of the slowly varying function y(A) 
on the sum considered will only be to add to the 
harmonic Ne, cos (a,t — y,) another slowly varying 
function which can do little to obscure it. 

If this sum be plotted against ¢, the result is thus 
approximately the sine curve Ne, cos (a,t — ,) plus 
something which is approximately constant. The 
amplitude of this harmonic can be found by measuring 
the total variation of height, V, of the graph which 
is equal to 2N¢ ; also, if the position of the maximum 
value is at ¢ = #,, we have y, = ajt;. 
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If this curve is not very approximately sine shaped, 
the presence of harmonics whose frequencies are 
multiples of a,/2x are indicated. If this is so, these 
harmonics can be determined by subjecting this curve 
to the usual process of harmonic analysis. 

The only difficulty about this method of analysis is 
to know when the value we take for 2m/a in forming 
the sum (7) agrees with the period of one of the 
harmonics present in f(?). 

If we use a value of a slightly different from a, the 
term ¢, cos (a,t — y,) will not be in the same phase 
in all the terms of (5); and if 2xNa,/a were equal to 
(N + 1)2x, the sum of all these terms would vanish, 
since their phases would be uniformly distributed over 
arange of 2x. So it follows that the range of variation, 
V, has a central maximum when a =a,. In fact, 
the variation is clearly according to the now familiar 
culve) — —*; or if we always consider V as positive 
(and it is not easy to tell when it should be considered 

: ; ‘sin x 
as negative) according to the curve y = | : 


Rs 

We accordingly work out the sum (7) for several 
different values of ¢ for each different value of a taken 
and plot the total range of variation V for changes of 
t against a (or 2m/a) and pick out each central maxi- 
mum. When drawing a graph through the points of 
this curve, unless we clearly see that one particular 
point lies right at the peak of one of the central 
maxima, we must re-work the sum (7), for the value 
of a which the graph indicates will be the highest 
point of the peak; for since the peaks are steep (the 
steepness increases with the time over which the 
observations extend), a point only a little way off the 
true centre may give a value of V much below the 
maximum value which is sought. 

Having found these peaks, we can easily write down 
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the harmonic term corresponding to each peak. 
Subtracting the sum of these terms from f(#), we can 
then examine the remainder and see whether it is 
constant, or some slowly varying (generally a decaying) 
function; or whether it consists of harmonic terms 
whose periods lay outside the region examined, or 
perhaps were missed, by taking too great a gap between 
the different values of a in the graph for V. 


An ItiusrratTivE EXAMPLE OF THE SEARCH FOR 
PERIODICITIES. 


We will suppose that some phenomena have been 
observed at 120 equal intervals of time (which we will 
call one day) with the following results : 


Day. f(é). Day. f(t). Day. Ft). Day. Ff(d). 
@ ut) ~ 3 62 24 93 14 
I 25 32 9 63 25 94 24 
2 28 33 17 64 25 95 34 
3 30 34 27 65 25 96 32 
4 30 35 37 66 17 97 29 
5 2 30 36 67 Io 98 22 
6 21 37 32 68 6 99 14 
Hi 14 38 25 69 5 100 4 
8 II 39 17 70 5 IOI 4 
9 9 40 7 7X 15 102 5 

1100) 9 41 Z 2 24 103 9 
Up 19 42 8 73 30 104 14 
iy 28 43 I2 74 34 105 23 
13 34 44 18 75 35 106 24 
14 Sy! 45 26 76 25 107 23 
15 39 46 26 aa, 15 108 22 
16 29 47 26 78 6 109 18 
17 18 48 24 79 fo) IIo 12 
18 10 49 21 80 —4 IIt 15) 
19 3 50 15 81 4 DE2 13 
20 ° 51 16 82 14 113 13 
21 8 52 16 83 24 114 14 
22 18 53 16 84 38 115 18 
23 28 54 17 85 39 116 14 
24 37 55 oT 86 33 LZ 13 
25 43 56 17 87 24 118 13 
26 30 57 16 88 14 119 13 
27 28 58 16 89 4 120 12 
28 18 59 16 90 =F 

29 7 60 14 gI —I 

3° =i 61 20 92 5 
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Here we see successive maxima at intervals of 10 or 
11 days, so there is certainly a component harmonic 
whose period is about that. We will test first of all 
for periods of 8, 9, 10, 11, 12, 13 days respectively. 

The work for the 8-day period is shown in the table 
below. We have only made use of the observations 
extending over Io periods, since the observations do 
not permit of us using a greater number in the case 
of the longer periods. 


TABLE A. 

19 25 28 30 30 29 21 14 
II 9 9 19 18 34 37 39 
29 To 10 3 fo) 8 18 28 
37 43 36 28 18 7 1 3 
9 a7), 2 SHE 36 32 25 wy 
77h 7 8 12 18 26 26 26 
24 21 15 16 16 16 17 21 
7, 16 16 16 14 20 24 25 
25 25 17 pce) 6 5 5 15 
ee ge 34 35 25 15 6 © 


202 211 200 206 I9t 185 178 188 


igeskanec, V, 16-211 — 178 = 32. 


Proceeding in this way for the other periods, we get 
the following results : 
Period 8 9 10 II 12 13 
ESN? 33 32 269 35 202 33,5 


No satisfactory graph can be drawn from these 
irregular figures, so we must test similarly for inter- 
mediate periods, say, of 8:5, 9:5, 10°5, 11°5 and 12:5 
days. ‘To test for 8§ days we write down in the first 
column the value of f(t) for 0, 8, 17, 25, 34, 42, 51, 
59, 68, 76 days respectively, and in the successive 
columns for 1 day greater than each of these numbers 
successively. The intervals between these numbers 


* In this last case we can only test over a range of g periods and we 
get V= 34. We therefore increase the result in the ratio 10: 9. 
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are 8 and g days alternately. Had a complete graph 
of f(4) been available for the whole of the 120 days, 
we could have picked out ordinates at intervals of 
8-5 days for this test, but ordinates at intervals of 
8 and g days alternately are quite as satisfactory, since 
the difference of phase from the mean phase is only 
2n/17, the cosine of which is -93. 

Similarly, if we were testing for a period of 83 days, 


ESeese 
Fl 7 ah 
Ca 
a 
ili ae 


eee \ 
CAST ao 


Fic. 24. 


we should write down in the first column the values 
of f()) forests 0,58, 16,925,533, 41 56cO ec cecas 
sively, where the intervals are 3,8, 9, S00 eum 

In this manner, we obtain the following results : 


Period 85 9°5 10'5 II'5 12°5 


Veg 58 112 167 184 116 

Fig. 24 shows the graph obtained by plotting V 
against the period in days. 

There is clearly one period of 10 days very approxim- 
ately and another about 12 days. Since the position 
of this latter maximum is a bit uncertain, we calculate 
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V for a period of 11} days, which gives V= 199. We 
then infer a period of about 11-9 days for this second 
peak. 

The heavy dotted curves show the secondary maxima 
corresponding to the peak for the 10-day period. 
This makes it clear that the distortion of the peak for 
the 12-day period is largely due to the first secondary 
maximum of the 10-day period peak. We may there- 
fore take 12 days for the second period. We cannot 
possibly expect to find the period to a greater accuracy 
than this ; for the difference is only 1 day in the length 
of 10 periods, which is only a difference of one-twelfth 
of a period in the whole range of the observations. 

We next plot the sums of Table A for the 1o0- and 
12-day periods against 7. 

For the 12-day period the sums are: 


Ue Sum. i, Sum. ioe Sum. 
fe) 269 4 181 8 gI 
I 291 5 138 9 139 
2 2606 6 92 Io 176 
3 241 7 89 II 235 


These when plotted give sensibly a sine curve of 
which the equation is 


y = 100 sin (= = =) 
is 


superposed on a gradually decaying function; but 
since each ‘‘sum”’ is the sum of ten items, the equation 
of the harmonic in f(#) is 


; ant It 
y = I0s1n (Z 1. *): 


Since the range of observations used in obtaining 
this curve—120 days—is exactly 12 times the 10-day 
period, this curve is not at all affected by this latter 
period. It is different, however, with the test for 
the 1o-day period. A range of 100 days only was 
used for this test, which is not a multiple of the 12-day 
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period. Since we have sufficient observations at our 
disposal, we re-calculate this test, using 12 ten-day 
periods instead of only 10. We then get the following 


figures : 
t. Sum. E: Sum. i Sum. 
fe) 75 4 310 8 187 
I 133 5 369 9 127 
2 I9I 6 310 
3 252 7 248 


On plotting these figures, we get the result shown 
in Fig. 25, which we recognise as a Class IV periodic 
function * whose equation is of the form 


400 


320 


NJ 


oe A ee a 
EVAR | 


240 


Fic. 25. 


ia 8C/ COS 23° eat cas hc 

=i Cos | 

4 x | “Ht 0 oe es 
On determining the amplitude, which is as near 

1§0 as possible, and dividing by 12, we find that the 

expression 


| eye) 27t 1 
12°55\=— 29 cos =— =| = (COS Oot 4 J cos Tot ae | 
nm” | ie 9 10 2 


* -, 


is a constituent of f(2). 


* Neglecting the gradually decaying function mixed up with it 
which is indicated by the ordinate for the roth day being less than 
the initial ordinate. 
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We next subtract the sum of this term and the 
12-day period term from f(t) to study the residue. 
The values of both terms are very easily read off from 
the graphs we have previously constructed of them. 
We do this for every fifth day, thinking that sufficient. 
The following results are obtained: 


ts Residue. t. Residue. 

fe} 10°3 65 5:0 

5 9:0 7° De) 
10 9:0 75 ee 
15 9:0 80 4°7 
20 8-7 85 4:0 
25 8-0 rele) 3°77 
30 Py 95 4:0 
35 7-0 100 4:0 
40 70 105 3:0 
45 6-0 110 33 
50 6:3 mas 3:0 
55 6:0 120 3°3 
60 fo) 


Since /(2) was only given in the first place to integers, 
these residues must be considered as uniformly 
decreasing. Plotting them reveals a curve of an 
exponential form which seems to decay to zero ultim- 
ately. ‘The time of falling to 1/e of its initial value, 
which is 10, is seen, from the graph, to be about 100 
days. 


Hence the expression 


: 100 2nt 
1o¢g7 °°! + Io sin (# +- =) —- Fo 5 — we (08 —— 
12 3 Mf 10 

++ : one ie eae ‘h 

10 J 


represents the given observations to well within the 
accuracy to which they are given: viz., half a unit. 


CHAPTER IX 
A BRIEF HISTORICAL SURVEY 


A.ttuoucu the purpose of this book is to impart 
such a knowledge of the subject that the reader shall 
be able to make practical use of his knowledge, a few 
historical facts about the development of the subject 
will probably be welcomed for their interest. ‘The 
history of our subject is very simple: * so much is due 
to Fourier and so little to other people. 

The fact that an analytic function could be repre- 
sented in what we now call a Fourier series was known 
in 1777 to Euler (1707-1783), who obtained the 
coefficients in the same manner as we obtained them 
on page 22. With an analytic function, of course, 
there is no question about the possibility of the 
expansion. 

Jean Baptiste Joseph Fourier (1768-1830) was, 
however, the first to assert that any finite artificial 
function could be represented, between any two finite 
limits, by such a series, and that the coefficients were 
given by the same formule as for analytic functions. 
Fourier’s assertions were at first emphatically denied 
by the leading French mathematicians of the day, 
including Lagrange, and his proofs were regarded as 
inconclusive. 

Fourier arrived at his results during the consider- 
ation of the solution of problems on the flow of heat ; 


* We are not now considering the complicated theoretical develop- 
ments of the theory of Fourier’s series which have been made in com- 
paratively recent years and with which this Practical Treatise is 
entirely unconcerned. 
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and his results, though published previously, are fully 
incorporated in his famous Théorie Analytique de La 
Chaleur (1822), one of the great classics of Mathe- 
matical Physics, which even now no mathematician 
ever reads without feeling a genuine admiration for 
Fourier’s genius, or ever refers to except in terms of 
praise. Fourier’s work and Fourier’s theorem will 
always remain fundamental in the solution of that 
great variety of important physical problems which 
can be formulated in terms of partial differential 
equations. No doubt Fourier’s proofs of his theorem 
leave a modern mathematician rather doubtfully 
convinced ; but Fourier’s work is a striking illustration 
of the fact that the intuitions of a genius are of 
infinitely more value to science than the proofs of a 
pedant. 

Fourier also gave the special application of his 
theorem to odd and to even periodic functions, and, 
further, the result now known as Fourier’s Integral 
‘Theorem, which is what the ordinary theorem becomes 
when the period is made infinite. In fact, in nearly 
every way, Fourier showed a much clearer compre- 
hension of the subject than many subsequent writers. 
Fourier alone, at first, grasped the fact that he had 
found a method of great importance. Applied to the 
analysis of analytic periodic functions only, the theorem 
was of no importance whatever: it was only an 
awkward way of effecting a transformation that could 
readily be effected by a little algebra, as in Chapter IV. 

Poisson (1781-1840) and Dirichlet (1805-1859) were 
among the first to give a proof of Fourier’s theorem 
by the direct summation of the terms of Fourier’s 
series, Dirichlet’s work being the more satisfactory, 
and rendering any further doubt about the validity 
of Fourier’s expansion, for functions which satisfied 
his simple conditions, impossible. 

Dirichlet proved that for any given value of ¢ the 
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sum of Fourier’s series was f(t), provided ¢ was not a 
point of discontinuity, and equal to 3{/(t— 0) + 
f(t + o)} at such a point. Dirichlet’s proof is repro- 
duced in a simplified and abbreviated form on pp.136—8. 
The fact that, at a point of discontinuity of a func- 
tion, the sum of the series was really indeterminate 
within certain definite limits, and that the range of 
such indeterminateness was 1-179 times the abrupt 
change in magnitude of the function, seems only to 
have been noticed comparatively recently ; though for 
ourselves we can hardly make out how it could have 
escaped detection in Fourier’s day that the series 


sin ¢ sin 2t 


oe So ee 
if Zz 


X sin x 
———dx can 
[e) 

take ; since everybody knew even then that the limit 


when ¢t = 0 of 


Af) + flat) + ft) +... + SO} 


when ¢= 0, can have any value that \ 


xX 
was | fodax, and one has only to multiply both 
oO 


numerator and denominator of each term of the series 
by ¢ to make the result instantly obvious. 

The subject of the harmonic analysis of a function 
defined by a set of points was not developed in any 
practical detail by these early mathematicians because 
the scientific observations requiring such treatment 
were not then extant. Lagrange (1736-1813), sub- 
sequent to his doubts over Fourier’s results, obtained 
the equation of a finite sine series passing through the 
tops of a number of equidistant ordinates; while 
Poisson extended this result to a finite series including 
both sines and cosines and so gave the fundamental 
formule for all harmonic analysis with equi-distant 
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ordinates. Lagrange’s and Poisson’s results are exactly 
the same as those given in our Chapter V. No doubt 
Poisson, had he required to undertake practical 
harmonic analysis, would have invented Schedules 
for calculating the coefficients almost indistinguishable 
from those we have given. Several formule for inter- 
polation for periodic functions, when the ordinates 
were not equi-distant, were given by Gauss (1777- 
1855). As far as we know, the complete list of these 
formule, with indications of the types of function 
to which they are respectively applicable, is now given 
for the first time at the end of Chapter V. 

One would think it impossible for mathematicians 
to pay much attention to the harmonic analysis of 
functions defined by a set of equi-distant ordinates 
without soon arriving at the expressions giving the 
coefficients of the different harmonics in terms of the 
actual Fourier constants of the function. It is some 
surprise to us, therefore, that we have not seen these 
relations in print earlier than A. Clayton’s article in 
the Fournal of the Institution of Electrical Engineers, 
Vol. LIX, 1920-21, p. 491; though it seems inevitable 
that several workers must have obtained these relations 
for themselves prior to this. 

Coming to the subject matter of Chapter VII, 
which is essentially due to Fourier, the late Lord 
Rayleigh (1842-1919) was the first to prove that if 


\ ft) cosatdt=A 
fats 
and fF) sinatd: = 8, 


Ce I (oa 
then \ fQ)edt = 4 (A2 + B2)da, 
——i0C ie) 


which is the extension to pulses, of the well-known 
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result that the energy of any periodic wave is the sum 
of the energy of its different harmonics. 

Finally, we must state that the very simple method 
given at the end of Chapter VIII for the practical 
determination of harmonic constituents in a non- 
periodic function is the same as that given in Whittaker 
and Robinson’s Calculus of Observations (1924). 
Although much attention has been given to this 
important practical problem, from the time of La- 
grange to the present day, the theory of the problem 
that we have given shows that nothing simpler can 
be obtained than this; while at the same time it shows 
that more elaborate methods are not worth the extra 
labour of applying, and are no more necessary than 
is the ideal diffraction grating, which only gives 
a single order spectrum, to the experimenter in 
spectroscopy. 
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